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ABSTRACT 


The electrical and thermal conductivities have been measured over a 
wide range of temperature of specimens of pure chromium, both in the 
ductile and embrittled states. It has been found that the low temperature 
conduction properties aré very similar to those of other transition metals 
investigated. At low temperatures the electrical resistance is mainly due 
to the scattering of electrons from the s-band to the d-band, while (s-s) 
scattering contributes about 60°, of the thermal resistance. 


§1. InTRODUCTION 


Wirt the discovery that very pure chromium can be prepared in a ductile 
form (Wain et al. 1954-55), it has become possible to obtain specimens of 
the metal suitable for the accurate measurement of its thermal and 
electrical conductivities. These measurements are an extension of work 
on transition metals previously reported (Kemp ef al. 1955, Kemp et al. 
1956). 

§2. SPECIMENS 

The specimens were supplied by Dr. H. L. Wain of the Aeronautical 
Research Laboratories of the Commonwealth Department of Supply in 
the form of rods, 8 cm long and 3 mm diameter. These had been prepared 
by forging high purity chromium in a sheath to form a bar of cross section 
11mm square, then grinding and electropolishing to size. <A full de- 
scription of the metallurgical properties of the metal are given in a series 
of papers by Greenway (1954-55). The purity of the metal is of the order 
of 99-998%.. 

The three specimens designated 1, 3 and 5 in the table, were supplied 
respectively in the cold worked condition, partially recrystallized and 
fully recrystallized. Recrystallization of ductile chromium leads to its 
embrittlement. After measurement specimens 1 and 3 were annealed in 
vacuo at 1050°c for 4 hours to become specimens 2 and 4. 


+ Communicated by the Authors. 
+ Now at the Division of Pure Physics, National Research Council, Ottawa. 
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The Specimens 
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Specimen Physical Po K,/T Kopo/L 
number State 10-7 ohm cm| W cm-! deg? | 10-§ W ohm deg 
it Cold worked 2-55 0-092, 2-35 
2 No. 1, annealed 1-81 0-135 2-44 
3 Partially 1-25 0-202 2-53 
recrystallized 
4. No. 3, annealed 0-90 0-275 2-48 
5 Fully 0-55 0-440 2-44 
recrystallized 
Theoretical Lorenz — — 2-45 
ratio 


§3. EXPERIMENTAL PROCEDURE AND RESULTS 
The thermal conductivity measurements were made in a cryostat 
described previously (White 1953), and the electrical conductivities were 
measured simultaneously with the aid of a galvanometer amplifier 
(MacDonald 1947). 
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Resistivity of chromium. 


Connections to the chromium were made by copper plating the ends of 
the specimens for a sufficient length to accommodate the end pieces and 
thermal and electrical connectors. These were attached by means of 
zine-cadmium solder, which is not superconducting down to 2°K, the 
lowest temperature at which measurements were made. 
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Thermal conductivity of chromium at low temperatures. 


Fig. 3 
[ n=2 
z 
oe Or 
eo 
v [ 
ef 
= L 
> - 
E 
2 
= = 
w 
a 
Ww 
a 
2 orl —o— SPECIMEN No! 
2 —w— SPECIMEN No.2 
? —e— SPECIMEN No4 
(= 
= —x— SPECIMEN No.5 
oO 
2 
= nl 
a ve 
4 
es 
oo; H4 Ne al a Ne ee) 
6 Te) 100 


TEMPERATURE (°K) 


Intrinsic thermal resistivity of chromium 


a79 


580 A. F. A. Harper et al. on the Thermal and Electrical 


The electrical resistivity of specimen 5 above 273°K is shown in fig. I. 
These results are similar to those quoted by Sully (1954, p. 92). 

The thermal conductivities and the intrinsic thermal resistivities of all 
the specimens are shown in figs. 2and 3. The intrinsic electrical resistivity 
of specimen 5 is shown in fig. 4. The high purity of specimen 5 is indicated 
by the low value (4:6 < 10-%) of the ratio, pp/p.273, of the electrical resistivity 
at lowest temperatures to its value at 0°c. 


Fig. 4 
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§4. Anaysis or Resuuts anp Discussion+ 


The electrical resistance is, to sufficient approximation, additively 
composed of the residual resistance p, and the intrinsic resistance pj due 
to the scattering of electrons by static imperfections and lattice waves 
respectively : 

P= Oc: Oi meet gese ck 0 id bien 


+ For a review of the theory see, for example, Klemens (1956). 
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po 18 independent of temperature and p; tends to zero at lowest tempera- 
tures. Thus p tends to a constant value at lowest temperatures. Values 
of py thus obtained are given in the table. The intrinsic resistance Pi, 
shown in fig. 4 for specimen 5, which had the lowest value of p,, varies as 


. [32 between 10 and 100°K. This indicates that most of the intrinsic 


resistance is due to interactions taking electrons from states in the s-band 
mto the d-band, for (ss) scattering leads to a 75 variation, while (s—d) 
scattering results in a 7 dependence (Kemp ef a. 1955). At about 30°K 
the intrinsic resistance is approximately given by 
ais 1G] ohimem. ae 12) 
If heat conduction by the lattice can be disregarded, as seems most 
probable in the present case, the thermal resistivity is similarly composed 
of a residual and an intrinsic resistance due to the scattering of electrons 
by imperfections and lattice waves respectively. 


WWW: . . c-. hae tes BD) 
where W, and py are related by the Wiedemann—Franz law 
esti gts ls 2b Shae tae tae Pee) 


and where L,=2:-45 «10-8 watt ohm deg~? is the Sommerfeld value of the 
Lorenz ratio. This relation should hold irrespective of the details of the 
electronic band structure. From theory W ; would be expected to vary 
as 7? for temperatures well below the Debye temperature (Klemens 1956) 
so that at lowest temperatures W; should be negligible relative to W, and 
K should vary as 7. This has been found to be so and in the table are 
given values of K,/7', the low temperature limit of K/7’, and of Kop,/T. 
The latter values do not differ significantly from L,, as is required by 
eqn. (4). 

Values of W; deduced from (3) for the various specimens have been 
plotted in fig. 3. Above 30°x all points fall on a single curve, which from 
30°K to 50°K, is a J? curve. Below 30°K there is a considerable scatter of 
the deduced values of W;. Below about 10°K W; is too small a fraction 
of the total thermal resistance to permit any conclusion to be drawn 
about the derived values of W;, but from 10° to 30°K there is a definite 
trend for these values to lie above the 7? curve except for specimen lL. 
Thus for all specimens 

W;=1-1x 10~-*7? cm deg W-! for 30°K<7'<50°K, . . (5) 


and for all except specimen no. | 

W,=4:3 x 10-4714 cm deg W-1 for 10°K<7T<30°K. . . (6) 
It is not clear whether the deviation from the 7” relationship below 30°K is 
significant. A careful evaluation of the experimental procedures has 
failed to reveal any source of error which could account for the deviation, 
If it is real it would not seem possible to account for it on the basis of the 
(s-d) scattering, especially as there is no corresponding anomaly in the 
electrical resistance. 
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If expression (5) is taken to apply to lowest temperatures 
PE 
In the case of transition metals the thermal resistance is due to (s—s) 


and (s-d) scattering, so that Wi=W(s, s)+W(s, d), similarly 


{Onis p(s, S) + p(s, d). 
Now at low temperatures 


W(s, s)ocT? and W(s, d)ocT”, while p(s, s)ocT® and p(s, d)ocT™ 


see Klemens (1956). For (s-d) scattering the Wiedemann—Franz law 
should be satisfied, so that p(s, d)=L,)7'W(s, d). Hence 


W(s, d) J 
Ws, d) +W6,8) 0 don ht qo 


and in the present case W(s,d)/W(s,s)=0-7. This may be compared 

with a value of 0-5 for the corresponding ratio in the case of palladium 

(Kemp et al. 1955) and 0-8 in the case of nickel and iron (Kemp ef al. 1956). 
The high and low temperature values of Wj are related by 

WiT,) @ 

: ASH’ ee oo! eee 

Wi(T's) T° eo 9) 


Li= 


From present value (6) of Wi(7,) and the high temperature limit of 
1-54 watt~1 cm deg for W(Z’,) (Sully 1954, p. 103) C=29, using for @ the 
value @,,=630°K, the Debye temperature derived from the low tempera- 
ture specific heat of chromium (Rayne and Kemp 1956). 

Theoretically, for a single band of N electrons per atom, C—64 N2° 
(Klemens 1954). If it is assumed that the ratio W(s, s)/W(s, d) is inde- 
pendent of temperature, this value should also apply to (9), because the 
single band theory should apply to the component W(s,s). Although 
there is gross disagreement (by a factor of order 4) between the theoretical 
and experimental values of C’ for the monovalent metals, in the case of 
palladium (Kemp ef a/. 1955) it was found that close agreement with 
theory was obtained if it was assumed that V—0-6 and a value of 1-506) 
was used for @ to allow for longitudinal waves only being able to interact 
with conduction electrons in an s-band satisfying the results of the Bloch 
theory. In the present cases if the same values of V and 6/6, are used 
the theoretical value is about 30°, less than the experimental one. When 
allowance is made for the uncertainties in the number of electrons in the 
s-band and in the appropriate values to use for @ and indeed for 6+, it 
may be concluded that, as with palladium, there is substantial agreement 


+ Other values for @p (Simon and Ruhemann 1927, Anderson 1937, Estermann 
et al. 1952, Wolcott 1956) are all lower than the value used here, the lowest 


being 418°K. A value of 538°K would give complete agreement between theory 
and experiment. : 


way Nie 
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between the theoretical’ and experimental results, in support of the 
conclusion that the Bloch theory yields better results when applied to 


(s-s) transitions in these metals than when applied to the monovalent 
metals. 
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On the Activation Energy of High Temperature Creep in Metals 


By Pavun FELTHAM 
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ABSTRACT 


The temperature dependence of the activation volume or © stress- 
concentration factor ’ q appearing in the expression 


de/dt=A exp (—-H/kT) sinh (qo/kT’) 


relating the equilibrium tensile creep rate de/dt, the temperature 7’ and 
the applied tensile stress c has been studied in the case of several pure 
metals subject to creep above about $7',,, where 7’,, is the melting point. 
On writing g=nb?, where b is the distance of closest approach of atoms in 
the lattice, the relation 


a= : Q i tt 
N=Ny EXP eal wie | 


was found to be valid in all cases, with V=«T',, and «=4-9 (cal/g atm)/°K, 
Le. Y is equal to about twice the latent heat of melting. The increase of 
q with temperature is ascribed to the increase in the mean spacing between 
adjacent logs in screw dislocations, and @ appears to be the activation 
energy for the conservative motion of jogs along dislocations. In face- 
centred cubic metals large values of n) are found to be associated with low 
stacking fault energies, and vice versa. 


§ 1. INTRODUCTION 


Ty a recent paper Mott (1956) suggested that in cases of high temperature 
creep in metals in which the steady, second-stage, or equilibrium creep 
rate de/dt obeys a law of the form 


de/di=A exp (—A/kT) sinh (qo/kT), 


where o is the applied tensile stress, H the activation energy of volume 
self-diffusion, and where A may be regarded as approximately constant 
over considerable ranges of the temperature, the rate determining 
mechanism may be the migration of segments of screw dislocations 
impeded by jogs. One would then expect the activation volume gq to be 
given by 

q=tlb?=in85, cee ee WD a OF | ele) 


+ Communicated by the Author. 
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where /(=nb) is the lengthof a segment of screw dislocation bordered on 
each side by a jog, and b is the Burgers vector. The jogs in the moving 
dislocations are assumed to form as the latter cut through other screw 
dislocation of the forest. Mott noted that 7 might depend upon the 
temperature ; as to the nature of this variability suitable data were not 
however available. The principal purpose of the present work was to 
obtain such data. 


§ 2. THe Errecr or MicrosrRUCTURE ON THE MAGNITUDE OF ¢ 

The creep data compiled by Kauzmann (1941) and the recent results 
obtained with aluminium by Dorn (1956) seem to show that in the steady 
state of creep g increases with temperature, but is independent of the 
applied stress. Seeger (1955 a,b), who considered the temperature 
dependence of the inter-jog spacing qualitatively for certain types of 
deformation, also believes that, irrespective of the distance travelled by a 
screw dislocation through the forest, an average distance between jogs, 
corresponding to a dynamical equilibrium, should tend to establish itself. 
Further, a preliminary compilation of g-values obtained from data on the 


creep of a number of metals (Feltham 1953, 1956), shows that 


(a) in pure polycrystalline metals g, measured in each case at the same 
fraction of the melting temperature 7’,,,, is of the same order of magnitude, 
e.g. with 7/T,,~0-6 g=70+10 (cal/g atm)/(kg/cm?) in lead, aluminium 
and zinc. 

(6) In annealed single crystals, in which essentially only one slip system 
operates during creep, g may be an order of magnitude greater than in 
polycrystals of the same metals flowing at the same respective temperatures. 

(c) Still higher g-values are found with single crystals of exceptionally 
high purity, e.g. g=2> 104 (cal/g atm)/(kg/em?) in the case of zine of 
purity 99-999°,,, flowing at room temperature, but 

(d) as the creep data on steels show (Feltham 1953), the effect of 
elements in dilute solid solution on the value of ¢ is much less pronounced 
in polycrystals ; ¢ decreases only slowly with increasing alloy content. 

(ec) In precipitation hardened dilute alloys q is about an order of 
magnitude smaller than in the corresponding pure polycrystalline parent 
metals. 

(f) Finally, ¢ does not appear to depend upon the grain size, provided 
of course the grains are sufficiently small in relation to the dimensions of 
the specimen to allow the latter to be regarded as truly polycrystalline. 
This conclusion is based on recent data on the creep of lead (Feltham 1956). 
Algo, the general law of the temperature dependence of q, indicated by the 
parallelism of the lines in the figure, does not seem to be affected by the 
temperature dependence of the grain size, for the parallelism was main- 
tained irrespective of whether the grain size was kept constant (as for 
copper), or was allowed to increase as the temperature increased (all 


remaining metals). 


586 Paul Feltham on the Activation Energy of 


§ 3. Tue TEMPERATURE DEPENDENCE OF q 


In view of the observations listed above, particularly however because 
of (a), it was decided to consider the temperature dependence of q¢ in poly- 
crystalline metals of comparable high purity (99-99%), thus attempting 
to avoid complications such as might arise with alloys, or with single 
crystals sensitive to handling or very small differences in impurity content. 

The results shown in the figure in which 4n (eqn. (1)) is plotted against 
T/T, where 7',, is the melting temperature of the respective metal, 
except in the case of the body-centred alpha-iron (actually 0-4%C steel), 
for which 7',, was taken to be the A,-point. (850°c). As the lines in the 
figure are drawn parallel, the equation of any one of them is given by 


; n=N exp (—Q/kT), NPA See ia 
with Q=«T',,, and «=4-9 (cal/g atm)/°K. For a given metal N is a 
constant. 
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The relation between n (eqn. (2)) and the ratio 7',,/7’, where 7’,, is the melting 
point, except for alpha-iron for which 7’,, has been taken as thé 
A,-point (850°c). The data have been compiled from the following 
sources: Kauzmann (1941) (alpha-iron); Dorn (1956) (aluminium) ; 
Tyte (1938, 1939) (lead and tin); Meakin (unpublished) (copper) ; 
Weertman and Shahinian (1956) (nickel) ; Feltham (1953) (gamma-iron). 


It should be noted that two separate lines are given for copper, depending 
on whether the applied stress o falls below or above a critical stress o,. 
In so far as the equilibrium creep must be regarded as a recovery process, 
it must be affected by variations in the mechanism of work hardening. 
Now, preliminary results obtained by Meakin (1957), and the study of 
the modes of work-hardening in polycrystalline copper (Feltham and 
Meakin 1957 a) and of the creep of lead (Feltham 1956) seem to point to 
the conclusion that the mode of hardening in the creep of copper at 
stresses below oc, is of the type treated theoretically by Friedel (1955), 
while above a, it is probably fine slip. The data relating to nickel were 
also obtained with o>o, (Meakin 1957), and this applies also to the 
remaining face-centred cubic metals. 
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— $4, Discussion 
The values of Q, of the latent heat of melting L, and of the ratio OL 
are given in the table. It can be seen that in the case of the face-centred 
cubic metals Q/L—2-0, i.e. the activation energy Q is small, just sufficient 
for “ melting ’ a group of about two atoms. The absolute magnitude of 
@, and the essential constancy of Q/L for all face-centred cubic metal for 
which data were available, suggest that the activation energy of the 
process described by eqn. (2) is that for the conservative motion of jogs 
along dislocations. The local derangement of atoms, being akin to melting, 
would however be expected to have an activation energy of the form 
Se LET Oy, edo eel be ea BB) 
involving the entropy of melting. It is probable, therefore, that in 
eqn. (2) , 
ean, OR OL Oe ea des o> onl) 
where, as is easily seen by substituting for V from eqn. (4) into eqn. (2), 
mg is the value of n when 7'=T7',,, and is obtained from the figure by 
extrapolation of the lines to the left. 


Sn|Pb| Al Cu Ni | Fe (gamma) 


2512-9] 4-6+| 66 | 8-5 8-8t 


+ Values uncertain. 


It can be seen that n, decreases in the order Al, Pb, Ni, Fe, Cu. Now, 
the stacking fault energies for aluminium, nickel and copper decrease in 
the order in which they are named (Seeger 1955a). On the basis of Seeger’s 
theoretical considerations one would also expect the stacking fault energies 
of the transition metals nickel and (gamma-) iron to be close together, 
similarly for the polyvalent metals aluminium and lead. Evidence in 
support of a high stacking fault energy in lead has also recently been 
obtained by Feltham and Meakin (1957 b) from experiments on the mode 
of work hardening of lead single crystals at low temperatures.” Thus i 
appears to decrease in the same order as the stacking fault energies. The 
full significance of the variability of n) from metal to metal is not, as yet, 
clear ; from Seeger’s (1955 b) treatment of the temperature dependence of 
the flow stress of close-packed metals n) would be expected to increase 
rather than decrease in the above order. 
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Small Angle Scattering from Deformed Metals+ 
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and R. D. LowpE 
Metallurgy Division, A.E.R.E., Harwell, Berkshire 


[Received January 29, 1957} 


RECENT experiments on the small angle scattering of x-rays of wavelength 
1-544 from pure metal foils have shown that the intensity measured 
in the range of scattering angles 1° to 6° is increased by cold working 
the specimen (Blin and Guinier 1953, Blin 1954). Blin (1954) has 
suggested that the effect is due to the formation in the metal of cavities 
with diameters of about 30 A, since the intensity of scattering was thought 
to be too great to be caused by dislocations (Dexter 1953). However, 
Beeman ef al. (1957) have shown that at least part of the scattering 
at low angles is due to rays which have suffered two successive Bragg 
reflections. Surface irregularities also appear to contribute to the 
intensity (Robinson and Smoluchowski 1956). These two effects may 
have been present in Blin’s experiments, and his interpretation is therefore 
open to doubt. 

Both these sources of error can be avoided by studying instead the 
low angle scattering of slow neutrons. Firstly, it is possible to use 
wavelengths sufficiently long so that no Bragg reflections occur ; secondly, 
for most metals the true absorption cross section for neutrons is very 
small, and therefore thick specimens can be used for which the surface 
effects are negligible. Such an experiment has been carried out on copper 
using a neutron beam (from the Harwell pile BEPO) which has passed 
through a bismuth filter (Egelstaff and Pease 1954). This filtered 
neutron beam had a mean wavelength of 8-74 and contained less than 
0-3°% of neutrons capable of Bragg reflection in copper. 

Polycrystalline copper specimens of commercial purity and 99-995, 
purity were used in the form of bars 0-6cm thick. Recrystallized 
material and specimens deformed by an extension of 10°, gave the same 
total intensity to within 0-003 barns steradian~' atom? over the range 
of scattering angles 2° to 12°. This scattering cross section is equivalent 
to less than one hundredth of the x-ray scattering observed by Blin. 
It follows that at least 99°, of the effect reported by Blin was due to 
processes other than true small angle scattering from defects in the volume 
of the metal. 

A full account of the experiment will be published later. 
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ABSTRACT 


Kairkwood’s variation method for the calculation of the polarization of 
an atom or molecule in an external potential V(r), uses a wave function 
) 


Yo{lted Very} 


where YW, is the wave function in the absence of the perturbation, and 
the summation is over the electrons. The energy is varied with respect to 
the parameter ». In this paper an extension of the method is presented 
using a wave function 


AES 


and varying the energy with respect to the function u(r). The optimum 
form is found to satisfy a second order differential equation which involves 
the single and pair electron density functions. 

This equation is solved numerically for argon, using the Hartree-Fock 
wave function. The form of wu turns out to be nearly quadratic, and shows 
that almost the whole of the induced dipole moment is due to the outer 
shell of the electrons. An estimate of the dipole-dipole energy is made 
using the function wu, which in the case of argon gives good agreement with 
that estimated by other means. 


$1. INTRODUCTION 


THe theory of the electronic polarization of atoms and molecules pro- 
ceeds by investigating the change in the energy and wave function under a 
perturbing electrostatic potential V(r), the nuclei being held in fixed 
positions. In the particular case of polarization by a uniform field E the 
potential has the form —E .r but other potentials may be of interest in 
dealing with problems of molecular interaction. It is possible to get a 
formal expression for the change in wave function due to V(r) by perturba- 
tion theory. This is not the most convenient method in practice, however, 
as it requires detailed knowledge of the excited electronic states of the 
unperturbed system (including continuum states). Most work has there- 
fore been based on variational methods. 
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a 


An important method was evolved by Kirkwood (1932) who used a 


variation function 
v= +n S Vir, onc es wp, ee ee IL) 


where Y, is the unperturbed wave function and the summation is over all 
the electrons. The parameter is varied to minimize the energy in the 
applied field. For an atom in a uniform field this leads to the expression 


a—(4iN){ fot (Sri 00865) % ark OS 


for the molecular polarizability, where N is the number of electrons in the 
atom. 

The chief disadvantage of a variation function of this type is its 
limited flexibility with only one variable parameter. Since Yo is multi- 
plied by a function of V only, no allowance is made for the different 
polarizability of different parts of the molecule, or for the screening of the 
inner part of the molecule from the external field by the polarization of 
the outer. Later work has attempted to correct this in various ways 
mostly based on changes in individual orbital functions (Buckingham 
1937, Sternheimer 1954). The aim of the present paper is to describe a 
method which is a simple generalization of the Kirkwood method and 
which gives direct information about the relative contribution of 
different parts of the molecule to the polarization. The variation function 


chosen is c 
F—= Woy 1+ Sur}, ee Boe EN 
where u(r) is a one-electron function. Since the modifying factor 
{1+ >,u(r;)} is symmetric in the coordinates of all electrons, the total 
function is still antisymmetric. The form of u(r) is then chosen to 
minimize the energy in the perturbed system. This leads to a simple 
partial differential equation. 


§ 2. GENERAL THEORY 


We shall use a wave function of the type (3) and minimize the energy 
for the Hamiltonian 


KH =H y+ > V(r,). . . : . . . : (4) 
It will be supposed that the perturbing potential V does not have the_ 
same symetry as %. 
Y} will satisfy the unperturbed equation 
Bp Gain . . . . . . . 2 (5) 
and will be taken as normalized. It will be assumed that no magnetic 


fields are present so that all wave functions may be taken as real. It is 
convenient to choose the zero of energy so that 


jr] = rea | Vt =0 4 t,o) 


the integration being over all electronic coordinates. 
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It can then be shown (Buckingham 1937) that the calculated second 
order change in & takes the form (using atomic units) 


62h | %S leraduiry Pyle +2 | bf | Suey | [s Vie, | ee ALUN 


This is the expression that has to be minimized by varying u(r). It can 
be simplified if we define two density functions 


px(r)=N | YW 2drdr,...dry ds, 


pls, r,)=N(N—1)| Wedradt,...dryds,ds,. . (8) 
the integration dr being over all but one or two of the N electrons. p,(r) 


is the electron density and p,(r,, rz) the density function in pair space. 
In terms of these functions, &, can be written 


Ey =} (grad u)?p,(r r)dr-+2o(r)u(r)de eat yC) 


where 


o(r)=pilr)V(r)+{V(r')palr, dr’. 2... (10) 


The integrals in (9) and ad are now over the coordinates of a single 
electron. If the function u(r) is varied in (9) 


36 .=|p, (grad wu . grad dw) dr +2 odu dr 
= | u{20—div (py ORO YE) Crees 2 at Ne. au uuh L} 


since p->0 sufficiently rapidly at infinity. Hence the energy &, is mini- 
mized if w satisfies the equation 
Cie; CCA ae oe ne oe ev os wae (EB) 
with appropriate boundary conditions. The value of &, is then given by 
either of the forms 
&,=—}| (grad U)*py dr=|ou Ge Ni Oe SP to relics) 


px(r) and o(r) are given functions if Y%) is known and eqn. (12) is then 
equivalent to the general electrostatic problem of determining the 
potential in a field of variable charge density and dielectric constant. 
In practice only an approximate function Wy will be available, but this 
can be used to find approximations for p, and c. 


§ 3. THe PoLarizaTIon oF ATOMS AND Ions 1n UNIFORM FIELDS 

The simplest application of eqn. (12) is to find the function u(r) for a 

spherically symmetric atom or atomic ion in a uniform field E. If 
this is measured in atomic units 

Pe COs weer! pe as os eee ee) 

using spherical polar coordinates with the direction of the electric field as 

polar axis. From the symmetry of the system o and w will also have the 

form 
o=FHs(r)cos0,u=Hg(r)cos@#. . . . . . (15) 


252 
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The Kirkwood approximation in this case corresponds to taking g(r) 
proportional to r. The equation for its optimum form is 


a 7) — “Ht =e. Spee eon eae) 


ee [ip pe =. 
redr\' lar 72 


The boundary conditions are that g=0 at r=0 since g cos@ must be 
uniquely defined, and that p,g?—+0 as r->oo so that the wave function is 
small at large distances. The polarizability « is then given by 


| Sorta. so P< oe ae ae 
0 


The simplest example is the hydrogen atom for which 
py! exp (—2r), saan y exp (= 2r): Poste See hes) 
The solution of (16) is then the familiar result 
g=—= (rb): os). et ee AS 
Equation (16) can be solved for other atoms and ions if estimates of 
the functions s(r) and p,(r) are obtained from Hartree-Fock wave 


functions. This has been carried out in detail for argon using the orbitals 
tabulated by Hartree and Hartree (1938). They give 


CONG) =arhrer? p(T )iae oo utc a ee Ce 
and P,,,(7) where P,,, is 7 times the radial part of the orbital, the angular 
part being normalized to (1/47). Defining 


PG OR oe Si Seon ee) 

(16) becomes 
aie dent d® 712 207-12 99 
Api ah peal ql \ ple ZUNE Ds eae) 


where ® is obtained from eqn. (10) and, after integrating over angle, 
becomes 
®=Ur—4 YY PyoPor | P..,Paridr. og t23) 
N=1,2,3 n’=2,8 “YO 
Equation (22) was solved numerically by taking F(0)=0 and adjusting 
#'(0) so that F(r)—0 as roo. A check on the solution is obtained by 
calculating « from both expressions (17). 

The resulting function g(r) is shown in fig. 1 together with the linear 
function g(r)=—0-6lr obtained by the Kirkwood procedure. The 
significance of the difference between these functions is best illustrated by 
the corresponding function M(r) defined so that M(r)dr is the contribution 
of electrons in the shell (7, 7+-dr) to the dipole moment in unit field. 

The total dipole moment is given by 


- CO 


U(r)dr=—|¥ Yr, cos 6, Pde sh fe eee ee Cee) 


~ 0 


= =H p11 cos bdr 
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where p,’ is the perturbed electron density, 
r 2 
pr (rs=N | Yo 1+ d9(";) cos a,\ dt,...dty 


=pilty) +24 {px(rs)g(r1) 008 0;+ | pal, t2)g("2) COs Oy dr'g}-+O(E?). 
: (25) 


In terms of the functions used in the calculation, the second term is 


= 


B oo 
= Ug= . 2 os Paola Gill | Pike wot) GOB U7 (26) 
J0 


77) n=1,2,3 n’=2,3 


Fig. 1 


A@ : 4°O 6-0 G.U 


(a) The function g(r) for argon. The continuous line shows the optimum form, 
and the broken line the best linear function, g(r) =—0-61r. (6) shows 
the behaviour of g(r). over the region of the ‘ inner shells ’. 
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Fig. 2 


(a) The function M(r) for argon. M(r)dr is the contribution of electr i 
shell (r, r+dr) to the dipole moment induced in unit field. reat 
tinuous line shows the function obtained with the optimum g(r), and the 
broken line that obtained with the linear g(r). (6) The behaviour if 
M(r) over the region of the ‘ inner shells ’. | . 
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so that ~ 
M(r=—3] 09-4 YS Pw Par 


R=1,2,38 n’=2,3 ( 


- 0 


PyP yt dr wee a) 
) 


M(r) is shown in fig. 2. 

From the general form of the functions it is clear that the Kirkwood 
method over-estimates the induced moment in the inner region of the 
atom but under-estimates it for the outer shell. The Kirkwood result 
is greatly improved if a variation function with two variable parameters, 


PATA or. 6) eh es es (28) 
is used. In fact the function —0-02r—0-49r?, obtained in this way, is a 


very close approximation to the optimum form. The calculated values 
of « are compared with the experimental value in the table. 


The polarizability of argon 


Optimum g(7’) 2-03>el0n2— em= 
g(r) = —0-61r 0:98 x 10-74 em? 
g(r) = —0-027 —0-497? 1-95>..10-"* cm? 
Observed 1-69 x 10-24 em’ 


$4. THe DirroLe—DreoLe INTERACTION OF Two ATOMS 


An estimate of the dipole-dipole term in the van der Waals interaction 
between two atoms at large separation may be made using the function 
g(r) calculated in § 3. 

Let R be the internuclear distance and (7,, y;,2;), (x,', y;',2,;) the 
coordinates of the electrons in the two atoms relative to their respective 
nuclei, with the z-axis along the line of centres. Then the dipole— 
dipole interaction is represented by the addition of a perturbing term to 
the Hamiltonian 

V=R* > = (x ay’ +y ys —222;) . . . « (29) 
i=1 j= 
where NV, and JN, are the numbers of electrons in each atom. 

If RF is sufficiently large for overlap between electrons on different 

atoms to be negligible, the wave function for the unperturbed system is 


La ye Vie . . . . ° ° . . ° (30) 
where ¥, and Y, are the ground state wave functions of the two atoms. 
Writing the perturbed wave function 


Y=YPo{l+u(r,, r;')} Pw an te 5 al) 
the second order energy takes the form 
N, N, 
é=1[¥,¥ 2 (grad w)?+ > (grad yu)® Pail 
i=1 j=1 


+2 ¥,PouV P,P ydr dry. Saka tee ane eS) 
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The variational function suggested is 


N, Ne Hide pf 1 ' ’ : 
PWS SS BOE +Yy ys — 22,2) . - - (33) 
j=1 j=1 i! j 


A being chosen to minimize the energy. Physically, this means that the 
field at atom 1 due to each configuration of electron in atom 2 is weighted 
as if atom 2 is polarized in some average field of magnitude AR, and 
vice Versa. 
With the function (33), 
6 = RB {3V?(4 By +083) + 3Ax 4% } oo heron ye ee (Oa 
where «,, %» are the polarizabilities of the atoms, and 


p=—| Mingiryidr, 2 8) 
0 
The minimum value of &, occurs when A=—«%5/(%;B.+a8,) and gives 
for the dipole-dipole constant 
By toe 
6 ee ee 
te? Bot ooP 
When g(r)=Ar, then this gives Kirkwood’s formula (1932) 
L425 ee 
C 2 38 . 37 
2B (04|N 1) V/ (Ha/ 2) ae 
For atomic hydrogen, where g(7)=—r—tr?, (36) gives c=6-36 a.u., 


compared with 6-00 a.u. given by (37). A more accurate variation treat- 
ment gives c=6-499 a.u. (Pauling and Beach 1935). 

For argon, using the optimum g(r) and the calculated polarizability 
the result is c=109 x 10~-® ergs/em®. This is expected to be larger than 
the true value since the polarizability is too large. In fact this value is 
close to the inverse sixth power coefficient in the Lennard—Jones 6-12 
potential obtained from viscosity data, which contains contributions from 
induced quadruple and higher moments. 


ACKNOWLEDGMENT 


One of the authors (P.S.) wishes to acknowledge the receipt of a 
D.S.I.R. Maintenance Allowance during the tenure of which this work 
was carried out. 

REFERENCES 
BuckincuaM, R. A., 1937, Proc. Roy. Soc. A, 160, 94. 
Hartree, D. R., and Hartren, W., 1938, Proc. Roy. Soc. A, 166, 450. 
Kirkwoop, J. G., 1932, Z. Phys., 38, 57. 
Pavutine, L., and Bracu, J. Y., 1935, Phys. Rev., 47, 686. 
STERNHEIMER, R. M., 1954, Phys. Rev., 96, 951. 


The Inter-Diffusion of Two Charged Particles (with particular 
reference to Ion-Exchange in Zeolites)+ 


By I. R. Bratrie 
Department of Chemistry, King’s College, London 


and 


D. R. Davizs 
Department of Mathematics, The University, Sheffield 


[Received January 2, 1957] 


SUMMARY 


The general importance of ion exchange is mentioned, together with the 
work of previous investigators. The mechanism whereby two ions may 
change places is examined, and it is concluded that in many cases two ions 
diffusing in opposite directions in a charged substrate will diffuse inde- 
pendently of one another apart from any space charge effects. An equation 
is developed relating the amount of material exchanged during a reaction 
with the square root of the duration of the experiment. Practical 
applications of this equation are briefly discussed, as is the effect of 
immiscibility of solid solutions within the exchanger. 


§ 1. IyTRODUCTION 
lon exchange may be represented by the simple equation 
ZA (s)+ B+ (aq.) = ZB (s)+ A? (aq.), 

where A and B are, for example, two unlike monovalent cations. The 
anionic framework (Z) of the solid exchanger (ZA) may be either organic or 
inorganic in nature. Both the kinetics of the exchange process and the final 
equilibrium position are of fundamental interest, although only the 
former will be discussed in this paper. In the case of zeolitic exchangers 
(aluminosilicates containing loosely held water), the material is usually 
available only as an irregularly shaped powder. Studies on the diffusion 
of ions in to and out of the zeolite are frequently carried out by im- 
mersing the solid in a well-stirred solution of known concentration. The 
mathematical treatment of the results is considerably simplified if the 
time is short enough, or the diffusion coefficient low enough, for the 
external surface area of the solid to be treated as a reasonably accurate 
representation of the cross sectional area into which diffusion has occurred. 
In this case the powdered zeolite may be regarded as a semi-infinite solid. 
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A further simplification is achieved if the concentration of the exchanging 
species in solution can be kept at a constant value throughout the duration 
of an experiment. The concentration at the interface may then be put 
equal to the equilibrium concentration at time infinity. The usual 
method of maintaining a constant concentration is applied by having two 
solid salts, each containing one of the exchanging species, in equilibrium 
with the solution. 

Quantitative data on ion exchange in zeolites have been obtained by 
Barrer and others (1953, 1954, 1955, 1956), although in many instances 
the presence of two immiscible solid solutions has complicated any 
theoretical treatment of the results. Conduction processes in anhydrous 
analcites and chabazites (M,O. Al,O,.4SiO,) have been examined by 
Beattie (1954). In certain cases these zeolites constitute an ideal system 
where the alumino-silicate network remains unchanged throughout an 
ion exchange. Further, the possibility of cation—cation interactions 
other than through electrostatic forces is unlikely in a zeolite such as 
analcite. Each ion may then be considered to be uninfluenced by the 
nature of the neighbouring cations, apart from any stoichiometric require- 
ments. A general understanding of ion exchange is essential for many 
purposes. Examples of ion exchange are to be found in _ biological 
systems, in soils, in water treatment, in the weathering of glasses, and in 
the separation of the rare earths. Fundamental studies on ideal ex- 
changes are likely to be of particular importance in the interpretation of 
more complex cases. A formal mathematical treatment of ion exchange 
has not previously been given in the literature. It is the purpose of this 
paper to examine mathematically the process of ion exchange in a system 
in which the normal Fickian diffusion would occur in the absence of any 
effects due to an induced electric field. 

There is a number of ways in which ion exchange may occur in crystals 


in general, and in zeolites in particular. These are considered separately 
below. 


§ 2. Direct ExXcHANGE 
This is essentially the simultaneous exchange of places by two ions A 
and B. It may be brought about in a number of ways, involving different 
numbers of ions. Any mechanism involving only two ions is unlikely for 
steric and energetic reasons. It may be noted that none of these processes 
is significant in relation to the conduction of an electric current. The 
exact significance of ion exchange rates from the point of view of individual 
diffusion coefficients is difficult to assess if simultaneous migration of two 

different species is a fundamental of the process. 


§ 3. INDEPENDENT MIGRATION 


A process involving the simultaneous migration of several ions is likely 
to be less favourable than diffusion by a defect mechanism. It has 
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recently been shown (Beattie 1955) that over part of the composition 
range both ions in a heteroionic crystal of analcite may migrate inde- 
pendently, apart from a certain necessary stoichiometry. In a case 
where ion exchange is occurring the two ions will diffuse independently 
apart from any space-charge effects. If two similarly charged particles are 
migrating in a structure, unless the diffusion coefficients of the two ions 
are equal, an electric field will be developed in such a direction that the 
rate of migration of the slower moving ion would be expected to increase, 
while that of the faster moving would decrease. In the steady state the 
two rates would be equal, with the measured rate of ion exchange de- 
pendent upon the diffusion characteristics of the two ions concerned. 

Mott and Guerney (1950) give for the migration of an ion in an electric 
field the formula 
0c, 
On 
where J,=ionic current, D,=diffusion coefficient, C,=concentration, 
H=field strength, v;=ionic mobility and x= distance in the direction of 
diffusion. 

Suffix 1 refers to ion 1. Where two ions are migrating in opposite 
directions an exactly analogous expression may be written for ion 2. 
The elimination of H between these two equations, plus the use of the 
Einstein relationship between diffusion and conduction, yields eqn. (2). 

of -- he = Stn (GAGE) |) Samat Sere) Meets) 

The ionic current at w is J,, while at (w+6z) it is J,;+(0J,/dxv) dx. The 
rate of accumulation between the two planes defined by x and (x+62) is 
therefore given by 


Dee De BUG ee ee Scot) 


eetod Hoa yer—= (COW Cilor. - 2.) se ae 4S) 
In order to proceed with this analysis it is necessary to make a simpli- 
fying assumption that the ionic flow in each direction is the same. This 
is equivalent to assuming the attainment of steady state conditions and 
can only be regarded as an approximation. 
Putting J,=—J,, and hence 0C,/d¢x=—0C,/dx, eqn. (2) becomes 
jee 0c, C*D,D, | (4) 
; Gl oly C*D,+C,(D,—D,) ‘ j ; ; 
with a corresponding expression for J,. Here C*=C,+C,= a constant. 
Using eqn. (3) we obtain the result 
1 ay _ af a ee ee ea 
For the case of the semi-infinite solid the boundary conditions are 
given below. 
C.=C,. att=0-tor «>0. 
C=C, at «=0 for t>0. 
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We first derive an approximate solution which is applicable in certain 
cases. Putting C*D,Dt=T, and 
C*D,+C;(D:—D,) 
C*D2+C (Di—Ds) 
then eqn. (5) may be written 


00, exp (—6,) 070, 


= exp (9) 


I~ CFD +p (D:—Dy) aa i 
The new boundary conditions are 
Oi = 0 at 70 tora 0: 
evan) email = 0, abiae==0 fOr te 0, 
If 6, is small, then this equation may be simplified to read 
OO ee 
OT = C*D,+C)(D\—D,.) ox? 


Equation (7) will be valid for the conditions 

(a) T' +0, when 6, + 0, 

(6) D, is approximately equal to D, and 

(c)eCr is nearo Cy. 

Although conditions (6) and (c) are not met with frequently, the 
solution of eqn. (7) may be of use in some particular cases. The appro- 
priate solution of eqn. (7) may be expressed in the well known and easily 
computable form 

6,=6,,/[1— erf {a./[C*D,+C,'(D,—D,)2V/T}). . . (8) 

It is of interest to note that Douglas and Isard (1949) approached this 
problem in a rather different manner. If J, were a function of a and C, 
only, then it would be possible to-integrate eqn. (4) under suitable con- 
ditions. The relationship assumed by these authors is essentially 
C=. exp (—ba)-La, where a and 6 are constants. This assumption 
enables J, to be put equal to some constant value over the range x=0 to 
x=d. Proceeding in this manner these authors obtain a solution appli- 
cable when C,,’=0. 

In this case 


ayer 2D,D, dD, D, ) 
ono a Bylo p (a) If]. 


where Q is the total flux that has crossed «=O after time f. 

An exact solution for all values of D, and D, can, however, be derived 
from the analyses of Fujita (1952 a,b) given in another context. Fujita 
has recently obtained a formal solution of the equation 


OG. 0 oC 


D=DiU= 0) Ea en 


where 
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D, and 4 are both constants. Writing 


pee a eh ; 
LEE (025 SETI 8 tm 
baer Di) te 

eee D)0pl te pio okat! Se 


and C=C,—C,’, eqn (10) becomes identical with eqn. (5), and the 
boundary conditions are now expressed in the form 

C==0, +>0, t=0 

eC gst, =U east 0, 
The solution given by Fujita (1952 a, b) is valid for 0<«<1, where «=AC). 
This analysis is discussed fully by Crank (1956), the final solution being 
written as two parametric equations : 


1 -0 
aoe ory | In 6?) 26} exp | | a a ao, | (14) 


“S 


it q i 
= acetone, (0,2—p In @,?)-¥ i,t | (15) 


fl 
p=—In (1—2)=2{ (6,2—w1n8,?)-12d8,, . . . (16) 


~ 0 

C/Cy=c and y=a2/24/(D,t). Now, 
PO PCs Oy) 

[C*D2+(Di—D2)Co] 
where all the physical quantities are positive. Since C,,’ must of necessity 
be less than C*, it can be shown that « can only take values between +-1 
and —o. A solution for — 0<«%<1 is therefore required for this paper. 
We find that the analysis by Fujita may be extended, with some modi- 
fication, to cover the additional range «<0. The final solution for this 
case may be rewritten 


(6+ (62—p In 62)12} 


and 


where 


ry 


i 
ee ven 
9 
xexp — {{ (G° = in 6 5°)-3" i,b—2 exp (—#/2) | sy ttle) 
0 
(so that y > 0 as 6+ 1 and y— + © asd 0), 


and 
ees tex (0,2—p In 6,2)? do ja] eg as 
“seq p42], Capa : 


where in this case 
ib 
p=+In (1—2)=2| (6,2—p In6,2)-H2d0,, =... (19) 
0 


The relationship between « and p in eqn. (16) is given in graphical form 
in the paper by Fujita (1952a, b), and also in Crank (1956), and the 


604. I. B. Beattie and D. R. Davies on the 


corresponding relationship in eqn. (19) may readily be deduced from this. 
6 

Only one numerical integration (of { (0,2—p In 6,?)-1? ds) is required, 
0 


since this expression occurs in all the final equations. It is then possible 
to write down values of y and c for appropriate values of @, giving a 
parametric solution of the original equation, for —co<«<l. 

It is frequently convenient to study ion exchange from the point of view 
of the total quantity of material exchanged during a reaction. For 
example, the change in weight of the sample may be followed. In such 
cases it is necessary to derive an expression for the total flux which has 
passed per unit area of the plane x=0 up to a time f. 

In both the above instances («<0 and 0<a<1) the required solution 


hes ( —C*D,D, (32) dt 
o LC*De FC, (Di—Do)] \ 0% Jose 


where Q is the total quantity of material passing through the plane x=0 


up to time ¢. 
/de _ {de. dé. dy 
de) pg \dB dy de) 9” 


Now 
Using eqns. (17) and (18) 
(=) 24? exp (B/2)62 
«=0 


oF Sees it et 20) 


~ [exp (—f)D I? 
while (14) and (15) yield 


is 


G __ —21/2 exp (—f/2)t-12 
Ti) yy Tex (B)—UD a? a 
The final expression for Q with «<1 is then given by 

Q oe 93/24 (C*D,D,)*7t1/2[0* D+ (D,—D,)C.’]¥? (22) 


p*/?(D,—D,) : 
where A—surface area of the interface. 
In the limit as D, > D, (i.e. « + 0, 8 + 0, and » > o) it may be shown 
that : 


p®(D,—D,)— +/(27). 
so that eqn. (22) then reads 
Q=2A(C,,’—C,’) vag) 2 ay SS eee 


where D=D,=D,. As expected, the expression is that found for simple 
Fickian diffusion in a semi-infinite solid. 
In certain simplified systems ion exchange across a membrane is of 
interest in the steady state. Equation (4) may be written 
7, ae (ee ) aC, 
\O*D,+C,(D,—D,)]° dx 


(24) 
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for the case of the steady state flow between the planes «=0 and a—d. 
Integration of eqn. (24) then yields the expression 


‘ (Di—D,)d°- \C*D,+C"s(D,—D,) 


since J, is a constant throughout the membrane. C’, and C’, are 
concentrations of C',, at x=0 and x=d respectively. 


§ 4. Dirrusion wirH LvMIscIBILITY PRESENT IN THE SYSTEM 


There are two main possibilities in the case of granular material, either 
diffusion of one ion is rate controlling in the outer layer, with diffusion of 
both ions independently present in the inner region, or conversely the 
migration in the outer layer is dependent on the diffusion of both 
ions, while diffusion in the inner region is rate controlled by one 
ion. Only the former case is at all amenable to calculation. There 
are effectively two fixed composition limits for the outer layer, the com- 
position at the interface of the two solid solutions, and also the com- 
position at the surface with the aqueous solution. Since the diffusion in 
this layer is controlled by one ion only, no electrical field will develop. It 
appears, therefore, that this is analogous to a tarnishing reaction and the 
equation of Danckwerts (1950) may be applied. 


§ 5. APPLICATION 


The major contribution in this paper is the formulation of an equation 
relating the quantity of material exchanged during an ion exchange 
reaction with the individual diffusion coefficients of the two migrating 
species. In the case of Fickian diffusion the solution obtained is 


_ 284A (C*D,D,)*07[0*D y+ (Dy—Dy)Cn'}!” (25) 
er oe a a 


The most useful practical applications of the equation are likely to be 
when the diffusion coefficients differ by orders of magnitude, and where 
@, —0 or: Ct. 

When D,>D,, eqn. (25) may be simplified to read 


23/24 [0*(C*—C,,')]2 
= mamas eo 


In the case where C,’=0, « is found to have the value C,,'/C*. The 
observed diffusion coefficient is then equal to that of the slower moving 
ion multiplied by a constant dependent only on concentration. Where 
Cy'=C*, « has the value D,(C.,’—C*)/C*D,, dependent on both diffusion 
coefficients. 

Where experimental results yield only values of @ as a function of v7, 
it will in general be impossible to obtain individual diffusion coefficients. 


V(Djft). 
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SUMMARY 


Single crystals of magnesium oxide have been subjected to ultra-violet 
x and neutron irradiation. The bands produced by these radiations, 
their optical and thermal bleaching and associated luminescence are 
discussed. Taking into account the data of other workers, schemes for 
ultra-violet activation and luminescence are proposed. The ultra-violet 
activation is attributed to the freeing of positive holes, perhaps with 
exciton formation as an intermediary process, and the subsequent 
trapping of both electrons and positive holes at negative and positive 
ion vacancies. The role of the various impurities is discussed and it is 
concluded that they do not directly give rise to the optical absorption 
bands observed. The luminescence scheme assumes an optical band 
width of over 7-0ev. It is then shown that the various luminescence 
photon energies correspond to values which are to be expected on the 
assumption that the optical absorption energies represent the heights of 
bands measured from the top of the valence band. 


§ 1. INTRODUCTION 


HIBBEN (1937) was the first to report that periclase (single crystals of 
magnesium oxide) could be coloured by light from the 2537 A mercury 
line. Subsequent to this Molnar and Hartmann (1950) gave details of 
certain absorption bands developed under x-irradiation and Weber (1951) 
showed that these bands could also be produced by additive coloration 
with the parent atoms. Day (1953) examined photoconductivity in 
neutron-irradiated crystals, and details of the luminescence and x-ray 
induced photoconductivity were given by Einstein (1954). 

The present series of experiments fills in some of the gaps. Details of 
the ultra-violet induced activation are given and some of the effects of 
neutron irradiation examined. These results together with those of 
previous workers are then discussed in the light of an elementary model. 


§ 2. EXPERIMENTAL DETAILS 


The magnesium oxide single crystals were obtained from the Infra-Red 
Development Company of Welwyn Garden City. Due to differences in 


+ Communicated by Dr. H. M. Finniston. 
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the optical absorption property between various specimens, a detailed 
spectrographic analysis was carried out to see whether these differences 
could be associated with differing impurity content. Typical analyses 
are shown in table 1. The specimens A and B show the order of variation 
typical of the crystals, and specimen C is of one kindly provided by Dr. 
B. Bleaney of the Clarendon Laboratory ; these figures are accurate to 
within a factor of two either way. An analysis and examination of about 
30 different crystals show that no apparent correlation existed. 


Table 1 


Impurity (parts per million) 


Ni | Be| Fe |Mnj| Al | Cr Ca Co Si 
Al 10) 107 100s 10 4100 10 | 500 10 10 
Bar LOOT SLOOsebO SLOG TN e107) 200 ee LOOR Aor 
C |.10 | 10 | 700 | 10 | 100 | 100 | 1500 | — | 700 


Thermal annealing of the crystals at 1200°c in vacuo did not make the 
absorption spectra more consistent and so the following method of 
standardizing the crystals was adopted. If the intensity of the light 
incident on one face and emerging from the opposite face of a crystal of 
thickness L be represented by J, and J then 


Inia plein KI, <9. 2 ee 


where p is the absorption coefficient and K is a constant related to the 
refractive index. Hence a plot of crystal thickness versus log trans- 
mission for any given wavelength should give a straight line of slope — 
and an intercept on the transmission axis of In KJ,. It was found that 
only certain of the crystals followed the straight line reasonably and gave 
a correct intercept of the transmission axis ; this intercept was calculated 
using the values of refractive index given by Strong and Brice (1935). 
The crystals which satisfied the relation were selected for the experiments. 

Some of the crystals used were cleaved pieces, but where special thick- 
nesses were required it was found necessary to grind the specimens on 
diamons impregnated lead and nylon laps, special precautions being taken 
to remove the lubricating layer subsequently. 

The neutron irradiations were carried out in the Harwell BEPO pile, 
and at Windscale where higher fluxes were needed. In BEPO three 
separate facilities were used : 

(a) Where the ambient temperature was 60°c—80°c and the fast-to-slow 
neutron flux ratio was 1: 1. 

(6) In the centre of a hollow uranium slug where the pile temperature 
was 20°c and the fast-to-slow flux ratio was 3: 1. 


Single Crystals of Magnesium Oxide 609 


(c) For irradiation times from 1 sec to 1 hour, in a pheumatic tube 
facility where the fast-to-slow flux ratio was 1: 1, and the temperature 
50°c-60°c. In this facility for short time irradiations the specimen does 
not have time to attain the pile ambient temperature, and so all specimens 
were raised to slightly above pile ambient before irradiation. 

Absorption measurements in the range 6 ev (2060 A) to 1-2 ev (1030 A) 
were made on a Hilger ‘ Uvispek ’ spectrophotometer. Measurements in 
the infra-red were carried out by Mr. R. Wood of A.W.R.E. Aldermaston 
on a Unicam double-beam recording spectrophotometer, and in the range 
5 ev (2470 A) to 7-5 ev (1640 A) on a vacuum spectrophotometer details of 
which will be given elsewhere. The thermal annealings were made in a 
vacuum furnace capable of being maintained within a temperature of 
+2°c. Light intensities for the bleaching experiments were measured 
with a Schwartz thermopile, and where necessary the light was mono- 
chromated by passing it through a Hilger D222 quartz monochromator. 
Light at wavelengths of 2537 A and 3650 A were obtained from special 
lamps giving predominant lines at these wavelengths. Where lumin- 
escence photon energies are quoted, it means that the brightness was 
sufficient to allow inspection on a Beck pocket spectrometer. 

The optical absorption figures given on graphs refer to the reflection 
corrected absorption. The actual transmission measurement made on 
the spectrophotometer is //I, where J) corresponds to 100°, transmission 
and I/T)=K exp (—pL) 
where K=16n?(n+1)+4 
is the reflection correction, n being the refractive index at the wavelength 
of measurement. If after irradiation the absorption coefficient changes 
from yu to »’, corresponding to a transmitted intensity change from J to I’, 
then —(p’—p)= {In (I/D) }/L. . 

This quantity has been called the * Irradiation Induced Absorption 
Coefficient ”’. The number of defect centres contributing to any optical 
absorption band has been calculated from the “ Irradiation Induced 
Absorption” versus “‘ Photon Energy” curves using Dexter’s (1956) 
formula 

N j=0°821 x 1017n(n?-+2)-2f-? f u(H) dz em 

where NV, is equal to the number of centres per cm*, # is equal to the 
photon energy in ey and f is the band oscillator strength. This has been 
taken as unity, so that Ny will represent a minimum. It has been 
assumed that the refractive index, n, is not appreciably changed by 
irradiation. It was checked that irradiation did not alter the surface 
seattering properties of the crystal to +1%. 


§ 3. Puriry or CRYSTALS 


As already stated crystals from different fusions showed variations in 
impurity content. These crystals always showed initial absorption bands 
with maxima at 4:3 ev (2880 A), 5-6 ev (2200 A) and 6-8 ev (1820 A). The 


Zn 2 
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heights of these bands varied slightly between the various erystals, but 
only the 6-7 ev band seemed to vary as the impurity content, becoming 
greater with larger amounts of impurity. This band was found to have 

an absorption maximum between 35 cm~! and 25 cm~! in the various 
crystals tested. The other bands which varied slightly between the 
crystals could be removed by heating at 1300°c in a vacuum better than 
10-> mm of mercury for 20 hours, but the 6-7 ev band still remained (fig. 1). 


Fig. 1 


5-6 ev. 
(84x10) 


ABSORPTION COEFFICIENT (cm™') 
vn 
*) 
Cee ee ae a ea one oe ) ere | a | 
_— . 4 


PHOTON ENERGY (€v) 


The optical absorption of an untreated crystal, with impurities 
corresponding to analysis A of table 1. 
- —— -: The optical absorption of a crystal heated in vacuo at 1300°c for 
20 hours, and allowed to cool to 100°c in 4 hours. 
— — — The suggested bands responsible for the profile of the untreated 
crystal. 


The effect of irradiations on crystals with varying impurity content 
will be mentioned later, but it can be said here that the amount of 
impurity did affect the maximum values of ultra-violet (2537 A) 
irradiation-induced absorption. In the crystals with the larger amounts 
of impurity the maximum attainable absorption was slightly lower. No 
compensating absorption band was observed anywhere in the region from 
the absorption edge to 0-2 ev (60 000 A) where transmission begins to 
drop rapidly. For x-ray induced absorption however there seemed to 
be no such differences, the maximum absorption being the same within 
the limits of experimental error, 

Using paramagnetic resonance techniques Bleaney and Orton (private 
communication) have established that iron, chromium and manganese 
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ions exist as Mn?*, Cr’+ and Fe?+ or Fe?+. In crystals which do not 
contain Fe*+ initially, ultra-violet irradiation produces Fe*+ ions which 
are then fairly stable, annealing on heating at around 600°c. In erystals 
which do contain Fe** initially the amount is either unchanged or de- 
creased by irradiation. The Cr?* ion is always present, and its amount 
may be changed either way by ultra-violet irradiation. These experi- 
ments make it probable that the impurities exist substitutionally but 
the charge compensation which would be expected in the case of the Fe?+ 
ion is not contained in the ions immediately surrounding the impurity. 
The crystals having the highest impurity contents did not fit the straight 
line of eqn. (1) and so were not used for the other irradiation experiments. 


§ 4. EXPERIMENTAL RESULTS 


4.1. Ultra-Violet Irradiation 

frradiation with light of 25374 (4:95 ev) wavelength produces bands. 
with definite maxima at 4-3 ev (2880 A), 4:95 ev (2500 A) and 5-5 ev (2250 A). 
In addition a shoulder appearing on the low energy side of the 4:3 ev band, 
may be shown to be due to a band having a maximum at 3-7 ev (3340 A), 
and a broad band having a maximum somewhere between 2-4 and 2 ev 
may be broken down into two bands having maxima at 2-35 (5250 A) and 
2-05 ev (6000 A). This irradiation-induced optical absorption due to 
ultra-violet light is shown in fig. 2, where the resultant curve is broken 
down into the various suggested bands. The minimum number of 
centres responsible for each band is also shown ; this has been calculated 


Fig. 2 
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The full profile represents the irradiation induced optical absorption after 
exposure of the crystal to light of 2537 A wavelength. This profile is 
broken down into suggested components, and these bands are shown by 


the lighter lines. 
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from Dexter’s formula using an oscillator strength of unity. This parti- 
cular curve does not represent fully developed absorption. The 4-95 ev 
band decays away in seconds at room temperature so that it does not 
show up clearly on the absorption profile. At around —50°c this decay 
is much less rapid, and it is possible to make a definite measurement of 
absorption. The other bands are unaffected by this lowering in tempera- 
ture, and show no corresponding changes. This observation applies only 
to the visible and near ultra-violet regions of the spectrum and low tem- 
perature measurements in the infra-red and vacuum ultra-violet have not 
been made. The growth of optical absorption with time of ultra-violet 
irradiation is shown in fig. 3. It has been observed in some cases how- 
ever that optical absorption has decreased after some minutes of irradiation 
and then it has increased again. Some thirty separate growth curves 
have been plotted and of these six showed abnormal variations under 
carefully controlled conditions. However fig. 3 shows the majority 
results, which agreed between different crystals to within +5%. 


Fig. 3 
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Growth of absorption coefficient at various photon energies, with time of 
irradiation with 10~* watts cm-* of light of 2537 4 wavelength. 


These bands may be produced by light of smaller photon energies, but 
the efficiency in producing optical absorption decreases with the energy. 
Day (1953) has produced ultra-violet activation of MgO in respect of 
photoconduction at 2-1 ev, by light of as low an energy as 3-9 ev (3120 A). 
This activation could not be reproduced in the present series of experi- 
ments for crystal which fitted the straight lines of eqn. (1). On crystals 
other than these, however, light of energy down to 3-5 ev (3500 a) has 
produced appreciable absorption, and activation in respect of photo- 
conductivity. Light of energies up to 7 ev (1760 A) has been used for 
producing optical absorption, and though only a small amount of light 
was available at these energies, the quantum efficiency was very high. 
At several photon energies below 6 ev, the final saturation of aibcot ption 
was found to depend only upon the total amount of irradiation incident 


on a given area, and was independent of dose rate. This also applies to 
photoconductivity (Day loc. cit.). 
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Hibben (1937) reported the complete decay in the dark of the visible 
bands. The present series of experiments showed that all bands at and 
below 4:3 ev initially decay quite rapidly in the dark, but stay fairly 
constant subsequently. For example an absorption coefficient starting at 
2 cm~! decayed to 1-1, 0-9, 0-86 em~! in 0-5, 18 and 66 hours respectively. 
Hence all measurements on 2537 A activated crystals were made after at 

least half an hour, when the greatest rate of decay had passed. 

Further irradiation with light anywhere in the region of absorption up 
to 4ev (3090 A) causes all dark-stable bands to bleach, the rate of 
bleaching being greater for the higher energy bands. This is shown in 
fig. 4 (a), (b) and (c), where the rates of bleaching with light of 6000 A and 


Fig. 4 
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Change in absorption coefficient at various photon energies due to (a) heating 
at 220°c, (b) bleaching at room temperature with 5-7 10-4 watts of 
6000 A wavelength light, and (c) bleach at room temperature with 

10-3 watts of 3650 4 wavelength light. 


3650 A wavelengths are compared with the rate obtained by heating at 
220°c, at which temperature the absorption changes are of the same order. 
If it could be assumed that the oscillator strengths of the various bands 
were the same, then this would mean that more centres are disappearing 
at higher wavelengths and indeed this is the interpretation made by 
Molnar and Hartman (loc. cit.) in respect of X-irradiation damage. On 
this basis it would appear that for 6000 A light or thermal bleaching 
‘(during most of the first two hours of treatment) the integrated drop in 
absorption coefficient of the 4-3 ev and 2-35 ev bands at any instant is 
approximately equal to that of the 5-5ev band. This is not true for 
3650 A irradiation however. 
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In no case has all the ultra-violet induced absorption been caused to 
decay completely by optical bleaching, though after a week’s irradiation 
under a high-pressure mercury are lamp all bands except those at and 
above 4:3 ev had disappeared and these latter had been reduced to below 
5°, of the saturation absorption. 

Thermal bleaching results are erratic and complicated. The absorption 
of all bands in the region 1-0 to 6-0 ev decreases initially. At temperatures 
below 100°c a partial recovery is sometimes observed after some minutes 
and then the decay continues. This behaviour has not been observed at 
temperatures above 100°c, and typical bleaching behaviour at 220°C is 
then shown in fig. 4. 

This thermal bleaching is accompanied by a luminescence, though at 
100°c the eye cannot detect it. At 220°C the luminescence last about 
two minutes. As far as can be ascertained with the aid of a small spectro- 
meter the luminescence spectrum covers much of the visible spectrum 
with a predominance of red. Unfortunately this luminescence was 
detected only when the rates of decay of nearly all the bands was a 
maximum, and hence it could not be uniquely associated with the decay 
of any one band from these observations above. 

If a crystal about $ cm thick be exposed on one face to 2537 A wave- 
length light and the induced absorption along a section normal to the face 
examined, it is found that there is a region of darker absorption, followed 
by a uniform coloration over the rest of the section. This region of darker 
absorption covers about 2mm and so all experiments in the growth and 
decay of ultra-violet induced bands were made on specimens less than 
2 mm thick. 


4.2. Neutron Irradiation 


Neutron irradiation produces bands with apparent maxima at 1-2 ev 
(10 300 A), 2-15 ev (5750 A), 3°5 ev (3530 A), 4-2 ev (2940 A), 5 ev (2470 A) 
and 5:5 ev (2250 A). This irradiation-induced absorption is shown in 
fig. 5 where the resultant curve is broken down into various suggested 
components. This curve resulted from an integrated neutron flux of 
2 LOM ay vats 

In addition to the very apparent bands, the absorption between 3-5 ev 
and 2-15 ev may be broken down into two further bands at 3-0 ev (4120 4) 
and 2-7 ev (4580 A). 


4.2.1. The 2:15 ev band 


The 2-15 ev band has a shoulder on its low energy side, which does not 
appear convincingly for doses lower than 1017 n.v.t. ; a dose of 102° n.v.t, 
makes the presence of this shoulder very definite however. For doses 
between 101? and 10° n.v.t. the shoulder rises and falls at the same rate 


as the 2:15 ev band maximum, and on thermal annealing both fall at the 
same rate. 
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a Fig. 5 
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The optical absorption induced by a dose of 2x10!7n.v.t. In the profile 
analysis it has been assumed that the bands induced by ionizing radiation 
are present. The inset figure shows the suggested constitution of the 
2-15 ev band ; the absorption coefficient scale has been made larger, and 
the spectrum corresponds to a dose of over 1018 n.v.t. 
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Curve A is of the unirradiated specimen. Curve B is after pile irradiation of 
41018 n.v.t. Curve © is after subsequent optical bleaching. Tem- 
perature of irradiation=21°c. 
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On optically bleaching with 3500 4 light, this band increases in its 
height, a typical increase being from 4-5 cm~? to 7 cm~tas shown in fig. 6. 
this irradiation is accompanied by an intense red luminescence which 
lasts longer than the band change for any given 3-5 ev irradiation rate. 
If a neutron-irradiated crystal which has been subsequently irradiated 
with 3500 A light to produce a greater absorption at 2-15 ey, is left in the 
dark for several weeks further absorption changes are observed. The 
absorption changes are shown in fig. 6. 

The thermal stability of the 2:15 ev band is very neutron-dose depen- 
dent. Its rate of decay at 230°c, after optical bleaching, is shown in 
fig. 7, and this corresponds to a dose of 2-2 101 n.v.t. Its behaviour 


Fig. 7 
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Decay of 2-15 and 3-5 ev bands with time, on heating at 230°c. The erystal 
had received a dose of 2-2 10!8n.v.t., and had then been irradiated 
with 3500 4 light before each anneal. The decay in a crystal not sub- 
jected to this latter treatment was found to be erratic. 


after a dose of 2-1 x 10” is shown in fig. 8. At 230°c there is a general 
instability in the band height, and this instability continues until a 
temperature of about 400°c when the absorption begins to rise. This rise 
continues until a temperature around 600°C is reached, when a rapid 
decay commences. Complete decay may be achieved in a few seconds at 
1200°c. This behaviour has been observed in two crystals which were 
separately irradiated in the Windscale pile ; here the ambient temperature 
was between 80°C and 100°c in the facility used, being slightly higher than 
the corresponding figures for the BEPO pile. 
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Irradiation with 25374 light at any stage in this thermal annealing 
process brought a slight rise in the height of the 2-15 ev band, but similar 
irradiation with 3500 4 light produced no detectable effect. 

The shoulder on the low energy side of the 2:15 ev band, which rises 
and falls with the band maximum is not due to the tail of the 1-2 ev band. 
This is evident from an examination of the position and width of the 1-2 ev 
band (fig. 5). Some experiments on crystals which were rejected under 
the criterion expressed in eqn. (1), had previously been examined and in 
some of these, two distinct bands around 2-1 ev had been produced by 
3650 A irradiation. This suggests that the shoulder may be the result of 
superposing two bands very near each other, as shown inset in fig. 5. 


4.2.2. The 3-5 ev band 


As shown in fig. 7 for a 101% n.v.t. dose this band initially shows quite 
different annealing behaviour to the 2-15 ev band, but after 20 min at 
230°C it settles down to a decay rate identical to that of the 2-15 ev band. 
Both of these decays are taken for crystals which contain oxygen excess 
before irradiation. 


Fig. 8 
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The isothermal annealing behaviour over a range of temperatures. Each 
isothermal is for a period of 2 hours. The curves are due to bands at 
2-15, 3-5, 1-2 and 3-0 ev reading from top to bottom. 


After a dose of 2x 102° however this band behaves in a very similar 
manner to the 2:15 ev band. The peak absorption rises until a tempera- 
ture of around 600°c is reached when a rapid decay commences (fig. 8). 
The general variability in the absorption coefficient during heating from 
100°c—400°c is also noticed for this band. . 

The effect of bleaching with 3500 A light is to cause this peak absorption 
to decrease over a period of several hours when 10°? watts em™ sec”? 
(5 10% photons em~? sec“) is falling on the surface of the crystal. The 
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shape of the band is not altered however, and on leaving the crystal in 
the dark for several weeks, the band peak partially regains its previous 
height as shown in fig. 6. 

A crystal which has been the vacuum annealed prior to irradiation will 
show differences in the rate of decay of the band maxima as compared 
with the results of fig. 7 for doses of up to 10!%n.v.t. This is partly 
because in an ‘ excess-oxygen ’ crystal the width of the oxygen band at 
4-2 ev measured | day after neutron irradiation is not changed sufficiently 
at the tail to affect the decay rate of the 3-5ev band. Hence the contri- 
bution of the 4-2 ev band tail to the 3-5 ev absorption may be eliminated. 
In the vacuum annealed crystal, however, the whole of the 4:3 ev band 
is irradiation-induced and so decays with the 3-5 ev band, and affects the 
latter. 


4.2.3. The 3 ev band 


This band is very apparent for neutron doses in excess of 101* n.v.t., 
_and for a dose of 2 x 10?° its annealing curve is shown in fig. 8. It appears 
to be slightly more stable than the 2-15 ev and 3 ev bands since it com- 
mences to decay at a higher temperature and its decay rate is not so rapid. 
Ina crystal subjected to a dose of 10!’ n.v.t. this band is barely apparent 
and may be optically bleached, until it is no longer apparent in the resultant 
spectrum. For a 2x 107° n.v.t. dose it is optically more stable, decaying 
by about 10% after prolonged exposure. 


4.2.4. The 1-2 ev band 


This band appears on neutron irradiation but differs from the other 
bands in that it does not show a consistent growth on annealing of heavy 
damage. Its behaviour is shown in fig. 8. Its reaction to heating up to 
600°C is extremely erratic, and at this temperature it commences to decay. 
The decay between 600°C and 800°c can however be reversed by the action 
of 2537 A light, but above 900°c the decay is irreversible and the band 
disappears completely after a few minutes at 1000°c. 

After a dose of 2x 10!% n.v.t. the band rate of decay at 300°c is very 
similar to the 3-5 ev band but there is a wide scatter of readings. At 
this dose the band may be 80°% bleached by exposure to a high intensity 
mercury-are source for about a day. 


4.2.5. Bands at and above 4:2 ev 


For doses of 10'*n.v.t. and upwards these bands absorb over 99-9°%, 
of the light incident at these photon energies ; there is an accompanying 
fluorescence which makes absorption measurements less meaningful. 
This fluorescence is at least partly associated with radiation damage, and 
so even the use of a magnesium oxide crystal as a spectrophotometric 
standard does not assist in interpretation. Boyd et. al. (1947) have 
observed absorption above 4-3 ev for small neutron doses; from their 
analysis it appears that the 5ev band may comprise two components, 
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one of which is very unstable decaying in a few minutes at 600°C, and 
the other following a first order decay process subsequently. 

All these annealing results apply to heating in a vacuum better than 
10-> mm of mercury. 

If the crystal has these annealings carried out in air, the results are 
similar to those mentioned above, up to 650°c. At this temperature the 
absorption decay rate lessens, and at 800°c a band at 4:2 ev begins to 
grow rapidly. Both the 2-15 ev and 3 ev bands will continue to decay 
however. These results apply to a crystal which has been vacuum an- 
nealed, prior to any irradiation treatments, to remove excess of oxygen 
in the lattice. 


§ 5. Discussion 


Attempts have been made to show thermodynamically whether 
Frenkel or Schottky defects exist in thermodynamic equilibrium in ionic 
crystals. Jost (1952) has made a calculation assuming the ion vacancy 
radius to be half the normal ion distance. This vacancy radius is 
critical however, and Mott and Littleton (1938) improved the argument 
by an exact determination of the vacancy radius. The application of 
these calculations to the alkali halides has shown that Schottky defects 
are expected to exist in equilibrium at all attainable temperatures, and 
experimental measurements have confirmed this. It is not so easy to 
make a similar estimate in the case of magnesium oxide however, since 
quantities such as the critical ion vacancy radius and lattice energy will 
depend upon the degree of ionic character, which is unknown. However 
following Mott and Littleton’s arguments, Yamashita and Kurosatawa 
(1954) have shown the presence of Schottky defects is not energetically 
unreasonable, and experimental evidence to support their predominance 
in magnesium oxide will be mentioned in this text. 


5.1. Non-Stoichiometry and Impurities 

On heating crystals of magnesium oxide in the vapour of oxygen or 
magnesium, bands appear at various photon energies. It is important to 
arrive at some conclusion concerning the types of centre giving rise to 
these bands. In the case of magnesium addition, bands are produced at 
2-35 ev, 3-6ev and 4:95ev. Bands having an identical profile with 
maxima at these same energies are also produced by potassium addition ; 
hence it is clear that the spectrum is not due to resonance transitions at 
magnesium or potassium atoms because the resulting spectrum would be 
different in each case. It is probable that the spectrum arises from the 
added magnesium atoms becoming thermally ionized and thereby giving 
up their electrons to the interior of the solid. These would become 
trapped at existing centres. An examination of the absorption spectra 
produced in various planes at right angles to the direction of diffusion, 
shows that the spectrum does not fully develop in the crystal interior, 
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unless heating is carried out for a considerable time. One explanation of 
this is that the spectrum is due not only to electrons trapped at existing 
centres but also to those trapped at new centres which take time to 
diffuse inwards from'the surface ; the electrons themselves would reach 
the centre of the crystal in a few seconds at 1200°c, because it is known 
from irradiation experiments that this type of centre is highly unstable 
at a few hundred degrees centigrade. The indication is that Schottky 
type defects are formed by the addition of magnesium to the crystal ; the 
absorption is due to electrons trapped at oxygen-ion vacancies. Varley 
(private communication) has suggested a reason for the electrons ap- 
parently not diffusing at a very fast rate into the crystal interior. After 
ionization of the added magnesium atoms, the electrons would be free to 
diffuse away. However the effect of a number of new positive ions near 
the surface would be the setting up of a space charge which would tend to 
localize the freed electrons, and keep them in traps near the space charge. 
As time progressed the ions forming the space charge would gradually 
diffuse into the crystal interior and so allow the electrons to do so too ; 
eventually the space charge would be homogeneously distributed through- 
out the crystal. Even on this scheme however the absorption is still due 
to electrons trapped at negative ion vacancies. 

In the case of oxygen addition to the lattice, it is hard to distinguish 
from the data itself whether the absorption is due to neutralized mag- 
nesium ion vacancies or to interstitial oxygen atoms. The oxygen atom 
is about 1-2 A diameter, smaller than the magnesium ion, and so size 
considerations favour a high mobility. Furthermore Weber has pointed 
out that the averaged frequency rates of certain groups of 2p—3s and 
2p—3d transitions occurring in oxygen atoms in the gaseous state is the 
same as that of the two bands formed by oxygen addition to magnesium 
oxide. He suggests that the transitions behind these two bands may 
represent the same transitions as occur in the gaseous state but perturbed 
by the crystalline field. 

The evidence from the production by ultra-violet irradiation of ab- 
sorption at the same photon energies as the oxygen band is also relevant. 
It will be shown later that the ultra-violet irradiation produces free positive 
holes or electrons and that the trapping of these gives rise to the same 
absorption. Now if the absorption is due to resonance transitions at an 
interstitial oxygen atom, the oxygen must have been a singly charged 
negative ion in the vacuum annealed crystal. But this ion should 
itself give rise to absorption (an oxygen atom has an electron affinity of 
2-2 ev for one electron) whereas none is to be seen. 

Furthermore it seems rather unlikely that Schottky and Frenkel defects 
should occur together. It will therefore be assumed that additive 
coloration results in the formation of Schottky defects. 

If the absorption is not due to the oxygen atoms as such, but rather to 
trapped positive holes generated as a result of the addition of oxygen, 
the question arises as to where the positive holes would be trapped. In 
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addition to the suggestion already made of a positive ion vacancy, the 
possibility of impurity ions acting as trapping centres must also be 
considered and in the ultra-violet irradiation case this does occur. But 
whereas the holes trapped at the iron or chromium ions may be removed 
by heating the crystal in air at 600°c, the ‘ excess oxygen’ bands are not 
changed by this treatment. Of the other impurities the Ca?+ and Be?+ 
ions already have the inert gas structure and so we would not consider the 
possibility of their trapping positive holes ; Mn?+ has been shown to remain 
stable under irradiation and crystals have been obtained without any 
detectable cobalt or nickel, but their physical properties have been un- 
affected. This leaves the aluminium and silicon impurities to be discussed. 

Crystals with varying amounts of these have been obtained, but the 
saturation value of the optical absorption of the two bands under con- 
sideration has not varied. 

It is therefore concluded that if these ions may trap positive holes, then 
they do not give rise to absorption spectra at the same wavelengths as 
those of the excess oxygen bands, The overall conclusion is that these 
bands are due to positive holes trapped at lattice inhomogeneities, and 
are not associated with impurity ions in the lattice. 


5.2. Ultra-Violet Irradiation 

The absorption spectrum arising from ultra-violet irradiation has a 
profile very similar to that arising from x-irradiation. Im the case of 
x-irradiation Weber has showed that a superpositioning of the absorption 
spectra produced separately by magnesium and oxygen addition to the 
lattice, follow almost precisely the absorption profile produced by 
x-irradiation ; this makes it highly probable that the absorption bands 
arise from the same types of centre in the two cases. A similar inference 
may be drawn here in respect of the ultra-violet induced bands. 

Ultra-violet irradiation in the range 2000 A—3000 A will in general cause 
only electronic displacements in a virgin crystal. From Day’s photo- 
conductivity data a more precise idea of the exact type of displacement 
may be deduced, the experimental facts are there : 

Ultra-violet irradiation in the range mentioned causes increased 
photon absorption at certain wavelengths. Part at least of the change 
is reasonably stable at room temperature in the dark. During irradiation 
the photocurrent decreases from its initial value by a factor of two, but 
subsequent irradiation with light of lower photon energies causes the 
photocurrent to increase by a factor of a thousand at these wavelengths. 

The ultra-violet activating photon energy is about 4 ev (3120 A) and 
the enhancement of photoconductivity occurs with peaks at 1-2 ev and 
2-1 ev (10 500 A and 5800 A) ; this region of enhancement, may be shifted 
during its production, by the application of an electric field and it moves 
towards the negative electrode. 

These experimental results mean that the 4 ev radiation has caused a 
positive charge to become free to move in the lattice. The trapping of 
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this charge gives rise to the centre which can be photoactivated by light 
of lower energy. 

The positive charge is presumably a positive hole which is released 
by the excitation of an electron in a valence state of an ion, to a level 
where it becomes independent of the ion. Since the fundamental absorp- 
tion edge lies at some photon energy greater than 7 ey, direct excitation 
of the electron to the conduction band by such a low photon energy as 4 ev 
is not reasonable. 

Hence the electron excitation must somehow be associated with a 
lattice inhomogeneity which contains a trap level somewhere in the 
forbidden band of energy levels. There is initial absorption at the ultra- 
violet activating energy, but this shows no structure. This absorption 
appears to occur as long tail in the absorption edge in which case it is 
possibly due to exciton transitions. Any excitons formed may wander 
freely through the lattice and will either decay after about 10~* seconds 
or become annihilated at imperfections or dislocations or in the region 
of an ion vacancy. The latter might occur if the exciton wave were to 
pass a trapping centre with an electronic state of lower energy than that 
of the exciton ; in this case the exciton electron could tunnel through to 
a stationary trapping state leaving a positive hole free to wander in the 
lattice. The occurrence of some such exciton mechanism is supported by 
the increase in efficiency in the process with increase in irradiating photon 
energy, and the fact that the photon energy responsible lies in the absorp- 
tion tail. The direct excitation of an electron into the conduction band 
is not considered reasonable since there is good transmission up to 7-2 ey 
in the vacuum-annealed crystal. 

Since the crystal is essentially in its virgin state there are three pre- 
dominant possibilities for the trapping of an electron, 

(i) at a negative ion vacancy, 
(i1) at an interstitial positive ion or oxygen atom, 

(iii) at an interstitial impurity ion or atom. 

Following on the discussion under “ Impurities” a centre of type (i) is 
preferred. 

The hole which is released will jump from one lattice site to another 
very rapidly, making about 101 jumps per second. The simplest 
trapping centre it might finish at would be a positive ion vacancy. An 
electron in a valence state of an ion adjoining such a vacancy, would have 
a higher energy than on an ion in a complete part of the lattice. And 
since a minimum photon energy of around 2 ev is necessary in order to 
cause the positive hole to be released, the state is presumably 2 ev or so 
above that of a normal valence electron. Hence we may represent the 
electrons adjoining a positive ion vacancy as being in a state 2-1 ey above 
the valence band. It will be shown that this simple scheme fits several 
other observations on ultra-violet activated crystals. 

Alternatively the absorption at the irradiating ultra-violet photon 
energies could arise from centres with widely overlapping wave functions, 


Single Crystals of Magnesium Oxide 623 


which would tend to destroy any fine structure. In this case the tran- 
sitions would have to be to the conduction band, and the positive hole 
left behind would have to be free to wander about it in the lattice. 

On either of these schemes the photocurrent caused by ultra-violet 
irradiation is due to the subsequent movement of positive holes released 
at negative ion vacancies to positive ion vacancies. It is observed 
(Day 1953) that this photocurrent decreases by a factor of 2 after con- 
tinued irradiation ; this is due to the ‘ using up ’ of the vacancies present 
and to increased optical absorption around 4 ev. 

The decay of photocurrent at 2-1 ev takes place at a very slow rate 
occupying several hours. This is presumably due to the thermal acti- 
vation required to take a positive hole up to a negative ion vacancy. 
Because of this the probability for recapture at a positive ion vacancy 
would be greater than for returning to the negative ion vacancy. 

This scheme would mean that the level of activation made possible by 
ultra-violet light is proportional to the number of defect centres in the ~ 
virgin crystals. This is experimentally supported by both photo- 
conductive and optical absorption data. In the photoconductive case 
the activation level at photocurrent saturation is a factor of two lower 
after careful annealing ; this treatment would tend to reduce the number 
of defects (Day loc. cit.). The optical absorption level attainable by 
4-85 ev (2537 A) irradiation is reduced, by the factor of 4, after very slow 
annealing in vacuo from 1400°c. 

Since the electron on an ion in the neighbourhood of a positive ion 
vacancy has a raised potential energy, we would expect luminescence to 
be associated with the 2-1 ev centre production. Nothing is visible to 
the eye, perhaps due to the high light background from the irradiating 
source. However, Hisenstein (1954) has observed a multiband lumin- 
escence spectrum on X-irradiation, with a long-time phosphorescence 
at 3-6ev. After this phosphorescence was over the luminescence could 
still be stimulated by irradiation at around 2-ley. It is expected that 
x-irradiation would produce at least all the phenomena caused by ultra- 
violet irradiation, and so we would expect similar holes and electrons to 
be trapped. On the scheme under consideration then, the multiband 
spectrum during x-irradiation would be partially due to the hole trapping 
at positive ion vacancies, as well as simultaneous annihilation of F-type 
centres by other positive holes. The later stimulation of luminescence 
at 3-6 ev would then be the result of releasing the holes by 2-1 ev photons, 
and their recombining with electrons to give both luminescence and the 
observed optical absorption decay. 

A possible explanation may be given for the enhanced decay of absorp- 
tion at 2-lev by light irradiation plus simultaneous heating. Light 
irradiation alone may cause the holes trapped at positive ion vacancies 
to be released. These holes can either go to some impurity ion or to an 
occupied negative ion vacancy if they are to be trapped or annihilated 
without returning to a positive ion vacancy. However if the trapping 
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or annihilating centre carries an effective positive charge, the released 
positive hole will need a further energy to enable it to jump the resulting 
potential barrier. But many of the inhomogeneities supposed to exist 
may carry on effective positive charge, e.g. any impurity ion which may 
be capable of increasing its charge, Fe?+ etc., or a negative ion vacancy 
with one or perhaps even two electrons. Now the additional effect of 
temperature upon the kinetics or such a system would be to increase the 
rate of purely thermal release of both positive holes at positive ion 
vacancies and electrons at negative ion vacancies. In addition it would 
assist the rate with which released positive holes could surmount the 
potential barrier prior to trapping or recombination with trapped electrons. 
Hence in thermal bleaching alone the factor limiting the rate is the 
effect which the temperature has upon the release probability of a hole. 
For optical bleaching the factor limiting the rate is the effect of tem- 
perature upon the combination probability of a positive hole with an ion 
in the neighbourhood of a negative ion vacancy. Levy and Varley (1955) 
have observed a similar enhanced rate of bleaching for vitreous silica. 


5.3. Cathodoluminescence 


Cathodoluminescence results have been reported by Gandy (1955), 
Saksena and Pant (1955) and Woods and Wright (1955). The process 
involved in the formation of this luminescence is that electrons bom- 
barding the crystal cathode, penetrate the lattice and strip electrons off 


Table 2 
Source of Photon energy of luminescence 

information (max.) (e.V) 
Gandy L-7 | 1:9} 2:1 [2:21 2-4) 2-7) S36 15a 
Saksena and 

Pane 2-05 | 2-2 2:7 3:7 | 51 
Woods and é é 

Wright 1-95 | 2-04 | 2:3 2°8 | 3-1 


ions. These electrons are given sufficient energy to reach the con- 
duction levels where they are free to move until trapped at some centre. 
Part of the transition to this centre is non-radiative, because once the 
electron has left the ion, a lattice relaxation takes place. The remainder 
of the transition to a state of lower energy will possibly be radiative, and 
constitutes the luminescence. 

Luminescence in MgO has been detected at various photon energies and 
these are detailed in table 2. The 5-1 ev band has been shown by Gandy 
and Saksena and Pant to be associated with excess oxygen and the 2 ev, 
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2-4ev and 2-7 ev bands have been shown by Gandy, and Woods and 
Wright to be associated with excess magnesium. 

In interpreting the luminescence energies we would expect to be able 
to tie each to a transition from the final energy state of the released 
electron to an empty defect state. In general this energy will be less than 
that required to get an electron directly into the conduction band, 
because of lattice relaxation; the defect centre to which the electron 
moves will already be in a relaxed position. Hence each of the above 
emissive energies taken from the minimum energy of a released electron, 
in the conduction band, should represent the energy of various states 
measured from the top of the valence band. 

There is some evidence (Johnson 1954) that the absorption edge lies 
above 7-3 ev, and if the luminescence energies are taken from 7-2 ev there 
is a reasonable correlation as shown in table 3. 


Table 3 


Absorption band ryt Tet 
ect) 1-2| 2-1-2-35 | 3 |3-5-3-7| 4-3 4-8-4-9 5-6-5-7 


Emissive energy ‘ 9.9 « 
(ev) 5-1 3:6 & | 2-4, 2-2, 2-1 L7 
eo seve 2-1 3-6 |42/ 4-8, 5, 51| 5-5 


energies (eV) 


The emissive energy of 2:7 ev is not accounted for. This is consistent 
with there being transitions between the two states at 4:95 and 2:35 ev 
associated with excess magnesium. The luminescence band at 2-1 ev 
could also be made up of direct trap to valence level transitions. 

Saksena and Pant (1955 loc. cit.) have proposed a scheme for the 
observed cathodoluminescence in the case of an excess magnesium 
crystal based upon a conduction band width of only 5-9 ev. It is difficult 
to see how this can be so when it is possible to obtain good transmission in 
a stoichiometric crystal up to 7ev. However some of the transitions 
which they propose would be valid for the radiative decay of excitons 
formed in the tail of the fundamental edge or for the band to band tran- 
sitions associated with negative ion vacancies or other excess magnesium 
defects ; in this sense the present proposals supplement theirs by showing 
that the observed data and their proposals for band to band transitions 
are not inconsistent with assuming a forbidden band width of over 7 ev. 

Garlick (1956) has pointed out that since the ratio of dielectric constants 
for low and optical frequencies is (9-7/2-95), and the high temperature 
activation energy for semiconduction is 2-lev, the absorption edge 
should be around 7 ey as suggested by optical studies, 
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5.4. The Band Scheme Picture 


The band scheme picture which emerges is shown in fig. 9. If the 
optical transition width for the forbidden band be taken as 7-2 ev, then 
the positions of the 2:35 ev, 3:6 ev and 4:95 ev bands are unambiguously 
fixed relative to the top of the valence band. This is because transitions 
to the 4-95 ev, 3-6 ev and 2-35 ev levels from the top of the valence band 
are numerically the same as 2:35 ey, 3-7 ev and 4-95 ev transitions from 
these levels to the conduction band at 7-3 ev. Hence it has to be decided 
whether these transitions observed represent the movement of electrons 
to or from these levels. But since they have been associated with the 
trapping of electrons at negative ion vacancies, their representing tran- 
sitions to the conduction band is preferred. Hence the action of ae 
violet light is to fill these levels. 


Fig. 9 
ABSORPTION EMISSION 


CONDUCTION BAND 


VALENCE BAND 


The suggested scheme for absorption and luminescence. All transitions are 
given in electron-volts. 


Day’s results make it certain that the 1-2 ev and 2-1 ev transitions are 
due to the excitation of electrons into levels at 2-1 ev, which are full in 
the unirradiated crystal. The instability of the 1-2 ev and 4-95 ev bands 
at room temperature may be associated. The thermal excitation of 
electrons into the 1-2 ev band frees positive holes which are trapped at 
the 2-35 ev level; this reduces the number of electrons which can make 
up the 4-95 ev transition from this level. 

The 4-2 and 5-6 ev absorption bands are associated with positive holes 
trapped at positive ion vacancies. Hence the transtion is one from the 
valence level of an ion further away from the centre, to a level very near 
the vacancy. In an unirradiated crystal the valence band level will 
become locally raised due to the perturbing effect of the doubly-charged 
vacancy. ‘This is expressed in another way by noting that the lattice 
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relaxation will effectively change the ion-ion spacing locally, which will 
cause both the conduction and valence energy levels to be raised. Hence 
we would expect a transition to the perturbed 7:2 ev level to take place 
only if the ion was in a high vibrational state at the appropriate moment 
of absorption ; this would allow it to make the transfer to the minimum 
energy state for the conduction levels. For this reason the transition 
from the 4:3 ev and 5:6 ev levels to the 7-2 ev conduction band level 
would normally be forbidden, and these transitions are not in fact observed. 
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ABSTRACT 


Gamma-rays of energy 2-594-L0-015 mev have been observed following 
the inelastic scattering in 1°F of neutrons in the energy range 3-2 to 
3:7 Mev. They are attributed to a level at 2-791-.0-015 Mev which 
decays principally by a transition to the second excited state at 197 kev. 
The cross section for production of these gamma rays has been measured 
as a function of neutron energy and the result compared with the 
predictions of the Hauser-Feshbach model for different values of the 
spin of the level. It is concluded tentatively that the spin of the 2-79 Mev 
level is either 9/2 or 7/2. 


§ 1. InTRODUCTION 


InTEREST in the levels of 1°F has received a considerable stimulus from 
the intermediate-coupling shell-model calculations of Eliott and Flowers 
(1955) and Redlich (1955) ; they have predicted a number of properties 
for the even-parity states, including the spins, which follow the sequence 
1/2, 5/2, 3/2, 9/2. Paul (1956) has recently shown that very similar 
predictions can be made on the basis of Nilsson’s rotational model 
(Nilsson 1955). 

The data obtained from a variety of experiments concerning the 
positions, spins, and decay scheme of the first seven levels of 1°F are 
summarized in fig. 1.t The 197 kev level has been definitely established 
as a 5/2+ state and the 1-56 mev level is almost certainly 3/2+. The 
first excited state, at 110 kev, is 1/2—: it has been suggested that this 
state is like a ?°Ne nucleus with a hole in the p,,, shell of its 16O core. 
The levels at 1-35 and 1-47 Mev have been found both by inelastic neutron 
and inelastic proton scattering. Their spins and parities have not been 
established and the assignments shown in brackets are a guess which 
has been used in some calculation to be discussed below: the guess is 


+ Communicated by the Author. 

{ For references see the papers on gamma-rays following inelastic neutron 
scattering (Freeman 1956), on gamma-rays following inelastic proton scattering 
eae al. 1956), and on inelastically-scattered proton groups (Squires 
et ab. : 
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based on the suggestion that these two states may correspond to a 
**Ne* nucleus in its first excited state, again with a p,/, hole. 

The interest of the present paper lies mainly in the next level, at 
2-79 Mev, which is predicted to be a 9/2+- state. This level has been 
observed in particle emission in several experiments, the most accurate 
being the recent measurement of an inelastically scattered proton group 
by Squires et al. (1956); they obtained the value 2-784--0-008 mev 
for this level, which was the only one observed between the 1-56 Mev 
state and one at 3-9mev. In the experiments to be described here a 
gamma-ray of energy 2-594+0-015 Mev has been observed, following 


Fig. 1 


2-791 


259 


Ss 
Ls 
$+ 


Level and decay scheme of 1°F. 
Observed level energies and gamma-ray transitions are given in Mev. 


inelastic neutron scattering, the threshold being at a neutron energy 
somewhat greater than 3 Mev. Since there are no other known levels from 
which it could arise, this gamma-ray is attributed to the 2-79 Mev level 
with the assumption that it is produced in a transition to the second 
excited state at 197 kev. This would give a value of 2-791+-0-015 Mev 
for the level, in good agreement with the measurement of Squires ef al. 


§ 2. EXPERIMENTAL METHOD AND RESULTS 


The observations were made with a ring-geometry arrangement 
previously described (Freeman 1956). Neutrons, obtained from the 
T(pn) reaction with an energy spread of 70 kev, were scattered by a 
teflon (effectively CF.) ring of outside diameter 10 cm, inside diameter 
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5-5 em and thickness 2.cm, placed with its axis along the direction of 
the proton beam and its centre 36 cm from the neutron source. Gamma- 
rays were detected by a shielded sodium-iodide crystal, 38cm in 
diameter and 2-5 cm thick, situated at the centre of the ring. The 
crystal could thus receive gamma-rays emitted from the ring within 
the angular range 45° to 125° with respect to the direction of incident 
neutrons. The neutron flux was measured with absolutely calibrated 
BF, long-counters. Gamma-ray spectra were recorded on a hundred- 
channel analyser, first with the teflon ring in position, and secondly 
with a ring of carbon, of mass chosen to produce about the same amount 
of elastic scattering. The spectra were then subtracted to give the 
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Spectrum of gamma-rays following inelastic scattering of 3-56 Mev neutrons 
in 1°. The broken curve is the spectrum of a RdTh source of 2-615 Mev 
gamma-rays, fitted to the 2-594 mev fluorine peak. 


spectrum of gamma-rays due to inelastic scattering in fluorine. A typical 
result obtained in this way for a mean neutron energy of 3:56 Mev is 
shown in fig. 2. The energy of the peak at 2-594 Mev was measured 
by comparison with a RdTh source of 2-615 Mev gamma-rays. The 
broken curve shows the shape of the RdTh spectrum, adjusted to fit 
the 2:59mev peak. The Compton edge due to 2-59 Mev gamma-rays 
is evident, and also the position of the peak corresponding to the loss 
of one annihilation gamma-ray following pair-production in the crystal ; 
the peak corresponding to loss of both annihilation gamma-rays is 
somewhat obscured by the rise of the spectrum at lower energies. Evident 
here are the peaks of 1:36 and 1-24 Mev gamma-rays following inelastic 
neutron scattering to the lower levels of 19F. 
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The cross section for production of 2-59 Mev gamma-rays in the angular 
range of observation, at a particular neutron energy, was obtained by 
the following methods: (i) calculation of the efficiency of detection of 
the Nal crystal, (ii) comparison with a RdTh source of 2:615 Mev gamma- 
rays of known strength, and (iii) comparison with the yield of 1:36 mev 
gamma-rays following inelastic scattering to the 1:56 and 1-47 Mev 
levels ; the absolute yield of 1-36 Mev radiation was measured at lower 
neutron energies, first by a comparison with the reaction ®*Fe(nn’)5*Fe* 
for which the cross section is known (Freeman 1955 a), and secondly 
by using the measurement by Day (1956) of the yield of 1:36 Mev gamma- 
rays at a neutron energy of 2-56 Mev. The total cross section for 2-59 Mev 
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Excitation of the 2-79 Mev level in 19°F. The points give the observed cross 
section for 2-59 Mev gamma-ray production ; they represent the total 
cross section for excitation of the 2:79 mev level to the extent that 
other modes of decay of this level can be neglected. The full curves 
are calculations based on the Hauser—Feshbach model for the spin 
assignments shown ; the broken curve is a complex potential calculation. 


gamma-ray production was obtained, with an overall error estimated 
to be 40%, assuming the angular distribution of the gamma-rays to be 
isotropic to within 10%. 

The yield of the 2-59 Mev gamma-rays was studied as a function of 
the initial neutron energy up to the maximum energy obtainable with 
the Harwell Van de Graaff machine. The result is shown in fig. 3. 
Ignoring the curves for the moment, the black circles show the observed 
variation in yield with neutron energy ; the errors given on the points 
are statistical. The ordinate scale gives the cross section for production 
of 2:59 Mev gamma-rays; the total cross section for excitation of the 
2:79 mev level may be somewhat greater if gamma-ray transitions to 
levels other than the 197 kev level are not negligible. The observations 
show that 2:59 Mev gamma-ray emission is the main mode of decay ; 
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transitions to the ground and first excited states are together less than 
10°, but somewhat more branching to the levels at 1-35, 1-47 or 1:56 Mev 
could have occurred without detection. It is estimated that the values 
of the cross section for excitation of the 2:79 Mev level could be up to 
50% higher than those given in the figure for 2-59 Mev gamma-ray 
production. However the main feature of the experimental results 
remains, namely the very slow rise in yield for the first 300 kev above 
the threshold (2:94 mev in the laboratory system), which suggests a 
large spin difference between the 2-79 Mev level and the ground state. 


§ 3. THEORETICAL COMPARISONS 


An attempt has been made to determine the spin of the 2-79 Mev 
level by making a comparison with the Hauser and Feshbach model 
(1952). The basic assumption for application of this theory is that, 
at the given excitation energy of the compound nucleus and with the 
energy spread used in the incident neutron beam, many levels of the 
compound nucleus can be excited, so that resonance and interference 
effects are not encountered. The cross section for formation of a compound 
nucleus with neutrons of angular momentum / is assumed to be of the 
form o{)= (21+ 1)7A?7,(EH), where the 77(H) are transmission coefficients. 
In the decay of the compound state, competition between neutron 
emission to the various accessible levels of the target nucleus is considered, 
and the cross section for excitation of any particular level is obtained by 
summing over all possible /-values, channel spins, and J-values of the 
compound state. For our calculations the spins and parities assumed 
for the lower 1°F states were as indicated in fig. 1. 

In order to test the applicability of the model to our case calculations 
were first made for some of the lower 1°F levels. In fig. 4 the upper 
curve shows the calculation for excitation of the 1-35 Mev level and the 
points represent the observed cross section based on measurements of 
the 1-24 Mev gamma-ray by which this level decays (Freeman 1955 b). 
The lower curve gives the calculated cross section for excitation of the 
1-47 and 1-56 Mey levels taken together ; it was convenient to combine 
these since each level emits a 1:36 Mev gamma-ray (see fig. 1) and the 
experimental points were obtained by measuring all the gamma-rays 
emitted from both levels. It is clear that in the neighbourhood of the 
thresholds for these lower levels the model is not applicable, there being, 
apparently, insufficient states in the compound nucleus. However, by 
the time one reaches the neutron energy range 3-0 to 3-6 Mev, where the 
level of excitation in 2°F is 10 Mev, the theoretical curves seem to agree 
fairly well with the experimental points. This may be fortuitous, but does 
encourage one to apply the calculations to the 2-79 Mev level. In fig. 3 the 
full curves give the predictions of the Hauser-Feshbach model for different 
spins and parities assumed for this level. The radius of the 19°F nucleus 
was taken to be 39x 10-! cm and the strong coupling approximation 
was used for the transmission coefficients 7. An examination was 
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then made of the effect of introducing the parameters of the complex 
potential model (Feshbach et al. 1954) into the T, coefficients ; the 
broken curve in fig. 3 shows the results of one such calculation, for 
a square well, with Vj,=42 and ¢=0-05 for neutron energies up to 1 Mev 
and Vj=40 and ¢=0-07 for higher neutron energies}. The conclusion 
is that the experimental results cannot be fitted by curves of the 
Hauser—Feshbach type. It seems probable, therefore, that the assumption 
of a statistical distribution of levels in the compound nucleus is unjustified 
and that relatively isolated groups of levels are influencing the shape 
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Cross section for excitation of the 1-35 Mev level in 1°F (upper diagram) and 
of the 1-47 and 1:56 mev levels combined (lower diagram). The points 
were obtained from measurements of the gamma-rays and the curves 
from the Hauser—Feshbach formula. 


of the cross section curve. In this connection it should be noted that, 
although the total cross section is quite smooth in this region (Hughes 
and Harvey 1955), Marion and Brugger (1955) have reported some 
resonant structure in the yield of the 1°F(n, «)1®N reaction for neutrons 
above 3-5 Mev in energy. 

In view of the discrepancies one cannot make a definite spin assignment 
from the analysis of fig. 3; the behaviour of the cross section curve for 
the first 300 kev above the threshold suggests a spin equal to or greater 


+I am very grateful to Dr. W. S. Emmerich of Westinghouse Research 
Laboratories for making available to me some of his calculations of complex 
potential model transmission coefficients. 


634 The Excitation of a 2-79 MeV Level in YF by Inelastic Neutron Scattering 


than 9/2 for the 2-79 mev level, whereas the trend at higher energies is 
more appropriate to a 7/2 spin. Tentatively we conclude that this 
level most likely has a spin of either 9/2 or 7/2; this would be consistent 
with the 9/2+ shell-model prediction. ; 
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ABSTRACT 


The effect of local deformation on the growth features on crystal faces 
has been studied on crystals of lead nitrate and potassium di-hydrogen 
phosphate. New growth centres form at these points and these can be 
initiated where desired on the face and if strong enough will dominate 
over the existing centres of growth. This can be done repeatedly, and 
it is suggested as a possible method for controlling flawed regions in 
artificially grown crystals. 


§ 1. [ytRopUcTION 


THE natural growth phenomena of lead nitrate crystals when growing 
from aqueous solution have been described in a previous paper (Williams 
1957). Following the work previously carried out on the local deformation 
of cadmium iodide crystals, Korndorffer et al. (1952), it was decided to 
carry out similar experiments on lead nitrate and potassium di-hydrogen 
phosphate crystals as these show clear layer formation on the growing 
faces. It was hoped that by this method some study might be made 
of the movement or creation of dislocations in the crystals. 

The method used to observe and record the crystal growth using 
transmission phase contrast has already been described in the previous 
paper. The deforming tool was a fine glass needle. This was attached 
to the base of the condensing system on the microscope and could be 
raised and lowered with it, to enable the point to be brought slowly into 
contact with the upper growing surface of the crystal while under constant 
observation. Optical adjustment was made with the indenter raised and 
out of the field of view and the crystal was observed under these con- 
ditions apart from the time during indentation. The indenter being 
made of glass did not affect the illumination appreciably. 


§ 2. OBSERVATIONS 


Figures 1-9 (Plate) show a sequence filmed of a crystal of lead nitrate 
which was subjected to deformation. Before deformation the crystal 
was growing by the initiation of growth layers at the lower left-hand 
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corner of the face under observation. The layer grows rapidly along 
the edges of the face and finally forms an enclosed pit (fig. 1"(PL)ye Phas 
behaviour would continue throughout subsequent natural growth. As 
described in the previous paper it is concluded that this is growth by 
normal nucleation and not based on dislocation centres since these 
would move to the centre of the face at a later stage in the growth and 
this does not happen. Also, no evidence of the overlapping of a layer . 
on one face on to another face has been found. 

The crystal was then deformed and fig. 2 (Pl.) shows the fine point 
in contact with the surface. It appears very large relative to the crystal, 
but in fact contact is only made in a small area. The ‘ point’ must be 
an irregular circularly shaped section of the filament from which the 
point is made and the part of the rim of this which first makes contact 
with the crystal surface acts as the indenter. That this is so is seen 
in the following fig. 3 (Pl.), which is that taken immediately after the 
removal of the indenter, the time of contact being made as short as 
possible. The three black spots mark the points of contact. The natural 
growth layers can still be seen following their original pattern. For 
transmission phase contrast the high edge of a step on the surface appears 
dark and the low side appears bright. The next two photographs, 
figs. 4 and 5 (Pl.), show that the points of contact becomes an active 
growth centre. The old natural growth layers complete themselves on 
the right-hand side of the face but the next natural layer from the left 
is overwhelmed by a larger step from the new centre. 

The new centre is complex at first, but later, fig. 6 (Pl.), it becomes 
steady, regularly forming layers in the form of closed loops. The 
sequence of photographs, figs. 6-8 (Pl.), shows the formation of one of 
these loops, which appears to be by the action of two screw dislocation 
- sources of opposite hand. In fig. 8 (Pl.), the area inside the loop is quite 
clear of cracks etc., which shows that the deformation which was produced 
was very slight, yet it was completely effective in making the point of 
indentation the initiation centre of growth layers. It might be expected 
that this effect would only be a temporary one, which would die out 
after the damage had been healed. However layers continued to form 
from this point, and it provides the source of the future growth of the 
crystal. This behaviour has occurred on all crystals studied. Once 
deformed the subsequent growth is entirely from the deformation centre. 
When natural growth is taking place from a centre within the face before 
deformation, the point of deformation still becomes dominant and the 
natural centre dies out. This is illustrated with a deformation centre 
in figs. 10 and I1 (PL), after the crystal previously discussed has been 
deformed a second time and a new growth centre induced to the left 
of the first one. In fig. 10 (PI), this has just been formed and the previous 
centre is still active, the loop on the right having been formed by it, 
but a few seconds later, in fig. 11 (Pl.), the new centre is dominant and 
remained so for the remaining growth. This repeated initiation of 
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growth centres, each in turn dominant could be carried out at will, and 
thus one could locate the centre of growth of a crystal at any desired 
point. 

The behaviour of the layers at the deformation centres often appears 
to give evidence of screw dislocation, in that they often show a spiral 
form. Figure 9 illustrates this on the crystal of lead nitrate already 
described. The interesting thing is the manner in which these shapes 
form. The steps on the crystal surface are formed by the bunching 
together of smaller layers formed at the complex centres as described 
in the previous paper. This also happens at the deformation centres. 
Figures 12-14 (Pl.) show another crystal of lead nitrate after deformation. 
Before deformation the growth layers formed by nucleation, as described 
for the previous crystal. After indentation the growth was complex for 
a few seconds, layers of the closed loop type forming. In fig. 12 (PI.), 
the left-hand edge of the loop in the process of formation is thicker than 
the right-hand end and this causes incomplete joining of the layer so 
that the thicker turns and dominates the screw dislocation of the opposite 
hand leaving a spiral form on the surface of the crystal. This persists 
strongly for some time, all the growth taking place from the dominant 
screw dislocation centre. The behaviour is periodic, however, and after 
a certain time, about one minute, the centre becomes quiescent and the 
growth dies out. No more layers appear to form, the existing layers 
grow out and the surface becomes flat. The centre as a whole then 
becomes active again repeating the procedure of forming initially closed 
loops ; then the left-hand arm builds a thicker step and a spiral form is 
again set up. This periodic behaviour was observed to repeat at intervals 
of a few minutes. This is explained by the fact that the centre is complex 
and at one time the elementary layers reinforce and bunch together 
preferentially in one area but cancel one another in another area, as 
described by Frank (1949) in his description of the behaviour of several 
centres within an area. Such an array of dislocations must have been 
generated during the deformation process. The combinations of the 
elementary layers depends on their thickness, and the separations of 
the centres in the local area. At certain instants they will combine 
sufficiently to cause spiral growth forms, at others they will cancel and 
the face will appear to have stopped growing with visible layers. 


§ 3. MuutTreLe DEFORMATION 


This has already been mentioned with the first crystal of lead nitrate 
described. Repeated deformation is also illustrated in figs. 15-18 (PL), 
showing the effect of deformation on a crystal of potassium di-hydrogen 
phosphate. Figure 15 (Pl.) shows the layer formation from the first 
deformation centre made in the centre right where the crack mark at 
the point of indentation can be seen. Figure 16 (Pl.) is one taken some 
minutes later when an indentation had been made lower and to the left 
in the face. The crack mark associated with this can be seen and also 
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that due to the first deformation. The new centre dominates the growth 
on the face, though the old centre is still active, modifying the layers 
passing over this region. Figure 17 (Pl) shows the effect of a third 
deformation in the centre of the face to the left of the first point of 
contact and at the same time of a fourth centre initiated at the top of 
the face. These are now acting together and the influence of the second 
centre is still seen. The last photograph shows a final deformation made 
at the lower edge which has dominated all previous centres. The crack 
marks of the old centres are absorbed into the body of the crystal. The 
growth centre can be changed at will in this way, the strength of the 
deformation determining whether the effect of previous centres will 
remain. 

In this way it might be possible to control the growth centres on the 
faces of artificially grown crystals. The growth in crystal faces seem from 
the evidence of the experiments and previous studies on natural growth 
to be associated with flaws or imperfections which are propagated through 
the body of the crystal as growth proceeds. To control the growth 
centres gives a means of producing unflawed regions in the crystal which 
will be formed by layer deposition away from the initiating centre of growth. 

When these experiments were commenced it was hoped that it might 
be possible to move dislocation centres of growth so that the movement 
and resulting slip traces might be seen. Efforts have been made to 
produce movement by stressing cadmium iodide crystals which exhibited 
spiral layer growth. These experiments have not been successful. A new 
centre is generally formed at the point of stress, but the spiral centres 
were not affected. On lead nitrate the original centres were generally 
overcome by layers from.the new centres formed, the original centres 
dying out immediately. On cadmium iodide also, after some time if 
strong enough, growth from the new centre will cover the crystal though 
the later stages of growth are often very complex on this crystal. For 
potassium di-hydrogen phosphate the new centres took a short time to 
establish themselves and did not always overwhelm the old ones. No 
slip was ever observed. Deformation did not appear to have any effect 
on the stacking fault lines, which were described in the previous paper. 
During the deformation local movement of dislocations was probably 
taking place but on a small scale probably immediately masked by the 
subsequent growth and bunching which occurs in the formation of 
visible growth layers. 
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SUMMARY 

The assumptions involved in two well-known formulae due to Griineisen 
are investigated and with these assumptions an explicit formula for the 
volume of a solid as a function of temperature is derived. Application 
is made to the case of KCl and the results compared with experiment 
with a view to deciding whether defects are present at high temperatures. 
Defects in the theory and a previous investigation by Fischmeister are 
discussed, 


§ 1. IyrRopuctTION 


Ir has been suggested that when a solid approaches its melting point 
lattice defects such as vacancies and interstitial atoms may occur and 
that these could be responsible for the fairly rapid increase in coefficients 
of thermal expansion « in such temperature regions. Lawson (1950) 
has assumed that « should be constant for a perfect lattice and therefore 
deduces that defects do occur in the cases of AgCl and AgBr. Gertsriken 
(1954), in determining the number of holes and energies of hole formation 
in a number of metals and alloys, has tacitly made the same assumption, 
as has Uno (1951) in the case of NaCl. No theoretical reasons for this 
are given and it seems clear that the variation of volume with temperature 
predicted by theory for a perfect lattice is not fully understood. Obviously 
this point must be clarified before any conclusions can be reached 
regarding the presence of defects. Fischmeister (1956) has attempted 
to show that the theories of Debye and Griineisen account completely 
for the expansion of the alkali halides so that no defects need be present 
but it is shown here in §3 that his method is too insensitive to be 
conclusive on the point. In §§2 and 6 the assumptions involved in 
deriving two well-known formulae due to Griineisen (1926) are investigated 
from both the mathematical and physical points of view and in $4 a 
relation is derived between the volume and temperature of a Debye— 
Griineisen solid involving the same assumptions. Calculations are carried 
out for KCl and the results are compared with experiment. Factors 
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which might be responsible for some of the wide discrepancy between 
theory and experiment are discussed and possible improvements to the 
theory are mentioned. 


§ 2. GRUNEISEN’S EQUATIONS 


Let ®(V,7) be the potential energy of volume V of a solid at 
temperature 7’ with all atoms at rest in their equilibrium positions 
and 7(V,7) the free energy of the vibrations of the atoms about these 
positions so that the total free energy of the solid is 
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Then the external pressure p on the solid is 
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and the equation of state of the solid may be put in the form 


pV+G(V, T)=vE(V, T) Pie, ete Ney eae oe 
where £ is the internal energy of the vibrations and G, y are defined by 
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No assumptions have been made so far although, as mentioned later, 
it may not be convenient or even possible to divide the free energy of 
the solid into two parts © and Z as defined above. Now assume 
(a) that ® is a function of V only, independent of 7’, so that G is also 
and (6) that (dy/d7');=0. Then, differentiating (2) with respect to 7 
at constant volume 
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But the specific heat at constant volume, coefficient of thermal expansion 
and isothermal compressibility are defined as 
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and hence (5) may be written 
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Alternatively, with assumption (a) only, consider the equation of 
state (2) when p=0. Let V, be the equilibrium volume of the solid 
when 7’=0 (suffix zero will always have this meaning) and expand G(V) 
as a Taylor series, putting V—V,—4 
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If it is assumed (c) that the atoms are at rest in their equilibrium 
positions at 7’=0, a point to be considered if quantum theory is 
introduced, then for equilibrium one must have 
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and hence G(V)) is zero. It may easily be shown (e.g. Roberts 1940) 
that 
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Thus (2) becomes 
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To a first approximation 4=yHx,; substituting this in the square 
brackets on the left-hand side the result obtained is 
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Thus assumptions (a@) and (6) are necessary in order to derive eqn. (6) 
and (a) and (c) to derive eqn. (10). Before considering the physical 
meaning of these assumptions a little more may be said about eqn. (10). 
For a monatomic solid treated as an elastic continuum by the Debye 
theory 
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where D(@/T) is the Debye specific heat function. At reasonably high 
temperatures this function is essentially constant so that, with the 
assumption that y is independent of temperature H/Q)=AT where 
A=—3(k/Vo)yRD(9/T) is approximately constant. Then 

Aten Alt 

Vee ATE 
and A/V, will increase faster than linearly with 7. This will be true on 
any theory provided that for the temperature range considered 
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OV 
where B is a constant (necessarily positive) and a, as defined by (9), Is 
positive. With regard to the latter Griineisen (1926) assumed for two 
atoms a distance p apart a potential of the form 
‘ a b “ 
=——+— APE CS, te ele) 
? pare 
and obtained «=(n-+-m-+3)/6>0. 
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§ 3. FISCHMEISTER’S INVESTIGATION 


Fischmeister (1956) has investigated the validity of eqn. (10) for the 
alkali halides when the Debye expression 6R7T'D(@/T) is used for the 
internal vibrational energy £ per mole of a diatomic molecular solid. 
By inverting the equation, he obtains 
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where a is a linear dimension of the solid. On the theory of Debye and 
Griineisen Q = Vo/y«p is independent of temperature so that, if @/(a@,;—4p) 
is plotted against 1/7'D(@/T), a straight line should be obtained, ignoring 
terms after the second on the right-hand side of (13). In each case, 
using a constant experimental value of @, Fischmeister finds that the 
experimental results for a@/(@,—4d ) do lie on a straight line and concludes 
that the dilatation of the alkali halides does vary in the manner predicted 
by Debye and Griineisen. There would then be no need to invoke the 
occurrence of lattice defects to explain the rapid increase of dilatation 
near the melting point. 

Two points in Fischmeister’s treatment may be thought to have 
influenced the results. The first is the approximation d)/(@;—d))=3V 9/4 
and the second is the use of a constant value for @ over the whole 
temperature and volume range. In the case of KCl the calculation was 
repeated by the present author, correcting the first approximation and 
allowing © to vary as V~°/? (see (19)) but no appreciable change in the 
graph was noted. 

There still remains the point that the temperature scale of these graphs 
is far from linear so that any deviation from a straight line in the important 
high temperature range would be difficult to detect. The question is 
whether this method of plotting the experimental results is capable of 
distinguishing them from the results of a theory which postulates no 
defects when there are in fact defects present. In order to test this, 
graphs were plotted of a /(ap7—d )) against 1/7.D(O/T) assuming (i) a 
constant coefficient of thermal expansion of 0-00014/°c for the range 
273<7'<1000°K and (ii) a coefficient «=6 x 10-°+1-:18x 10-7 for this 
range, the actual figures being so chosen that a valid comparison 
could be made with Fischmeister’s graph for KCl. The former is the 
type of assumption made by Lawson (1950) for a perfect lattice and the 
latter corresponds to assuming that a theory for a perfect lattice will 
give « increasing linearly with 7’ for medium and high temperatures 
and any deviation from this is due to the occurrence of lattice defects. 
The results are shown in fig. 1 and it is seen that the Fischmeister method 
shows clearly that a Debye—Griineisen solid could not have a constant 
coefficient of thermal expansion but it fails to distinguish between the 
actual solid and one which has « linearly dependent on 7’, i.e. there 
could be an occurrence of lattice defects at high temperatures in the 
actual solid. 
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For this reason it was decided to try to obtain a direct relationship 
between volume and temperature for a solid, based on the theories of 
Debye and Griineisen, from which a graph of « against 7’ could be plotted 
for direct comparison with experiment. ji 


Fig. 1 
(50 
100 
50 
{ 
| 2 5 4 5 


108/TD(O/T) 


Expansion of KCl. Experiment (Eucken and Dannohl 1934) ; 
66 Assuming «=0:00014/°; * x Assuming «=6 x 10-°-+1-:18x10-77'/°. 


§ 4. THEORY OF THERMAL EXPANSION 
When p=0 the equation of state (2) becomes 
G(V)=y#. Ad eae? CP) 
If we assume an interatomic potential of the form (12), it is easily shown 
(Griineisen 1926) that for a monatomic solid 


-y r\ 17 oP if mA (pe m/3 Lies ; 
G(V)=I (Fr), = a5 (7) -(R)"] . aa) 


where A=aN("+)3g, N is the number of atoms per mole and 
assumption (c) of § 2 has been made. Provided it is assumed that no 
distortion of the lattice structure occurs as the volume changes, o,, 
depends only on the type of structure and on m, being defined by 


Om= > (r/p)™, where V=Nr? and the summation is over all pairs of 
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atoms in volume V, each pair being counted once. On Debye’s theory 
for a monatomic solid treated as an elastic continuum # and y are 
given by 


(16) 


FP ae 
Vv 


me Jo OXp (Av/kT)—1’ 


Olnv,, be 
=- (FEF), er a G7) 


v, lees: ae Vee CVE, | ee, 


c, and c, being the velocities of longitudinal and transverse waves 
respectively. It is not a good enough oe ae to assume ¢,, C; 
independent of volume as this leads to y=}, a value generally regarded 
as much too small for any solid. Griineisen (1926), by considering the 
motion of one atom with all others fixed, the interatomic potential 
being of the form (12) with nm, has obtained v as a function of V’ 
for small amplitude vibrations and thereby deduces y=(n-+-2)/6. As a 
first approximation this seemed sufficient for present purposes and 
then (17) and (18) give 


where 


Vm=CoVoreV—mie, |, 49) 


From (7) and (8) it is easily shown that A=9V,'"+*)3/m(n—m)ky 
substituting (15) and (16) in ry. one obtains finally 


yy | (“” (n+2)/6 x3 dx _ 
yr De ym yr =T* F et] (20) 
where 2=hAvy/kT; y=ViV, an 
20 he ee 
ae ennmsneatin Sel Sa Sena: Pa Vee us 9 
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§ 5. CALCULATION FOR KCl 

It was decided to apply the above theory to KCl since it was one of 
the substances considered in Fischmeister’s paper and the theory should 
be reasonably applicable to it for the following reasons. ‘ 

(i) It is an ionic substance with the electrons tightly bound to the 
nuclei in the crystal so that electronic effects not considered above 
should be small. 

(ii) An interatomic potential of the form (12) with m~1, n~8 seems 
to be generally accepted as a reasonable assumption for KCl; one notes 
incidentally that n> m. 

(iii) The above theory is strictly applicable only to monatomic solids. 
However, if one assumes a potential (12) between like ions and one of 
the form ap~-”+-bp-" between unlike ones with the same values of a, b. 
m and n, and one assumes that the two kinds of ions are of approximately 
equal mass as is the case for KCl, then eqn. (20) may be derived as 
before provided N is taken as the number of atoms in half the molar 
volume. 
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For the actual calculation the measurements of Durand (1936) 
of the elastic constants of KCl were extrapolated to 0°K to give 
(€11)9= 5-07 X 10™ 5 (€y9)9=0-6 X 1011 dynes/em?. Substituting these in 
the formulae 


3 3(1— 3(1—2 
k= ————. : Cra ( a) : Co aS ( o) 5 og tee ; (22) 
Cyy+2Cy9 (l+-o)Kp 2(1+o)Kp Cri tCie 
and assuming a value 2-07 g/cm’ for the density p, the following values 
were obtained 
Ko=4:8X10-Y% cm4/dyne;  (€,)p>=4:94X 105;  —_(c,)9>= 3-28 X 105 cm/sec. 
Using a mean ionic weight of 37-28 g a value of 18-01 cm? was obtained 
for Vj and then assuming m=1, n=8 gives 
J =5-25 x 101° deg. ; K=344 deg. 

For various values of y from 1 up to 1-12, eqn. (20) was then solved for 7 
by trial and error. A graph was plotted of In y against 7’, from which 
the coefficient of thermal expansion can be deduced since 
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The variation of « with temperature is shown in fig. 2. It is seen that 


Fig. 2 
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Thermal expansion of KCl. 
—_— Theoretical ; @-@ Experiment (Eucken and Dannohl 1934). 


it is of the same general shape as that obtained for C, on the Debye 
theory as might have been expected. Moreover, it exhibits an upward 
curvature at high temperatures which does suggest that the Debye and 
Griineisen theories as applied here will account for this phenomenon, 


« 


646 G. C. Fletcher on the 


which is in fact observed, without postulating the presence of lattice 
defects. However, there is less cause for satisfaction when the curve is 
compared with experiment (Eucken and Dannohl 1934) as is obvious 
from fig. 2 and it is worth while considering various factors which may 
affect the theoretical curve at least for medium and high temperatures. 
If one puts J’=(n—m)J and y=1+z in the left-hand side of (20) and 
expands by the binomial theorem, it takes the form 


|e  48-m--o > 
fg ae.|, 


For the right-hand side ¢’—1 was expressed as a series and integration 
carried out term by term, giving a series valid for 7>Ky ('»’S, Le. 
7 +350°K in the present case. Again putting y=1+2 and expanding 
by the binomial theorem eqn. (20) may finally be reduced to a polynomial 
in z, whose coefficients are power series in Q/7’. In the present case it 
was found possible for 7>2K=700°K to neglect most of the terms 
resulting in the high temperature approximation 
a WY ea ei ecg Gd a x 
ae Ve T 6 eelirery . . . . (23) 


and 
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Small changes in the individual values assumed for m and 7 are unlikely 
to affect the sum n+m involved in this and in 

dn =6—ntm KS 

eee ee 

oT By) es 
The other important factor involved is K%/.J’ and from (21) this is found 
to be equal to N(m-+2)koh/2Vo. This does not depend critically on », 
which is certainly quite large, nor is there likely to be much error in 
the value assumed for V5. However, judging by Durand’s own estimates 
Of (C44): (Cr2)9 and Kp at 7’=303°K, it does appear that the value assumed 
here for x) is probably too low. A larger value would increase both « 
and dz/07' at high temperatures, improving the agreement with experiment. 


§ 6. PHystcaL NATURE OF ASSUMPTIONS 

Some consideration will now be given to the physical meaning of the 
assumptions made by Griineisen in deriving (6) and (10) which are also 
involved in the derivation of (20). 

(a) If the volume of a solid is kept fixed and the temperature varied, 
some of the heat added or taken away will change the energy and therefore 
the distribution of the electrons in the lattice. This will result in an 
alteration of the interatomic potential and hence of the total potential 
energy & of the solid. Thus in general @ will not be a function of volume 
only although this may be regarded as a valid approximation for certain 
solids at temperatures where electronic effects are generally regarded 
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as negligible. When it is necessary to consider the electrons specifically, 
e.g. in metals at low temperatures, it may be preferable to abandon 
the concept of a separate potential energy © and treat the total free 
energy F as a single quantity either by classical or quantum methods. 


(6) The total internal energy of a solid is defined by 


U=0+h=F418—F—7 (55) 
? 
so that H=7 —T(07 /éT), and by (4) 


O70 i 
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It is very difficult to get a physical picture of y from this definition 
or to appreciate the physical meaning of the assumption that dy/d7'),,=0. 
On the Debye theory (25) reduces to the usual well-known formula and, 
if a wave of frequency v,, has wavelength \,, and velocity c,, this is 


olny olnc olnaA 


™m m m 9f 

abun. oi din Se) 
The minimum wavelength A,, must be of the order of the distance between 
nearest neighbouring atoms in the lattice so that it is reasonable to 
assume A,, oc J'1/3 and the second term then contributes the constant 
value of } to y. As regards the first term, if V is increased, interatomic 
coupling is decreased and therefore c,, is decreased and this term is 
positive, giving y>4. More than this it is difficult to say without 
considering a specific model. The only reasons why y should vary with 
temperature at constant volume are the consequent changes in the 
interatomic forces, already considered above, and changes in the amplitude 
of the atomic vibrations thereby producing a smaller or larger mean force 
on one atom due to its neighbours which would affect c,,,. 

For different types of wave, e.g. longitudinal and transverse, a r,, 
and y may be defined for each; in general they will be different from 
one another and will show different temperature variation as the 
corresponding c¢,, are different functions of c,, and c,, for instance 
(see (22)). Again, 0 In J,,/¢ In V’ should be ~} for waves in all crystallo- 
graphic directions but @ Inc,,/¢ In V and hence y will vary with direction 
in general. Finally for crystals with atomic groups of much stronger 
internal coupling than that between groups the group dimensions will 
vary much less than the crystal dimensions. Thus y's which may be 
defined for vibrations within a group should be smaller than those defined 
for vibrations of the crystal as a whole. 

(c) Quantum theory shows that a simple harmonic oscillator has some 
vibrational energy even at absolute zero, so that the equilibrium con- 
dition there should be (@F/0V),_9=0 and not (d0/dV’), )=0, which 
would only be true if all the atoms were at rest. This point does not 
appear to have been considered before and is probably unimportant 
when relatively crude theories of the thermal behaviour of solids are 


involved. 
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§ 7. Discussion 


It has been shown that two well-known formulae due to Griineisen 
for the thermal expansion of a solid depend essentially on three 
assumptions. These same assumptions have been made in deriving 
a direct relation between volume and temperature, based as closely 
as possible on the theories of Debye and Griineisen. Application to 
the ionic solid KCl has shown that, although the experimental observations 
on the thermal expansion of KCl may be explicable without postulating 
the occurrence of lattice defects at high temperature, Debye—Griineisen 
theory does not give a good general agreement with experiment. This 
contrasts with the conclusion of Fischmeister, whose method, however, 
is shown to be too insensitive to be conclusive. It is realized that the 
present theory is very crude in that the Debye frequency spectrum for 
a solid is not regarded as more than a poor approximation and also that 
the Debye model of an elastic continuum has been mixed with the 
Griineisen atomic model but the author hopes to extend the method to 
more realistic models of solids by use of the lattice theory of Born and 
von Karman and of plasma theory, particularly in the case of metals. 
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ABSTRACT 


The inhibition of plastic deformation by dislocation locking promotes 
brittleness. The influence of manganese in iron is studied in detail 
and it is concluded that the ductility conferred by the manganese is 
associated with weakened dislocation locking as well as with grain refine- 
ment. The weakening of the locking appears to be due to interaction 
between nitrogen and manganese atoms. 


$1. InrrRopucTION 


Ir now seems fairly certain that the fracture of polycrystalline metals 
originates in cracks formed ahead of arrays of dislocations piled up 
against grain boundaries (Zener 1948, Mott 1953, Petch 1953, 1953, 1956 a, 
Stroh 1954, Greenwood and Quarrell 1954). Brittle (cleavage) fracture 
probably occurs when the crack produced by a dislocation array is 
sufficiently large to spread catastrophically across the specimen by its 
own concentration of the applied stress. On the other hand, ductile 
fracture would seem to arise when the crack produced is too small to 
spread in this way ; further deformation, with the formation of further 
cracks, is then necessary to complete the fracture. 

If plastic deformation occurs in the region of high stress ahead of 
a discloation array, cracking may be suppressed. Alternatively, if a 
crack does form, it will be a short one when its formation is accompanied 
by plastic deformation, since some of the energy of the array is then 
dissipated in plastic work. Thus, plastic deformation gives protection 
against fracture, or, if this does occur, ductile fracture is favoured. 
Since dislocation locking influences the occurrence of plastic deformation, 
such locking will also influence the occurrence of fracture. 

If the dislocation locking is sufficiently strong, complete brittleness 
in a tensile test will result. Less strong locking may be wiped out 
by deformation prior to fracture, but even then the onset of plasticity 
will be delayed to a higher stress and the plastic strain at fracture will 
be reduced. 


+ Communicated by the Authors. ; 
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The possibility that the dislocation locking is wiped out by prior 
deformation is less in fracture from a notch, since the strain rate is 
higher. In such fracture, dislocation locking will reduce the amount of 
plasticity associated with the spread of the crack from the notch and so 
will increase the crack velocity. This, combined with the locking ahead 
of the advancing fracture wil! increase the probability that it eventually 
changes to the cleavage type. 

The effect of dislocation locking partly, possibly mainly, explains the 
greater brittleness of the body-centred cubic metals compared to the 
face-centred cubic ones, for the locking is particularly strong in the 
body-centred cubic structure. Enhanced ductility follows removal of the 
locking impurity. Thus, chromium is ductile in the absence of nitrogen. 
Similarly, the low temperature ductility of iron is increased when the 
dislocations are unlocked by a prior small deformation at room temper- 
ature (fig. 1). 

Fig. 1 
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The influence of prestrain at 18°C on the ductility of mild steel at —196°c 
(Patwardhan and Petch, unpublished). 


Technologically, the most important case of brittle fracture is that of 
mild steel. There, one of the principal safeguards against the occurrence 
of this type of fracture lies in the presence of sufficient manganese in the 
steel (Barr and Honeyman 1947). This has never been fully explained. 
Partly it is due to the finer grain size produced by the manganese and to 
the lower carbon content that can be used to give the same ‘ maximum 
stress * but this does not appear to be the complete answer. The present 
paper presents a study of the possible association of this effect of 
manganese with an effect on the dislocation locking. 
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§ 2. MEASUREMENTS ON MANGANESE [RONS 


The approach adopted involved an examination of the effect of Mn on 
the lower yield point o) yp, and the fracture stress o¢ of polycrystalline 
iron at 18 and —196°c for a wide range of grain diameter 1. This 
permitted an evaluation of the influence of Mn on the constants oy. & and 
k* in the relationships 


G1 y p.=%+Al- 1”, Sg ences Men bis (Ly) 


or=0,tk*I-12, 5 St AA a, ee ae 


The form of (1) is consistent with the propagation of a Liiders band 
from grain to grain by the shear stress ahead of an array of dislocations 
held up by the grain boundary and proportional in length to 1. The 
form of (2) is consistent with the production of fracture by the tensile 
stress ahead of such an array. This relationship applies to both the 
cleavage and ductile fracture of ~-iron and k* has practically the same 
value in both cases (Petch 1956 a). 

The magnitudes of k and k* are a measure of the stress required ahead 
of the array for the propagation of the yield or for cracking. Thus, if 
Mn modifies the dislocation locking, it should alter /. 


Table 1. Composition of the Irons, Weight °, 


Mn © Si IP Ss N Al 
0-02 0-04 0-12 0-03 0-01 0-006 0-03 
0-25 0-03 0-12 0-03 0-01 0-006 0-03 
0-47 0-04 0-11 0-03 0-01 0-006 0-02 
0-91 0-04 0-14 0-03 0-01 0-006 0-02 
1-90 0-05 0-17 0:03 0-01 0-008 0:02 


| 
Ni, Cr, Mo and V all <0-02. 


The constant co, is a measure of the stress opposing the motion of a 
free dislocation through the grain and consequently opposing the compres- 
sion of the dislocation array (Cracknell and Petch 1955, Heslop and 
Petch 1956). 

Five irons were examined in which as far as possible the only composi- 
tion variation was in the Mn content. 

The general experimental procedure for the determination of the 
yield point and fracture stresses was the same as previously employed 
(Cracknell and Petch 1955, Petch 1956a). The variation in erain size 
was obtained by the same types of heat treatment as before and the 
specimens were all tested in the annealed condition. Tensile specimens 
with lin. gauge length of 0-125 in diameter were used and the testing 
machine was driven at a cross-head velocity of 0-013 in./min. As 
previously, the fracture stresses were corrected to zero plastic deformation 
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The influence of Mn on the fracture stress and lower yield point of «-iron 
at —196°c. 
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The influence of Mn on the lower yield point of «-iron at 18°c. 
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at —196°c and to the average plastic strain at 18°c (Petch 1956 a, b). 
Correction was also made for departure from uniaxial tension, if necking 
occurred. 

Figure 2 shows the dependence of oj yy and of of upon 1-1/2 for four 
of these irons at ~196°c. Figure 3 contains o, yp. at 18°c, and the 
values of k, k* and oy at 18° and —196°c are given in figs. 4 and 5, 

It will be seen that the initial additions of Mn have no effect. After 
that, o) imereases, the fracture stress slope k* is unaffected and the 
yield stress slope / decreases considerably at —196°c, but the effect is 
small at 18°c. 

The reduction in k at —196°c increases the gap between o) yp» and 
of and a corresponding increase in ductility is observed. For example, 
the reduction in area at fracture was about 20°, with /~! in the range 
6-7 mm~?? for the low manganese irons, but 45%, for the 1-9°% Mn iron. 


§ 3. Tae INFLUENCE OF Mn on k anv k* 

The initial lack of any effect of the Mn is readily explained. The 
Mn combines preferentially with the carbides not in solution in the «-iron, 
and only when these are satisfied is there any appreciable rise in the Mn 
content of the ~-iron itself. 


Table 2. Grain Sizes of the Irons after similar Treatments 


Mn (%) 0-02 0-25 0-47 0-91 1-9 
Grains mm? 980 1350 1300 5450 9200 
ee (mn ") 5:6 6-1 6-0 8-6 oS 


| 


The lack of an effect on the fracture stress slope k* is also to be 
expected, since k* depends on the rigidity modulus and the surface 
energy and neither of these should be much affected by small amounts of 
Mn. With some irons, manganese raises the fracture stress and increases 
the ductility by the prevention of intergranular brittleness, but there is 
no sign of this in the present case. 

The really interesting point is the alteration in the yield stress slope /. 
It is clear that at —196°c manganese makes the propagation of yield 
from grain to grain at the lower yield point easier. This might arise 
from an effect of Mn on the distribution of sources ahead of a pile-up 
of dislocations, but there seems no reason for expecting this. The more 
probable explanation of the reduction in & is a weakening of the dislocation 
locking. This then would seem to be the cause of the increased ductility 
produced by manganese at —196°c and constant grain size. 

In an impact test on notched specimens, the high strain rate in the 
notch makes the test near room temperature comparable to a low temper- 
ature tensile test. Thus, a lowering of the brittle-ductile transition 
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temperature at constant grain size by the Mn was expected and was 
observed in measurements using Charpy-type test pieces of 5/16 in. 
diameter (fig. 6). | 

The beneficial effect of Mn on the ductility at constant grain size is 
reinforced by the grain refinement that is also produced. This refinement 
is illustrated in table 2, which gives the grain sizes obtained in these 
irons after the initial ingots had received similar treatments in forging 
and rolling down to bars and had finally been held at 900°c for two hours 
and furnace cooled. 
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The influence of Mn and of grain diameter | on the transition temperature. 


$4. THE INFLUENCE OF Mn ON THE DISLOCATION LOCKING 


Normally, the locking in an annealed iron should be mainly due to 
nitrogen, since the concentration of this in solution is greater than that 
of carbon during the temperature range where a condensed atmosphere 
exists. 

One way in which Mn may weaken nitrogen locking is by displacement 
of some of the nitrogen atoms from the interstices they occupy at a 
dislocation in a plain iron to neighbouring interstices. Damping measure- 
ments on iron—manganese alloys show that the free energy of a nitrogen 
atom is lowered by 0-12 ev when the metal atoms that surround it 
include a Mn one (Dijkstra and Sladek 1953). This binding of carbon 
to an Fe—Mn interstice is sufficiently great compared with the 0-5 ev 
binding to a dislocation for some scatter of the nitrogen atoms about 
their normal positions at a dislocation to result, with some consequent 
weakening of the dislocation locking. 

The extent of the weakening that could arise in this way is uncertain, 
because of the uncertainties in the calculation of the interaction energy 
between nitrogen and a dislocation. However, fig. 7 can be used to 
illustrate what may happen. Figure 7 shows the interaction energy for 
a nitrogen atom in the various octahedral interstices below a positive 
edge dislocation lying along [211] in a (011) plane with the Burgers 
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vector in the [111] direction. The equilibrium position is taken to be 
as shown and the calculation follows Cochardt e¢ al. (1955). The letters 
a, b, c, indicate the directions of the tetragonal distortion of the interstices, 
and the interstices where one of the metal atoms in the direction of 
tetragonal distortion lies within the core of the dislocation have been 


neglected. 
Fig. 7 
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The interaction energy in ev for a nitrogen atom in the octahedral interstices 
below a positive edge dislocation in «-Fe. 


According to this calculation, the A interstice with extension along 6 
has the greatest interaction energy and the N atoms should normally 
condense there. But, there are a number of other interstices with 
interaction energies greater than ~0-45 ev and the N will prefer these if 
they are Fe/Mn interstices. 
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Consideration of a model shows that there are 16 metal sites per A 
interstice, occupation of any one of which by a Mn atom will cause 
displacement of the N atom from A. To a first approximation, the 
resultant contribution of the displaced atoms to the dislocation locking 
will be zero, so that, with Mn present at a concentration of ¢°%, the 
strength of the locking will be reduced by a factor {1—0-16 (c—0-25)} 
when 0-259, Mn is removed by the carbon. [If the slope & is directly 
proportional to the strength of the dislocation locking, / will be reduced 
by the same factor and this gives a 28° reduction at 2% Mn. Thus, it 
appears that an appreciable effect on / could arise in this way, although 
this conclusion may well be wrong if fig. 7 is not a good Spey 
to the interaction energies. 

Another possible explanation of the weakening of the locking by Mn is 
that nitrogen is removed from the dislocations and replaced by carbon 
and that carbon locking is weaker than nitrogen locking. 

Manganese lowers the solubility of-nitrogen in «Fe (Dijkstra and 
Sladek 1953) and the amount left in solution might become too low for 
dislocation locking. This does not seem to be so, however, because the 

9°, Mn iron strain-aged at 35°c at about the same speed as the 0-25°%, Mn 
iron and strain-ageing at 35°C is mainly due to nitrogen. Another 
possibility of the replacement of nitrogen at the dislocations arises 
because carbon does not appear to have any preference for the Fe—-Mn 
interstices (Fast 1950). If carbon and nitrogen are both present, the 
lowest energy state may be one with the carbon atoms at the dislocations 
and the nitrogen atoms in the Fe—Mn interstices away from the disloca- 
tions. 

Thus, replacement of nitrogen locking by carbon locking seems possible, 
but the suggestion that the latter is weaker has not previously been made. 
However, it has been found in concurrent work that deoxidation of a steel 
with 0-07% Al also reduces k at —196°c (Codd and Petch, unpublished). 
Now Al combines with nitrogen as well as oxygen and, with this small 
amount of Al. replacement of nitrogen locking by carbon locking seems 
more probable than scatter of the N atoms about their normal locking 
positions. This suggests that carbon locking is im fact weaker than 
nitrogen locking. 

Experiments to distinguish between these two possible mechanisms of 
lowering k by Mn have not yet produced a definite result. 


§ 5. THe INFLUENCE oF Mn ON oy 


Cracknell and Petch (1955) have obtained an expression for the 
resistance to dislocation motion due to interaction with random solute 


atoms. This gives 


T aes / 
yu ES 1 io is 


o.= —— = | —_ ln =; 
0 (a) | 3 Ar NLP 
where N is the number of solute atoms per unit volume, 2~0-5 is a 
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constant, js is the rigidity modulus, is the Burgers vector, A is the 
constant in Cottrell’s expression for the interaction energy and 7ry~0. 

Treating the resistance due to random Mn, C and N atoms as inde- 
pendent to a first approximation, (3) can be used to estimate o, due to the 
Mn. 

From the x-ray data of Walters and Wells (1935), the degree of misfit « 
of manganese in iron is 0-02, so A is ~3x 10°?! dynes cm? and ay is 
~1 ton. in.? for 1-6°% Mn in the «-Fe. Thus, (3) underestimates og due 
to Mn in the same way that it underestimates o, due to interaction with C 
and N (Cracknell and Petch 1955). 

The temperature dependence of co, due to Mn can be estimated as for 
carbon (Heslop and Petch 1956). If 7; is the shear stress required for 
dislocation motion without thermal help and 7 is the actual value required 
at room temperature, this calculation gives t=0-37;, so, in contrast 
to carbon, some temperature dependence of co, due to Mn would be 
expected. Figure 5 suggests that a little temperature dependence 
does exist. 
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ABSTRACT 

Magnetization measurements on the hexagonal and the face-centred 
cubic phases in several ferromagnetic Co-rich alloys in these systems have 
been made. By the use of suitable extrapolation procedures. the magnetic 
moments of both phases in each system at the absolute zero have been 
estimated. In the CoMn alloys there is a significant difference in the 
dependence of moment on composition between the two phases. This is 
not the case in the CoCr alloys, where the rate of change of moment with 
composition of both phases is the same. 


§ 1. InrrRopuctTIon 


In pure cobalt and in a number of cobalt-rich binary alloys, the equili- 
brium phase at low temperatures has a hexagonal close-packed structure, 
while at higher temperatures a face-centred cubic phase is stable. Both 
phases are ferromagnetic. Previous work (Myers and Sucksmith 1951, 
Crangle 1955) has shown that in pure cobalt, a definite difference exists 
between the magnetic moments of the two phases when extrapolated to the 
absolute zero of temperature. With this difference in view, the present 
work was undertaken to establish whether any difference in magnetic 
behaviour on alloying occurs between the two phases. 

The two alloy systems investigated here, cobalt—manganese and cobalt— 
chromium, are of special interest because of the proximity of both solute 
elements to cobalt within the first transition series. Also a certain amount 
of information on the hexagonal phase was already available (Sadron 
1932, Farcas 1937), and their metallurgical properties seemed suitable. 

The fact that a long extrapolation is necessary in order to estimate the 
magnetization corresponding to the absolute zero from measurements made 
at the lowest temperatures at which the cubic phase is stable has hitherto 
deterred precision work on this phase in cobalt-rich alloys. The only 
previous measurements on the phase were in pure cobalt, depending on the 
reduced magnetization~temperature curve having the same shape over its 
available temperature range as that of nickel, on which low temperature 
data exist ; and on some cobalt-copper alloys, depending on _ the 
magnetization at zero temperature varying linearly with composition. 
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Preliminary observations showed that neither of these conditions 
obtained in the cobalt-manganese system, and that the method previously 
used for pure cobalt was not applicable to the cobalt-chromium system. 
A new extrapolation technique was therefore devised, which is applicable 
in favourable cases. 


§ 2. MpraLLURGICAL PROPERTIES OF THE ALLOY SYSTEMS 


2.1. The Cobalt-Manganese System 


According to Hashimoto (1937), the solid solution of manganese in 
cobalt occurs over a wide range of compositions, extending to at least 50% 
at room temperature, and to not less than 50% at higher temperatures. 
The temperature at which the hexagonal phase transforms to the face- 
centred cubic on heating falls with increasing manganese content, from 
about 400°c in pure cobalt to 0°c at about 25°, of manganese. The Curie 
temperature of the cubic phase is shown by Hashimoto to fall from over 
1100°c in pure cobalt to room temperature at rather less than 40°, of 
manganese. 

2.2. The Cobalt-Chromium System 


Chromium is soluble in cobalt up to more than 30%. 

In the phase diagram reproduced by Hansen (1936), the temperature of 
the hexagonal to face-centred cubic phase transition increases with 
increasing chromium content. This diagram also indicates that the Curie 
point decreases uniformly with increasing chromium, so that alloys 
containing more than about 20°, of chromium are not ferromagnetic at 
and above room temperature. Thus the temperature range over which a 
ferromagnetic face-centred cubic phase is stable decreases rapidly with 
increasing chromium content. At 6% of chromium the cubic phase only 
exists at temperatures above about ihren quarters of the absolute Curie 
temperature, compared with pure cobalt, where the phase change occurs 
at a temperature of about half the Curie temperature. 


§ 3. PREPARATION AND CHEMICAL ANALYSIS OF THE ALLOYS 

Alloys having various compositions in the cobalt—manganese and cobalt— 
chromium systems were made by melting together appropriate amounts 
of cobalt and the respective solute elements in an induction furnace. 
Each melt weighed about 20g. The metals from which the alloys were 
prepared were of high purity and the crucibles were of pure recrystallized 
alumina. ‘The melting was carried out in vacuo, or in cases where serious 
metal loss by evaporation was expected, in pure argon at reduced pressure. 
To ensure adequate mixing, alloys were kept molten for up to ten minutes, 
and also each alloy was subsequently homogenized by heating at a tem- 
perature just below its melting temperature for about half-an-hour 
immediately after solidification. 

All the alloys were analysed chemically for the amount of the solute 
elements, and also spectrographically for the presence of impurities. 
The analyses are shown in table 1. 
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hatte S$ 4. EXPERIMENTAL Mrruops 


Measurements of magnetization in strong magnetic fields were made 
using experimental methods which have now become standard in this 
Jaboratory, and the established extrapolation techniques (covering the 
‘variation of magnetization with applied field) were used to obtain values 
of the spontaneous magnetizations of the alloys at temperatures up to 
their respective Curie points (see Crangle 1955). About nine different 
applied field strengths were generally used, between about 3 and 15 
kilo-oersteds. The | Curie temperatures were estimated by plotting graphs 
of the squares of the spontaneous magnetizations against temperature, 
‘and extrapolating the parts with steepest slope to the temperature axis. 
‘The probable error involved in the various Curie point determinations ‘is 
‘less than +2°c. 


Table 1. Analyses 


Cobalt—Manganese + Cobalt—Chromium { 
% Mn by weight} 1-87 | 3-58 | 4:91 9-87 | 17-48 == = — 
% Cr by weight | — -= 2-62 | 4:44 | 5-20 
°% Mn atomic 2-00 | 3:83 | 5-25 | 10-51 | 18-52 — — — 
°% Cr atomic = = — — — 2-96 | 5:00 | 5-85 


+ Spectrographic traces (<0-01°%) were detected of Cu, Ni, Fe, and Si. 
t Spectrographic traces (<0-01°%) were detected of Cu, Fe, and Si. 


In order to estimate the spontaneous magnetization of the face-centred 
cubic phase at the absolute zero (c,), an extrapolation procedure has to 
be followed, since the phase is not stable at low temperatures. The new 
procedure which has been developed is as follows. 

Measurements on an alloy of composition c give values of the spontaneous 
magnetization o, at various temperatures 7’. When the Curie temperature 
6 of the alloy has been measured, the reduced temperature 7'/@ corres- 
ponding to each o,7 is thus known. 

We may write 


Op—% (O7/Fo)c. T/0° 


Since o, is unknown, (7/09), 7/6 Will in general be unknown also. But 
it is related to the reduced magnetization of the pure solute element at 
the same reduced temperature by 


f 


(a7/Fo) T/6— (o7/o)o, ro 4 


_ where J expresses the deviation from the Law of Corresponding States at 
_the particular reduced temperature and composition being considered. 

A increases with increasing 7'/@ in the lower temperature range : that 
is, 4 increases as (07/09), 7¢ deviates from unity. Since the reduced 
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magnetization only falls gradually with increasing temperature when 
7/6 is small, it is assumed that it is possible to write 

A=« {1—(ay/09)o} +8 (—(er/90)0°+ - + - 
where « and f are constants, and to consider only the first term of this 
expansion when {1—(c7/09)9} is small. 

Thus we have that 

LGA L(o7/Fo)o, ryt +4 § =a]. 

It is thus to be expected that when (cp/c9)9 7/ is close enough to unity, 
o, will vary linearly with (c7/o9)9, 7,. By plotting 7, for the alloy as a 
function of the reduced magnetization of the pure solute element at the 
same reduced temperature as that to which o, applies, it should then 
be possible to estimate the true co, for the alloy by extrapolation of o, to 
its value where (¢7—99)o, 7/9 becomes unity. 


Table 2. Reduced Magnetization Data for 
Face-centred Cubie Cobalt 


T/@ op/ Go T/é O71 Fy 
0-493 0-934 0-949 0-467 
0-508 0-930 0-961 0-418 
0-536 0-919 0-971 0-385 
0-579 0-899 0-977 0-344 
0-624 0-879 0-980 0-331 
0-647 0-867 0-985 0-295 
0-673 0-849 0-987 0-271 
0-699 0-831 0-990 0-241 
0-753 0-787 0-992 0-211 
0-795 0-747 0-993 0-193 
0-836 0-695 0-996 O05 
0-873 0-641 0-999 0-120 
0-895 0-602 
0-911 0-571 
0-932 0-518 


Measurements on the magnetization of pure face-centred cubic cobalt 
exist down to about 400°c (7'/@=0-5) (Crangle 1955), and reduced data from 
them are given in table 2. Intermediate values of (op/o9)) were obtained 
for the various 7'/@ by graphical interpolation, with extrapolation to 
T'/6=0-4. The known magnetization data for pure nickel, which follows 
the same reduced curve as cobalt (Myers and Sucksmith 1951) were used 
as a guide to this extrapolation. 

Graphs of o as a function of (o7/o)y for all the present series of alloys 
are shown in fig. 1. The quality of the straight lines obtained over the 
range shown justifies the procedure used. In the case of the 2-96°% CoCr 
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alloy, the onset of the phase change from hexagonal to cubic on cooling 
was detected solely by the appearance of a definite departure from the 
main line which characterizes the cubic phase, at a temperature agreeing 
well with that determined previously by direct methods. In all the other 


Fig. 1 
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Extrapolation to the magnetization at the absolute zero for the face-centred 
cubic phase. The start of the phase transformation in the 2-96%, CoCr 
alloy is marked. 


alloys, the phase transformations were detectable also from marked 
discontinuities in the magnetization—temperature curves. 

In the case of the hexagonal phases of both alloy systems, considerable 
magnetic hardness was found, and in order to be able to estimate intrinsie 
magnetizations, extrapolations to saturation conditions were necessary. 
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Weiss and Forrer (1929) used a method which involved the extrapolation 
to infinite field of a graph of magnetization against the reciprocal of the 
field strength. Néel (1948) gave a theoretical justification for plotting 
against 1/H when the magnetic hardness is caused by cavities or other 
non-magnetic imperfections in the specimen, and against 1/H? when the 
hardness arises from magnetocrystalline anisotropy. 

During the present work on the hexagonal phases, applied fields ranged 
between about 7 and 15 kilo-oersteds. It was found that the CoMn 
alloys gave better straight lines when 1/// was used, and in the case of 
CoCr, 1/H? was better. The uncertainty in which parameter is correct 
affects the final estimates of spontaneous magnetization at the absolute 
zero. Values from the 1/H? extrapolations are of the order of 1° lower 
than when 1/H is used. At temperatures as far below the Curie point as 
those at which the hexagonal phase is stable in these alloys, any error 
arising out of the assumption that the intrinsic magnetization is effectively 
identical with the spontaneous magnetization will be insignificant com- 
pared with the uncertainty in extrapolating to saturation conditions. 

Values for the spontaneous magnetization of the various alloys in the 
hexagonal phases were estimated by making linear extrapolations of 
graphs of spontaneous magnetization against the square of the absolute 
temperature (fig. 2). 
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Extrapolation to the magnetization at the absolute zero for the hexagonal phase. 


§ 5. RESULTS AND Discussion 

It may be seen from fig. | that there is no difficulty in carrying out the 
‘extrapolations to o, for the face-centred cubic cobalt-manganese alloys 
with an error of well under -+1-0 ¢.g.s. units. For the cobalt—chromium 
alloys. however, the error is rather greater, due to the longer extrapolations. 

In the case of the 18-52%, CoMn alloy, the magnetization measured at 
any given temperature up to about 540°C appeared to be somewhat 
‘sensitive to heat treatment, depending on the rate of cooling. These 
effects are possibly due to an order—disorder reaction and further work 
is in progress to elucidate this point. The data presented here for this 
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alloy refer to the state produced by cooling quickly. No dependence of 
magnetization on heat treatment was found in any of the four alloys of 
lower manganese content, or in any of the cobalt-chromium alloys. 

: The results obtained for the hexagonal phase of the cobalt-manganese 
system agree well with those obtained by Sadron (1932) up to at least 10° 
of. manganese, as may be seen from fig. 3. No data are plotted here for 
the 18-52% alloy because it was not possible to obtain the hexagonal phase 
in this alloy without some retained cubic phase being also present and 
the magnetization measurements were therefore displaced towards 
higher values. Since Sadron’s data show no sign of this effect, it seems 
_ likely that there was some metallurgical difference between his alloys and 
the present ones, possibly of grain size, which causes different behaviour 
at the change of phase. 

For the hexagonal phase of the CoCr system, there is agreement on the 
rate of change of moment with composition between the present measure- 
ments and those of Farcas (1937). Where composition overlap exists, the 
values of Farcas are apparently about 0-05 Bohr magnetons higher than 
those made here. However, this difference may not be real, largely 
because of the possibility of errors inherent in the method used by Farcas 
for estimating oj. Fareas apparently only made measurements at room 
temperature, and obtained values of co, by multiplying the saturation 
magnetization at room temperature by a factor estimated from knowledge 
of the respective Curie temperatures of the alloys. 


Table 3. Collected Data for CoMn and CoCr alloys 


CoMn CoCr 
Composition, % atomic 2-00 | 3:83 | 5:25 | 10-50 | 18:52 2:96 | 5:00 | 5°85 
Curie point (@), °K 1350 | 1295 | 1243 1100 836 1275 | 1198 | 1160 
Magnetic moment per atom in o 
,cubie phase (Bohr magnetons) | 1-69, | 1-64, OO es Os 0-98, 1-53 «|. b:42. oy 136 


Magnetic moment per atom in 
‘hexagonal phase (Bohr 


_ nagnetons) L625 {oO ,) lta 1-24, | 0-924 1:53 1:38 | 1:34 


+ Probably anomalously high because of the presence of some retained cubic phase. 


The numerical data for all the alloys are collected in table 3, and the 
magnetizations at the absolute zero (after conversion to Bohr mag- 
netons+) are plotted against composition in fig. 3. 

This work clearly shows that the rate of change of magnetic moment 
with alloy composition in some composition ranges of the cobalt—manganese 
system is sensitive to crystal structure, but that no equivalent effect occurs 
in the cobalt-chromium system over the composition range over which 
observations on both phases are possible. 


+The Bohr magneton was assumed to correspond to 5586 ¢.g.s. units per 
gram molecule. — 
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The theoretical treatments of the electronic properties of transition metal 
alloys which have been made in the past cover a wide range of diversity, 
and the conclusions to be drawn from the present observations depend 
strongly on which picture is accepted. The opposite extremes are @ 
simplified band theory, and theories based on the assumption that for 3d 
electrons in elements of the first transition series, a Heitler-London 
approximation is more appropriate. 


Fig. 3 
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Magnetic moments of the alloys plotted as a function of composition. For the 
hexagonal CoMn alloys, the triangles signify Sadron’s measurements, 
and the squares the present data. 


In the simplest application of the band theory to these alloys, it is 
assumed that the 3d and 4s bands are common to all atoms present, and 
that electrons are contributed uniformly by both solvent and solute 
elements. The magnetic properties of any given alloy are supposed to be 
determined by the distribution in the 3d band between states of opposite 
spin of the electrons available. | 

On this theory, the dependence of rate of change of magnetic moment 
with composition on crystal structure would indicate that different 
schemes of electron distribution are stable states in the two phases. 
That is, that the energy density of states curve for the 3d band in cobalt— 
manganese alloys is structure-sensitive. 

It would be expected from this simple theory, especially since manganese 
and chromium are neighbours in the periodic table, that in the systems 
CoMn and CoCr, alloys with the same total number of 3d-+-4s electrons 
would possess the same magnetic moment. This does in fact occur in the 
hexagonal phases of the systems. In the cubic phases, however, there is 


no agreement. No simple explanation of this disagreement appears 
possible within this framework. 
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In the kind of theory which is at the opposite extreme, it is assumed that 
magnetic moments are localized characteristically on the various kinds 
of atom present in a given alloy. The magnetization of the whole is 
made up by compounding the effects of all the atoms. In the simplest 
application, the effective magnetic moments of the atoms of the solute 
elements may be estimated by extrapolating the graph of alloy moment 
plotted against composition, to 100% of solute. On this basis the 
folowing magnetic moments per atom (expressed in Bohr magnetons) 
may be obtained : 


Manganese in face-centred cubic cobalt 
0-5°% Mn — 0:25-0:25 
Above 5°% Mn —2-9 +0-1 


Manganese in hexagonal cobalt 


0-16% Mn = 2-8-40-2 

Above 16% Mn (Sadron’s data) —4-2. 
Chromium in face-centred cubic cobalt —5-1-L0-5 
Chromium in hexagonal cobalt —4:7+0-3 


The errors quoted here are outside estimates, based on maximum errors 
in the magneton numbers of +0-01 in cubic CoMn alloys, --0-02 in all 
alloys in the hexagonal phase, and -+-0-03 in cubic CoCr ; the lengths of 
the various extrapolations are taken into account. 

The negative signs imply that the solute atoms are aligned in anti- 
ferromagnetic fashion with respect to the magnetization of the solvent. 

In the case of the CoMn alloys, if the present simple picture is accepted, 
the variety of the extrapolated moments given here must mean that the 
effective moment of the manganese atoms depends on their detailed environ- 
ment. In the face-centred cubic structure, a given atom has twelve 
nearest neighbours and six next nearest neighbours at a distance 1-4 times 
further away. Thus the composition at which, on the average, there is 
one other manganese atom included in the eighteen nearest or next 
nearest neighbours of a given manganese atom is about 53%. This is 
close to the composition at which the bend in the curve of moment against 
composition for the cubic phase occurs. The presence of Mn—Mn inter- 
actions may then be responsible for a change in the effective moment of 
manganese. This argument is not very convincing, however, for in the 
hexagonal phase the bend found by Sadron was near 15°, of manganese, 
where about one in six of the atoms surrounding a given manganese atom 
is another manganese atom. In the unit cell of hexagonal cobalt, the 
deviation from ideal closest packing is slight (c/a is actually 1-624, compared 
with 1-633 ideally). The nearest six and the six next nearest neighbours 
are all at very nearly the same distance away, all twelve forming a 
configuration very similar to that of the twelve nearest in the face- 
centred cubic structure. It is difficult to believe that a Mn—Mn pair made 


668 Magnetization of Cobalt-Manganese and Cobalt-Chromium Alloys 


up of atoms which are nearest neighbours would have properties very 
different from one made up of next nearest neighbours. The arrange- 
ment of the six further neighbours in the hexagonal structure is quite 
different from that of the next nearest in the face-centred cubic, and they 
are now as much as 1-7 times further away than the nearest neighbours. 
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SUMMARY 


Stacking faults on the (111) planes of several face-centred cubic metals 
and alloys have been introduced by cold work, and estimates of the 
stacking fault probability ~, have been obtained from changes produced 
in the Debye-Scherrer spectrum. The faulting probability increases on 
alloying, from one plane in 300 in copper, to one plane in 25 for some high 
solute content alloys containing zinc, aluminium, tin or germanium. Both 
neutron irradiation (5x 101% n.v.t.) and ‘ quenched-in’ vacancies have 
little significant effect on the faulting parameter. 

Lowering the temperature of deformation increases the faulting 
probability but in copper the faults ‘ anneal-out ’ at room temperature 
after several hours. Line broadening analysis shows that the dislocation 
density is increased on alloying and also by lowering the deformation 
temperature. It is suggested that the faulting may be accounted for by 
the regions between separated half dislocations. Ribbons of stacking 
faults as given by extended dislocations and ‘infinite’ stacking faults 
produced by splitting of dislocations at high stresses are both considered 
and estimates of fault energy made. 

Some work on segregation of solute atoms to faults is reported with 
particular reference to the alloys silver-gold and copper—aluminium. 


$1. [yTRODUCTION 


RECENT experiments have shown that solute atoms (Smallman 1955) 
and irradiation damage (Smallman 1957) produce similar changes in the 
preferred orientation normally assumed by pure face-centred cubic 
metals. It is thought that these texture changes may be related to 
alterations in the stress-strain characteristics of the pure metal ; alloying, 
irradiation or quenching all resulting in a yield point followed by a low 
rate of work hardening. Overshoot of the active slip plane is also 
observed on alloying by Piercy e¢ al. (1955) and slip line studies indicate 
that it is also likely to occur in irradiated (Blewitt and Jamison 1952) or 
quenched (Maddin and Cottrell 1955) material. 


+ Communicated by the Authors. 
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It has long been realized that stacking faults may play an important 
part in these work hardening phenomena since latent slip plane hardening | 
could be accounted for by the difficulty of the dislocations on the second 
set of slip planes cutting through the faulted material on the first slip 
system (Barrett 1952a,b). Consequently an examination has been 
made of the influence of solute atoms and point defects on the stacking 
fault tendency in copper. 


§ 2. EXPERIMENTAL 
2.1. Method 


Paterson (1952) has shown that stacking faults on the (111) planes of 
a face-centred cubic metal produce a shift of the diffraction peaks as _ 
well as broadening them. It is difficult to obtain information about 
faulting from line broadening measurements alone as they are affected 
by both particle size and strain. However the peak shift is affected 
solely by faulting so that a direct determination of the stacking fault 
probability can be made. The peak shift of any one diffraction line 
(426) expressed in degrees is given by Warren and Warekois (1955) 


A20== tan @*cos* b270-4/3a/ata. 1 Eee eee] 


where ¢ is the angle between the reflecting normal and planes containing 
deformation faults, « the stacking fault probability (assumed small) and 
h,= |h+k+l|. The shift is positive or negative depending on whether 
hz is 3n-++1 or 3n—1 (n is an integer) ; when /,=3n there is no shift. 

For the reflection (111) and (200) the effect of stacking faults is to move 
the (111) peakyt to higher and the (200) peak to lower angles. Substituting 
in eqn. (1) this change becomes for the (111) and (200) reflections 


A (26 999 — 2834,) = — 245 1/3(2 tan Aoo9-+ tan 0,,;)/277. . . (2) 


For metals examined with copper Kz, radiation this approximates to 
A (20g 20 44) a ee 


As pointed out by Warren and Warekois (1955) the Paterson analysis 
assumes that the faults occur independently and on only one set of 
(111) planes. In severely cold worked metals the faults probably occur 
on more than one set of (111) planes. Consequently the significance of 
« obtained from eqn. (3) is not clearly defined but is considered to represent 
the sum of the probabilities for the different sets of active (111) planes. 


2.2. Procedure 
The stacking fault probability was determined from a change in the 
(111) and (200) pair reflections. Pairs from the other reflections (220), 
(311), (222) and (400) were also examined to ascertain that the shift was 


7 Three components of the (111) diffraction peak, (111), (111) and (111) are 
shifted while the fourth component remains unaffected. 


Stacking Faults in Face-Centred Cubic Metals and Alloys 671 


due to faults but the results quoted are from the (111) and (200) separation 
since this gives the most accurately measurable estimate of «~. The 
diffraction pattern was obtained with a Geiger counter diffractometer and 
monochromatic Ka, radiation. To minimize positioning errors the peaks 
were manually recorded on a single run. 

Difficulty was originally experienced in estimating the peak positions 
with any certainty and two methods were unsuccessfully tried before the 
one finally adopted. Using the standard Geiger counter arrangement 
a difficulty arises due to the overlap caused by the incident radiation being 
a mixture of Kx, and Kx, components. The observed diffraction peaks 
had to be resolved graphically and the resulting values although indicating 
faulting were insufficiently accurate to give reliable quantitative data. 
A second method was to determine the position of the centre of gravity of 
the diffraction peak shape and use this as a measure of the peak position. 
The uncertainty in this method is associated with the long tails to the 
diffraction lines, especially in cold worked samples, which profoundly 
influence the value of the centre of gravity. . 

The method eventually adopted depended on the elimination of the 
Ka, component of the incident radiation from the diffractometer. This 
was achieved by using a Guinier asymmetrical quartz monochromator set 
to give at its focus resolution of the Kx,«, doublet. It was then possible 
with a sharp slit to cut out the Kx, component, and the resulting diffrac- 
tion pattern is for the copper Kx, wavelength alone. With this experi- 
mental arrangement it was possible to determine the peak of the line to 
+4 minute in 20. 


§ 3. RESULTS 
3.1 
3.1.1. Znfluence of Solute Atoms 


The pure metals copper and nickel and solid solutions of copper with 
nickel, zinc, aluminium, tin and germanium have been examined. 
Spectroscopically pure metals were used (Johnson Matthey and Co. Ltd.) 
and alloys made by melting in vacuo or purified argon. The stacking 
faults were produced by deformation by filing and the stacking fault 
probability obtained by comparing the peak separation for a filed sample 
in both the annealed and cold worked conditions. The results are 
summarized in table 1. Figure 1 shows a plot of the stacking fault 
probability against atomic composition obtained from eqn. 3, for some 
copper alloys. The (26g9)-20,,,) separation decreases as a result of cold 
work and this effect increases in different alloys with increasing solute 
content. The faulting becomes very profuse in highly alloyed materials ; 
for example, in copper the value of «=3-3x 10 corresponds to an’ 
average of 1 stacking fault in every 300 (111) planes, but the addition of 
17 atm.% aluminium increases this frequency to 1 in 25 (111) planes. 
Similarly the addition of zinc, tin or germanium has a significant effect 
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but the effect of an addition of nickel is negligible. The values of the 
copper-zine alloys are in good agreement with those recently published 
by Warren and Warekois. 


3.1.2. Influence of Deformation Temperature 


The influence of deformation temperature on the stacking fault para- 
meter was investigated by examining a sample of pure copper deformed 
at liquid nitrogen temperature. The copper filings produced at —196°c 
were transferred to a pressure die also at liquid nitrogen temperature 
and pressed to a flat disc. This disc was then transferred to the diffracto- 
meter also maintained at —196°c. The measured peak separation 
indicated an increase in faulting probability by a factor of three to four 
compared with a sample deformed at room temperature. The specimen 
was allowed to ‘warm’ to room temperature and the stacking fault 
probability and line breadth measured after various intervals of time. 
The stacking faults completely annealed out together with a major part 
of the line breadths. (See also § 4.) 


Table 1. Values of Stacking Fault Probability from Peak Shift Data for 
Various Copper Solid Solutions 


Deformation Atm.% A (26599-20111) 


Material temp. Solute (minutes) axle 
Copper Room temp. 0 1 3:3 
Copper Liq. nitrogen 0 34-4 11-7+ 13-3 
Nickel Room temp. 0 $ 1-65 
Copper—Nickel | Room temp. 50 13 5 
Copper-zine Room temp. ae ; a 
Coo ee: 5-7 3 34 1011-7 

Se eee Room temp. ; ; : 2 ve 
2:5 14+2 5+ 6-6 
Copper-tin Room temp. 4-5 4} 14:5 
+ 71 10 34 
2-5 13 5 
Copper-— a 
eee ee a Room temp. oe ee - a 
Z 5 (=> oil 


3.1.3. Effect of Annealing on Stacking Faults 


The effect of annealing on the faulting produced by cold work was 
investigated in the copper-17 atm.% aluminium alloy by heating 
specimens for half an hour at successively higher temperatures 100°, 150°, 
200°, 250° and 300°c. The faulting effects annealed out in the range 
200—250°c together with a major portion of the line broadening. 
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3.2. Irradiation and Quench Damage 
3.2.1. Direct Effect of Damage 


It has been suggested by Kunz and Holden (1954) that the interstitial 
atoms produced by neutron irradiation agglomerate into planes producing 
stacking faults parallel to (111) planes. If we assume that the most 
likely site for an interstitial atom in the face-centred cubic lattice is 
that denoted by (3, $, 4) or (4, 0, 0) then an agglomeration would change 
the stacking sequence from the normal ABCABCABC to ABCACBCABC 
in the faulted structure. This is shown diagramatically in fig. 2. 


Fig. 1 


a xX 10° 


10 20 30 40 50 
Atm.% solute 


Stacking fault probability versus composition curves for «-solid solutions of 
copper with nickel, zinc, aluminium, tin and germanium. Full circle 
values for zinc are taken from Warren and Warekois. 


Fig. 2 


Agglomeration of interstitial atoms in the face-centred cubic lattice. 


To examine this point some annealed copper filings were irradiated to a 
total flux of 1019 n.v.t. at an ambient pile temperature of 75°C. Examina- 
tion of the irradiated filings showed no detectable change in peak eels 
tion indicating that any shift is less than } minute of are in 20, There 


2424 
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was, however, a small increase in strain broadening compared with the 
unirradiated metal. ‘Quenched-in’ vacancies also had a_ negligible 
effect on the stacking fault parameter as was shown by examining a 
specimen of annealed copper filings quenched from 1000°c, 


3.2.2. Indirect Effect of Damage 

The effect of neutron irradiation on the stress-strain curve, preferred 
orientation and slip line characteristics, indicates a close analogy between 
irradiation damage and solid solution hardening. To carry this hypo- 
thesis further, samples of copper irradiated to a total flux of 101° n.v.t. 
and 5x 1019 n.v.t. respectively were deformed by filing at room tempera- 
ture. The resulting peak separation for the specimen irradiated to 
101° n.v.t. was the same as for the unirradiated sample filed at room 
temperature. The specimen irradiated to 5x 101° n.v.t. showed a slight 
increase on some measurements but the overall averaged peak separation 
was insignificant (1’+4’). Quenching a sample from 1000°c before filing 
at room temperature also failed to show a change in the stacking fault 
probability outside the experimental error. 


3.3. Line Broadening 


To correlate the changes in stacking fault probability with dislocation 
density the line profiles from the pure metals and alloys have been 
analysed in terms of particle size and strain broadening using the Cauchy 
line broadening relationship discussed by Williamson and Hall (1953). 
If 6, is the total line breadth made up of contributions from particle 
size (<) and strain (€) then, by the Cauchy relationship 


B,=B.+Be 
or B,=/e cos 6+- 2 tan@. -. J+)". Veal. ant4y 


Thus on a plot of 8, cos 6/A against 2 sin 6/A where f, is the correct physical 
breadth in radians, the intercept corresponds to the reciprocal of the 
effective particle size and the slope the effective strain. It is unlikely that 
the values of particle sizes obtained by this method are absolute owing to 
the assumptions involved, but are probably accurate to +50%. However, 
other methods of particle size determination have similar limitations and 
are more time consuming (Williamson and Smallman 1955). Figure 3 
shows a line broadening plot for the copper—aluminium series of alloys. 
The particle size and strain values for these and many other metals and 
alloys are summarized in table 2 

The particle size and strain values in the [100] and [111] directions are 
denoted by «{hkl] and [hkl] respectively ; both change with increasing 
alloy content, the particle size decreasing and the strain increasing. The 
particle size in the high alloy content materials becomes extremely small 
(about 100 A) with the exception of copper-nickel for which significantly 
there is little evidence of stacking faults. One surprising feature of the 
particle size results for the materials containing appreciable faulting, is 
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that the value of the dimension in the [111] direction is at least twice that 
in the [100] direction. The particle shape indicated by such values seems 
most unlikely. Moreover, the absence of such anisotropy in materials 
containing low stacking faults, e.g. copper and copper-nickel indicates 
once again that the effect is due to stacking fault broadening. These 


Fig. 3 


250 


150 


B cos 6/Ax 104 


| Gi) 00) (222) (400) 
0 Seta thal Dat ahs Ve 5 eo ees a 


02 0.4 0.6 0.8 1.0 12 
2 sin 6/AX 10-8 


Line broadening analysis applied to copper—aluminium alloys. 
O 5-7 atm.% Al, x11 atm.% Al, A 17 atm.% Al 


indications are substantiated by observations during recovery (see § 4). 
Thus since stacking faults give rise to broadening similar to particle size 
broadening, the measurement of domain size from line breadths in 
materials containing stacking faults are incorrect. It is possible, how- 
ever, knowing the stacking fault probability « to estimate the contribution 
to the particle size broadening, since Warren and Warekois (1955) show 
the effective particle size due to faults «, , to be 

€g y= 2d41,/3% cos d. ve cee dar (D>) 
From this equation it can be seen that the stacking fault contribution to 
the particle size broadening in the [100] direction is much greater than 
that in the [111] direction. This is partly because ¢ is smaller for the 
[100] directions, and partly because only three out of four (111) planes 
are affected, whereas all (200) planes are affected by faulting. While 
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eqn. (5) refers only to the effective <,, determined from line shape 
experiments, a similar trend may be expected for the effect on the line 
broadening. In the case of copper and copper-nickel filed at room tem- 
perature the contribution to particle size broadening by stacking faults is 
negligible, but the particle sizes observed in the alloys may be com- 
pletely accounted for by fault broadening (see column 8, table 2). For 
example, from the measured values of « for copper—17 atm.% aluminium 
the stacking fault ‘ particle-size’ in the [100] direction, eg, [100] is 61 A, 
while that in the [111] direction, <,, [111], is 145 A. 


Table 2. Particle Size (c) and Strain (€) Values for some of the 
Copper Based Alloys 


Atm 
- | {111} | [200] | €[111]| €[200]] esp [111 
Alloy Treatment | % i k ive i ] a a vt i el ] 
solute 
Filed at : 4 
Copper roomten 0) 1000} 1000) 2 5-5 1665 
Copper Bees IW 450} 200/ 6 | 11 465 
Irradiated = =e 
Copper ‘5 10a t 0 20000|20000} 10-1 | 10-1 — 
Copper-nickel | Flledat | 59 | i000] 1000/ 5:5 | 9 1120 
room temp. 
Copper-—zine + 27 250} 100) 4 10 280 
Copper-tin 3 71 250 EL On ered 13 170 
17:0 170 70| 7 12 145 
pa : 11-0 | 300] 150] 6 | 11 200 
5:7 600 250) 4 6 520 
Copper-— ; 
Germaniin 5 7-2 300 125) 6 10 150 


| 


The strain values obtained are also markedly anisotropic as is expected 
from the differences in Young’s modulus for these directions. However, 
the anisotropy of strain cannot be accounted for by any simple isotropic 
stress model as has often been assumed in the past (Smith and Stickley 
1943, Warren and Averbach 1952), since the strain ratio &[200]/é[111] 
for say copper—27%, zinc is about 2-5 whereas the inverse ratio of Young’s 
modulus is 3-5. This anisotropy is probably related to the non-uniform 
distribution of dislocations. 


§ 4. RECOVERY AND SEGREGATION 


The results in §3 show that pure copper when deformed and examined 
at low temperatures has the faulting parameter of a highly alloyed metal. 
Whether this is due to a decrease in faulting energy or to the amount of 
deformation that can be introduced and retained, was examined from 


Stacking Faults in Face-Centred Cubic Metals and Alloys 677 


measurements of the line breadths and stacking fault probability at liquid 
nitrogen temperature and after intervals of time at room temperature. 


Table 3. Recovery Results for Pure Copper 


Time (mins.) at aa €[111] €[200] 111 £(200) 
room temp. A A 10? x 108 
0 75 450 200 6 ll 

10 85 — — — — 

30 100 — — — — 

60 120 — — — — 

150 150 — — — — 

300 200 — 600 — 4:5 

600 300 — — — = 

5000 > 600 1000 1000 0-6 0-8 


The results in table 3 show that as the stacking faults ‘anneal’ out, so 
the particle size increases suggesting that the particle size broadening is 
due almost entirely to faults. However, the total change in broadening 
cannot be attributed to the removal of stacking faults alone since the 
strain breadth also undergoes a marked change. A remarkable feature 
of the recovery process is that the final breadth is not that exhibited by a 
specimen deformed at room temperature, but has a much smaller strain 
value of about 8x10-4. Figure 4 shows the line shape of the (400) 


Fig. 4 


Intensity— 


se ° 
16° 117° 118° 119 


26— 


Line shape of 400 reflection, showing recovery in pure copper, A filed at —196°c, 
B filed at room temperature ; C filed at —196°c and recovered 3 days 
at room temperature; D annealed at 520°c for 20 minutes. 
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reflection from copper after various stages of recovery. The breadth 
after advanced recovery is intermediate between a room temperature 
deformed specimen and an annealed one. An attempt was made to 
produce this state of recovery in a room temperature deformed sample. 
The closest similarity was shown by a specimen aged at 140°c for 3 hours 
resulting in a similar strain value of 7x 10~* but a larger particle size of 
5 < 10-° cm. 

A similar recovery study was made with copper—17 atm.% aluminium 
and 50/50 atm.% silver-gold alloys. It was noticeable that recovery in 
both-alloys was less marked. The faulting parameter for the 50/50 silver— 
gold alloy is shown in table 4. It is doubtful whether the exact magnitude 
of the increase in probability « from 9-2 to 14-7 with ageing is significant 
as the change was not completely reproducible. However, a decrease in 
parameter such as that in copper was never observed even though slight 
recovery of line shape was detectable. 


Table 4. Recovery at Room Temperature of 50/50 
Silver—Gold Alloy 


Treatment «x 108 a 
Filed at room temperature 2-0 500 
_ Filed at —196°c 9-2 108 


| Filed at —196°c then ‘ annealed ’ 


' Lhr. at room temperature Bah ° 

Filed at —196°c then ‘ annealed ’ ll 
10 hr. at room temperature oy 
Filed at —196°c then ‘ annealed ’ 9.2 108 


17 hr. at room temperature 


It is possible that the effect may be associated with segregation of solute 
atoms to faults, since the solubility in the close packed hexagonal fault 
will be different from that in the cubic matrix, causing the faults to widen, 
as suggested by Cottrell (1953). The electrical resistance measurements 
of Aarts and Jarvis (1954) on this alloy show an increase when the specimen 
deformed at —196°c is allowed to warm to room temperature. If 
faults are responsible for electron scattering, as proposed by Broom 
(1952), the unusual faulting probability values would be consistent with 
this. The possibility of the segregation of solute atoms to faults on low 
temperature annealing was also examined in copper-17 atm.% aluminium 
by this method and by low angle scattering technique, but no change was 
detected. 


§ 5. Estimarus or Srorep ENERGY AND DistocatTion Dgnstry 


Values of stored energy for some of the alloys examined have been 
calculated from the relation 


Vi SiSi2cr, - parties eeleteeeenr 
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where # is Young’s modulus and S? is the mean square strain. The 
mean square strain may be estimated from &, the slope of the line 
broadening analysis, since ss 

eae Seg eat or ni sa (7) 


where A is a constant dependent on the strain distribution and of order 2 
for a Cauchy strain distribution, and 1/27 for a Gaussian. In table 5 an 
intermediate value, A equal to unity has been used. Values of about 
1 cal/g. are obtained which are in good agreement with measurements 
made calorimetrically. 


Table 5. Stored Energy and Dislocation Density Values for some 
of the Materials Examined 
(ek Se sks ake) Se EL ee te, 


; Atm.°% Stored Maximum 
Material ae S Treatment Energy | dislocation 
solute ; 

(cals/g) density p, 
Copper 0 | Filed at room temp. 0-16 pO ROES: 
Copper 0 | Filed at —196°c 0:3 1-5 x 101" 

Copper 0 | Irradiated LelOns 8 x 108 

Copper-nickel 50 | Filedatroom temp. 0-42 Ee 10e 
Copper-zine 27 ‘| Filedatroom temp. 0-51 Ix 10% 
5-7 | Filedatroom temp. 0-21 5x 101! 
Copper—aluminium { 11-0 Filed at room temp. 0-64 1-510" 
(17-0 | Filed at room temp. 0-75 13x10" 
Copper-tin vi! Filed at room temp. 0-83 2-0 x 10!” 
7:2 | Filed at room temp. 0-5 oo 10% 


Copper-germanium 


Dislocation density values have been estimated from the method 
outlined by Williamson and Smallman (1956). This depends on a know- 
ledge of both strain and particle size and consequently cannot be fully 
utilized here because an exact value of the coherent domain size is masked 
by the overwhelming fault broadening. The dislocation density p, 
calculated from the strain broadening is 

ipa eT Doers We earn, aah Ree a nah (es) 

where & is a constant of order 124, b is the Burgers vector and F is a 

strain energy factor. The dislocation density p, calculated from the 

particle size is 

(pe Sh ell 0 seared es pee) 

The inequality of eqns (8) and (9) indicates the probable dislocation 

arrangement and in the case of a piled up array F equals n, the number of 
dislocations in the pile-up, and 

Pe ee fo. ev yonete heirs 5 (10) 


Taking a reasonable estimate for the particle size of about 10~° cm it is 
apparent from eqns. (8) and (9) that the dislocation arrangement is piled 
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up in all cases ; the pile-ups being larger for the alloys and also for pure 
metals on low temperature deformation. 


§ 6. Discussion 


In all the alloys examined the faulting probability is low for dilute 
solutions and increases with increasing solute concentration. It has been 
suggested (Barrett 1952 a, b) that the reason for this is that the difference 
in energy between hexagonal close packing and face-centred cubic packing 
is reduced as the «-phase boundary is approached. This appears to be 
borne out by fig. 1 since the faulting probability increases more rapidly 
the smaller the range of «-solid solution. Difference in atomic size 
between copper and the solute cannot alone account for the variation, 
since germanium although equally effective in producing faults has a 
much smaller misfit value than tin. Arranging the solutes in order of 
increasing valency Ni, Zn, Al, Sn, Ge the atomic misfit values 
{(1/a)da/dc} are 0-025, 0-062, 0-063, 0-28 and 0-062 respectively. A 
correlation with all the factors affecting primary solid solution is 
probably more significant, but tells us little of the origin or form of the 
faults involved. 

It is of interest to note that if the origin of the stacking faults is the 
region bounded by separated half dislocations the dislocation density is 
sufficient to account for the observed values of «. With this assumption 
a represents the fractional areas of the close-packed planes contained 
within these dislocations. Thus if p is the dislocation density andr the 
width of the extended dislocations the stacking fault probability is given 


Nias oeprd: Ba hae nl 


For the case of pure copper deformed at room temperature «=3 x 10-3, 
pP.—3 X 1011 cm and d=2-55 x 10-8 em, then r is about 15 atom spacings. 
The equilibrium separation of partial dislocations is given alternatively 
by the relation (see Cottrell 1953) 

T= [107/24 ry... as asus <p oe ae SL) 
where « is the shear modulus and y the stacking fault energy. Thus 
eliminating r between eqns. (11) and (12) gives a relation connecting the 
dislocation density, faulting probability and faulting energy 


a—=parp/24a/Qry. . 2.7, 1. 6 (13) 
It is difficult to test this formula in practice since measurements are made 
on powders to obtain high deformations and also because line shifts in 
bulk specimens may be due to residual intergranular stresses (Greenough 
1952). However fig. 5 shows « plotted against p for pure copper, deformed 
at room temperature and liquid nitrogen, and for some of the alloys. The 
slopes of the graphs give 
Yoopper=20 ergs CM™*, YEopper-Niexel==35 ergs cm-2 and y Alioye=z 10 ergs cm=?. 
Although these results are in good agreement with measured (Fullman 
1951) and predicted (Seeger 1955) energy values it must be remembered 
that the physical significance of the faulting probability « is uncertain. 
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The non-linearity of the plot for the copper—aluminium series of alloys 
probably indicates the variation of stacking fault energy with alloying. 
Reasonable objections to the above hypothesis may be raised on the 
grounds that other evidence at present available points to the existence of 
wide faults. The oscillation photographs taken by Barrett (1950) of 
supersaturated copper—5% silicon alloys show clearly sharp rods in 
reciprocal space. This indicates the existence of wide stacking faults 


Fig. 5 


a < 108. 


psX 10-11 


Plot of stacking fault probability versus dislocation density for x copper, 
/\ copper-50% nickel, © copper—27°% Zn, (| copper—aluminium. 


rather than ribbons of faults as is expected from extended dislocations. 
As mentioned in § 2, « is derived on the assumption of infinite stacking 
faults and the diffraction pattern from finite stacking fault ribbons is 
not known. Similarly, the electron diffraction experiments of Hirsch 
et al. (1955) on gold foils show the existence of streaks indicating once 
again wide faults. It is possible, however, that the diffraction charac- 
teristics of wide faults may be produced by conditions which are expected 
to occur within a heavily deformed metal. For example, if a pile-up of 
dislocations against a barrier occurs then the addition of many stacking 
fault ribbons could lead effectively to rods in reciprocal space and the 
apparent existence of wide faults. On the other hand wide stacking 
faults may be produced similar to those shown in stainless steel by 
the electron microscope (Hirsch, private communication). The high 
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stresses present in the metals and alloys examined would certainly be 
sufficient to completely separate the partial dislocations. A stress o acting 
normal to the Burger’s vector [110] will repel the two partial dislocations 
completely provided 


o. Do >y; 

that is 
| OS 24/OV/. ae te ee 
For copper, when y=20 ergs cm~? and p.=4 x 1014 dynes cm? the stress o 
is about /160 dynes cm~?. This corresponds to a strain of 6x 10° 
which is less than the values indicated by the line broadening measure- 
ments. It is also significant that a measurable ~ is only found when the 
strains are of the order of magnitude necessary to split dislocations 
completely. In copper annealed for 5000 minutes at room temperature 
the strain is much less than that required by eqn. (14), and in this case « 
is no longer measurable. An order of magnitude estimate of the stacking 
fault energy y is possible from eqn. (14), provided the stress o for which 
aiszeroisknown. Extrapolating to zero « the results for copper deformed 
at liquid nitrogen and also at room temperature, gives £=4x 10~° and 
even this averaged strain yields y about 14ergscm~?. It is difficult to 
obtain more than an order of magnitude of y since the line broadening 
measurements only give averaged strains. It is also possible that some 
splitting of the dislocations might occur near the head of a pile-up (Friedel 
1955) when the average stress is lower than that given by eqn. (14), 
although from the line broadening analysis it appears that in recovered 
copper the dislocation distribution is unpiled. | 

Recovery at room temperature, of a specimen deformed at liquid 
nitrogen temperature, allows a major portion of the broadening to 
“anneal out ’ for the pure metal. The final ‘ recovered ’ breadth is larger 
than that of a well annealed sample but smaller than that of a specimen 
deformed at room temperature. This indicates that the lower the 
deformation temperature the easier it is to recover at a given temperature. 
Moreover, the values of strain 7x 10-4 and_ particle size 10-> cm observed 
indicates that the dislocation distribution is essentially random. This is 
of interest since in copper deformed and examined at room temperature 
the dislocation distribution is piled-up. It therefore appears that the 
barriers causing a pile-up of dislocations at liquid nitrogen temperature 
are dispersed on ‘ annealing ’ at room temperature. This process seems 
quite reasonable since the phenomenon of work softening discussed by 
Cottrell and Stokes (1955) relies on the combined influence of temperature 
and stress to unlock the barriers at the head of piled-up groups. As 
discussed above, with the dislocation density and arrangement known to 
exist in the recovered specimen any faulting would be negligible. 

The inability of neutron irradiation or quenching to produce peak 
shifts directly, tends to indicate the absence of extrinsic faults, since, 
although the Paterson analysis is not worked out for this specific fault; 
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the effect on the diffraction pattern is expected to be similar to that for 
deformation faults. In addition, the low strain broadening introduced 
by irradiation, about 10~¢, indicates that few deformation faults in excess 
of those in the annealed state are produced and are therefore well below 
the limits of detection. Considering the indirect effect of irradiation, 
i.e. irradiation followed by deformation, it is believed that irradiating to 
a total flux of 5x 1019 n.v.t. will produce mechanical behaviour analogous 
to that of an alloy containing 10°% zinc (estimated from preferred orienta- 
tion changes determined by the authors). Moreover, it is unlikely that the 
faulting energy is lowered as with alloying; therefore any increase in 
faulting is on the borderline of detection. 
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ABSTRACT 


6149 complete 7-~-e events have been measured in nuclear emulsion. 
3021 events were in ordinary Ilford G5 emulsion and 3128 were in G5 
emulsion containing twice the usual amount of gelatin relative to halide. 
,The angular correlation between the initial directions of motion of 
u-meson and electron, after possible edge effects have been eliminated, 
is of the form 1—(0:149+0-033) cos@ in ordinary emulsion and 
1—(0-190-L0-033) cos @ in diluted emulsion. This demonstrates 
conclusively the breakdown of parity conservation and charge conjugation 
in both steps of the 7—u-e process. It is shown that, with the aid of 
other data, the partial quenching of the u-meson’s initial spin orientation 
may be eliminated from these data and the true correlation of the 
fundamental process obtained: it is 1—(0-283-+0-060) cos @. This 
asymmetry is consistent with the prediction of the two-component 
neutrino theory with 


(| fa P+) fel) fefat* +h fy *)=0-85+0-18. 


§ 1. INTRODUCTION 


THE abi of Lee and Yang (1956) that parity may not be conserved 
in £-interactions has been tested in a series of experiments, notably those 
of Wu et al. (1957) and Garwin et al. (1957). In particular several groups 
of observers have examined the 7—u-e decay chain using nuclear emulsions 
(Friedman and Telegdi 1957, Castagnoli et al. 1957, Biswas et al. 1957, 
Bhowmik et al. 1957). If parity and charge conjugation are both 
violated in both links of this chain then the angular distribution of 
electrons from y-mesons of spin $ will be of the form 1-+-« cos 8 where @ 
is the angle between the directions of emission of w-meson and electron. 
While indications of violation of parity and charge conjugation appear 
unambiguous in some experiments, the results so far obtained using nuclear 
emulsions are not so certain. 
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A specific objective of our emulsion experiment was to determine the 
coefficient « which measures the asymmetry in order to compare it with 
the predictions of particular models for the 7--e decay. The only 
detailed model so far proposed is that based on the 2-component neutrino 
theory of Salam (1957), Lee and Yang (1957), Landau (1957) and 
Touschek (1957). This model predicts an extreme value of «=—i 
if electrons of all energies are detected without discrimination. It also 
suggests that the asymmetry in the electron emission should be strongly 
dependent upon the energy of the electron. Now while from the experi- 
ment of Garwin et al. one may infer a distribution of the form 1+-A cos 6 
with considerable accuracy the asymmetry coefficient so determined 
cannot be immediately related to the desired « of the fundamental 
process. This is because : firstly, the angular resolution so far achieved 
is poor and one must know how the electron detection sensitivity varies 
as a function of position across the face of the detector before a proper 
folding-in of this resolution function can be achieved ; secondly, the 
sensitivity of the detector varies with the electron energy in a manner 
which is not easy to determine with certainty ; thirdly, the degree of 
polarization of the 4.-meson beam in terms of that for the fundamental 
process isnot known. These difficulties are overcome by using emulsions : 
the angular resolution is very precise ; all electrons are detected without 
regard to energy ; since the complete 7—-e chain is observed it is sure 
that the polarization of the »-meson is as high as is allowed by the 
fundamental process. It appears at first sight, however, that these 
advantages are useless because the asymmetry observed for u-mesons 
stopping in nuclear emulsion is appreciably lower than that observed 
for other substances such as carbon (see Garwin et al. 1957). It is 
therefore apparent that some ‘ quenching’ or loss of memory of the 
initial .-meson spin orientation takes place in emulsion and so it seems 
that measurements in emulsion, despite the advantages listed above, 
cannot yield the true asymmetry coefficient « of the fundamental process. 
We shall show, however, that this is not so, and that the emulsion data 
may be analysed to give the desired true asymmetry coefficient if we 
make use also of the relative apparent asymmetry coefficients observed 
at Columbia for carbon and nuclear emulsion. We do not make use of the 
Columbia asymmetry measurements directly in obtaining our result 
but rather use them to deduce the degree of quenching taking place in 
nuclear emulsion. 

The aim of the present experiment is therefore two-fold : to establish 
clearly the 7-y-e asymmetry in nuclear emulsions ; and to obtain a value 
for the true asymmetry coefficient « integrated over the whole electron 


spectrum. 


§ 2. EXPERIMENTAL METHOD 


The Columbia workers observed a weaker asymmetry in nuclear 
emulsion than in carbon and certain hydrocarbons. They also found 
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that this weaker asymmetry did not depend on the time, after the 
stopping of the »-meson in the emulsion, at which they observed the 
electron. This last observation suggests that some u-mesons stopped 
in emulsion are ‘ quenched ’ very quickly and lose their initial polarization 
while others retain their spin direction for at least some microseconds. 
It is natural to associate the first class of mesons with those that stop in 
one component of the emulsion, probably the halide, and the second 
class with those stopping in the other component, the gelatin. This 
argument suggested to us that it may be possible to observe a stronger 
asymmetry if the content of gelatin in the emulsion were increased 
relative to halide. Such an enhanced asymmetry is obviously highly 
desirable if emulsions are to be used for similar studies of non-conservation 
of parity for rare particles such as the K-mesons. We therefore carried 
out the experiment using both ordinary Ilford G5 emulsions and identical 
emulsions in which the proportion of gelatin to halide is increased 
by a factor two (twice diluted or x 2emulsion). The two sets of emulsions 
were prepared at the same time from the same batches of halide and 
gelatin and they contained the same additives—stabilizers and so on. 
In order to achieve a good grain density on the minimum-ionizing electron 
tracks and also to minimize distortion the emulsions were not stripped 
but were exposed on their glass backing. The ordinary G5 emulsions 
showed a density for fast electrons of 30 blobs per 100 um; the density 
for the <x 2 emulsions was 27 blobs per 100 pm. 

The emulsions were exposed to a flux of about 1047+ mesons per cm? 
at the Liverpool University synchrocyclotron. The incident pion energy 
was about 35 Mev. The emulsions were contained in a steel box with 
a lining of mu-metal. The magnetic field at the position of the emulsions 
was measured with a flipcoil to be 0-15 gauss. This magnetic field is so 
weak that its effect in depolarizing the »-mesons before they decay is 
completely negligible (the precession is by 1:6° in one mean life). 
It is sufficient, however, to effect a very large depolarization of any 
triplet ‘muonium’ atoms (u+-e~) which may be formed (a precession 
of about 160° in one mean life). It is nevertheless likely that such atoms 
would be depolarized very quickly by internal magnetic fields in so 
complicated a structure as nuclear emulsion. 

Before processing, the total thickness, emulsion plus glass backing, 
was measured at several points on each plate using a dial gauge. These 
total thickness measurements were repeated at the same points on the 
processed emulsions and the thickness of the processed emulsion itself 
was measured at the same points under the microscope immediately 
afterwards. By combining these measurements we know the original 
thickness of the unprocessed emulsion at several points on each plate and 
so know the shrinkage factors accurately for the processed emulsions. 
The original thickness of the ordinary G5 emulsion was about 610 um ; 
that of the x2 emulsion was about 480 um (although both nominal 
thicknesses were 600 jm). 
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The emulsions were scanned using air objectives under low magnification 
and 7-1 vertices were sought. The apparent ending of each -megon 
was examined under the same low magnification and if it appeared to 
end in the body of the emulsion the whole event was examined under 
high magnification using oil objectives. It is felt that the p-mesons 
ending more than about 30m from either surface (in the unprocessed 
emulsion) would be detected as ending in the emulsion in the low power 
examination. We consider only complete 7--e events lying within 
a single emulsion and so we are sure of the initial direction of motion 
of the .-meson whose ultimate decay we also observe. | 

The y-mesons frequently suffer strong scattering but because they are 
very non-relativistic their spins are not rotated appreciably by these 
(purely electric) deviations. The angle 6 which we want to measure is 
therefore that between the initial directions of the .-meson and decay 
electron tracks. 

The readings taken for each event were: the angles (read to 1°) of the 
p-meson and electron tracks projected in the plane of the emulsion ; 
and depth measurements (read to 1m) for the 7—y vertex, the point 
along the u-meson track whose projection on the plane of the emulsion 
was 50 um from the z— vertex, the u—-e vertex and the point 504m in 
projected distance along the electron track ; in addition depth measure- 
ments were made for the top and bottom surfaces of the emulsion at 
each p.-e vertex so that the shrinkage factor was continuously monitored 
and we knew the closest approach of each p—e vertex to an emulsion 
surface. Sometimes smaller projected distances than 50 um had to be- 
used for the y.-meson track because of scattering. Smaller projected 
distances were also used for a few very steep electron tracks as a precaution 
against the effects of distortion. 

From these measurements and from the known original unprocessed 
thicknesses cos @ was computed. The computation was programmed 
for EDSAC, the electronic digital computor at Cambridge University. 
Because of our knowledge of the original emulsion thicknesses the 
determination of cos @ is rather accurate and introduces no appreciable 
error into our final analysis of the asymmetry. 

Measurements were made on 3021 complete 7—y-e events in ordinary 
G5 emulsion and on 3128 complete events in the twice diluted emulsion. 


§ 3. ANALYSIS OF THE RESULTS 


These events should not be analysed directly because one may suspect 
a bias in favour of backward-moving electrons for those 4-mesons which 
stop very close to either emulsion surface. This bias has two origins : 
firstly, in many cases the decay electron can be seen even under the low 
power used for scanning and, if the electron is emitted into the body of 
the emulsion, this may draw attention to events in which the ~-meson 
stops so close to one surface that it might have been rejected as having 
left the emulsion had the decay electron been emitted, unnoticed, through 
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the surface of the emulsion ; secondly, when the p-meson is identified 
as stopping in the emulsion, but close to a surface, there is a greater 
chance of failing to find the electron, under oil, if it is emitted out of the 
surface than if emitted (generally backwards) into the body of the 
emulsion. Because of this possibility of bias we analysed the results 
in two groups according to the distance of the -e vertex from the nearer 
surface. The first group, called the ‘centre group ’, is those events for 
which the p-e vertex is further than 50m from either surface of the 
(unprocessed) emulsion. The second group, the ‘surface group’, is the 
remainder. The figure 50m was chosen because it was felt to be a 
generous upper limit to the distance from a surface at which a stopping 
u-meson might. be mistakenly rejected as having left the emulsion in the 
low power scan to find the events. A warning sign that more detailed 
analysis is needed would be that this surface group of events show an 
anomalous apparent preference for backward emission of the electron. 

In the ordinary G5 emulsion the centre group consisted of 2789 complete 
events, while 3 u-mesons which ended in the same region of the emulsion 
did not give decay electrons which could be identified with certainty. 
The corresponding figures for the twice diluted emulsions were 2731 
and 11. These numbers of missed electrons are so small as to have a 
completely negligible effect on the reliability of the results. 

The best way to analyse the results was considered carefully. The 
method chosen was to assume that the distribution was of the form 
1+A cos @ and to compute the mean value of cos 6, which is equal to 
A/3. It may be shown that, when the number J of events is large and 
when A is small, the answer given by this method approximates very 
closely to the maximum likelihood estimate and that the standard 
deviation in A so deduced is +/(3/N). 

The results of this analysis are presented in table 1. 


Table 1. Analysis of the a-u-e Angular Distribution in the form 
1+A cos @. Centre Group means Events whose p-e Vertex is 
more than 50 ym from either Surface of the Unprocessed Emulsion ; 
the Surface Group is the Rest 


Emulsion Group Events A 
G5 centre 2789 —0-149 +-0-033 
G5 surface 232 —0:198+0-11 
2 centre 2731 —0-190 +0-033 
Sc, surface 397 —0-090 +0-087 


We see at once that the asymmetry is established with virtual certainty. 
The second observation is that the events of the surface groups do not 
show a stronger backwards asymmetry than the bulk of the events in the 
corresponding centre groups. This, together with the generous nature 
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of the 50.m thickness of the surface zones} and the fact that we miss 
only ¢% of the electrons from the centre zonet makes us feel confident 
that the values of A displayed in table 1 for the centre groups are not 
biassed. These A values for the centre groups we shall use later, rejecting 
all events in the 50 um surface groups. 

We also note that the asymmetry observed for the twice diluted 
emulsions is a little stronger, though scarcely significantly so, than that 
for the ordinary G5 emulsion. This is at any rate consistent with the idea 
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Angular distribution of the 7—--e decays observed in ordinary G5 emulsion 
(top graph) and in twice-diluted emulsion, containing only half as much 
halide relative to gelatin (bottom graph). The errors shown are the 
square roots of the numbers of events represented by the points and the 
full lines correspond to the asymmetry coefficients displayed in table | 
for distributions of the form 1+-A cos 6, where @ is the angle between 
the initial directions of motion of the p-meson and electron. The 
events shown here are only those of the centre groups for which the 
p-e vertex is at least 504m from either surface of the unprocessed 
emulsion. 


+ Even when the analysis was repeated for surface zones of only 25 wm no 
significantly stronger backward asymmetry emerged. 

+ In a very few cases what appeared to be a decay electron was displaced 
by a few microns from the end of the »-meson track. Scattering of the electron 
through a large angle very close to its origin is most unlikely and it is conceivable 
that we are here witnessing the diffusion of a neutral muonium atom prior 
to its decay. We are informed by Professor G. Bernardini that a similar 
phenomenon has been noticed by workers in Rome. 
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that the quenching of the »-meson spin orientation takes place chiefly 
in the halide. In the figure we display the angular distributions them- 
selves, grouping the events in 0-2 intervals of cos @. The distributions 
are well-behaved and quite consistent with the assumed form of 1+-A cos 0. 

For comparison the results of other workers using emulsions (see § |) 
(in all cases ordinary G5) are gathered in table 2. 


Table 2. Asymmetry Coefficients for 7—u-e Decay as Measured by other 
Workers using Nuclear Emulsions 


Authors Events A 
Bhowmik e¢ al. 1562 —0-:08 +0-05 
Biswas et al. 2003 —0-095 + 0-044 
Castagnoli et al. 1028 —0-222 + 0-067 
Friedman and Telegdi 1300 —0-:12 +0-05 


§ 4. Tor TRUE ASYMMETRY COEFFICIENT 


Our first objective, the demonstration of a clear asymmetry, is reached 
and we now consider how we may infer the true asymmetry coefficient 
of the fundamental z-p-e distribution: 1+acos@. This we are 
obviously unable to do without calling on some other data from which 
we may determine the fraction y of the u-mesons which are quenched in 
the emulsion and lose their polarization before decaying. 

As we remarked earlier the Columbia experiments may not be interpreted 
directly for a variety of reasons; they may be used, however, to give 
a quite accurate indication of the relative initial quenching in carbon 
and G5 nuclear emulsions. If, furthermore, we may assume that there 
is little quenching in carbon we may use the Columbia data to determine 
the quenching y in emulsion, which, again following the Columbia data, 
we take to be rapid compared with the lifetime of the «.-meson. We 
would then have the true asymmetry coefficient «—A/(l—y) where A 
is our apparent coefficient of table 1. The assumption that quenching 
is weak in carbon is likely to be a good one because just the same apparent 
asymmetry is observed for a variety of substances at Columbia and also 
at Liverpool (Cassels et al. 1957) including carbon and copper which 
certainly differ widely in their structure. It would be very surprising 
if the quenching were appreciable but accurately the same in such different 
substances. It is also known (Cassels e¢ al. 1957) that with very good 
accuracy no loss of memory of the u-mesons’ spin orientation in carbon 
takes place over periods as long as 4 microseconds. Since carbon is a 
pure substance the most reasonable conclusion is that the quenching in 
it is small. We make this assumption. 

The Columbia workers measure the emission of electrons in a fixed 
counter as the -mesons precess in a magnetic field. This emission shows 
a succession of maxima and minima as the magnetic field is changed 
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because of the asymmetric emission of the electron relative to the 
-meson’s spin orientation. Call the observed peak-to-valley ratio for 
this emission Re when the j.-mesons are stopped in carbon and Re when 
they are stopped in nuclear emulsion. One may now easily show that, 
irrespective of the bad geometry, lengthy time gate and unknown energy 
response of the detecting equipment, and of the unknown degree of 
polarization of the j.-meson beam, 
Re—1Re.+1 
Seen eel Bei 
e c 

We use the results Re=1:86+0-07 and R.=1:40-+0-07 to deduce 
Vos=0°45-40-11. 

Before considering this result we turn to the problem of determining 
Vx , the quenching in the twice diluted emulsion of the second half of 
the experiment. This we may do, continuing our hypothesis that the 
quenching takes place chiefly in the halide. We do not need to assume, 


however, that all «.-mesons which stop in the halide are quenched. We 
find 


Yx2__ 1—0-86(1—p) 
Yas 2—0-86(2—p) ” 

In the expression 0-86 is the proportion by weight of halide in ordinary 
G5 emulsion and p is the mass stopping power of halide relative to 
gelatin. p may be deduced from the standard datum that the number of 
particles ending their tracks in halide and gelatin in ordinary G5 emulsion 
stand in stopping ratio 7:3. This gives p=0°38 so y, 9/yg;=0:77 and 
Vx g=0°35+40-10. The final result is, however, extremely insensitive 
to this assumed stopping ratio and changes in it over a factor of two 
produce changes in « small compared with errors from other sources. 
We have increased the error in the above value of y,., to allow for this 
uncertainty in the stopping ratio. 

We now use these values of y to correct our observed centre group A 
values of table 1 to find the true asymmetry coefficient «. We display 
the results in table 3. 


Table 3. The Observed Asymmetry Coefficients A, the Quenching 
Factors y, and the Deduced True Asymmetry Coefficients « for 
the Two Types of Emulsion. The Mean Value of « Quoted Allows 
for the Fact that the Error in y is Common to the Two Emulsions 


Emulsion A y a 
G5 —0-149 +0-033 0-45 +.0-11 —0-271 0-080 
x2 —0-190 +0-033 0-35 +.0-10 —0-292 + 0-067 


Mean value «=: —0-283 +0:060 
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Our best value «=—0-283-+0-060 has been worked out allowing 
for the common errror in y (which is more serious than those in the 
A values). 

We now consider the effect of errors in our assumptions. Even if as 
much as 20% of the ,.-mesons are rapidly quenched in the carbon of the 
Columbia experiment, instead of none as we have assumed, our final 
result for « would be made more negative by less than the above standard 
error. Again, even if it is unplausibly contended that the quenching 
in nuclear emulsion is complete for the gelatin phase instead of zero as 
we have assumed then our final value for « becomes more negative just 
to the limit of its present standard error. The two values for « derived 
from the two different emulsions are then in much poorer agreement, 
however, and this in itself argues for the correctness of our more reasonable 
assumption that the gelatin phase does not quench strongly relative to 
the halide phase. Finally the result is very insensitive to the assumption 
that some p-mesons are quenched immediately in nuclear emulsion and 
that the rest retain their memory until they decay. Even if we make the 
diametrically opposite assumption, which is certainly contrary to the 
evidence of the Columbia experiments, namely that there is no immediate 
quenching and that all their lowered asymmetry in emulsion is due to 
a gradual loss of polarization exponentially with time, our result for « 
merely becomes more negative by 0-02 which is well within the error 
limits. This fortunate insensitivity is due,to a happy choice of the time 
gate in the Columbia experiment. We see that our result for « is 
surprisingly insensitive to the assumptions made in deriving it. It is 
consistent with the extreme value «=-—4 allowed by the two-component 
neutrino theory and relaxation of any of our assumptions moves the 
experimental result into better agreement with that extreme theoretical 
value. Alternatively we may accept the two-component theory} and 
use our result to determine its parameter € (Lee and Yang 1957) which 
measures a relationship between the axial vector and vector couplings, 


E=(faP+lfeP) Mfeta*+haty*). 
We find €=0-85-L0-18. 
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ZIMMERMAN AND PaEvskKy (1955) have reported that the scattering 
cross section of cold neutrons in aluminium calculated on the incoherent 
approximation is about 10% lower than their observed value at 
temperatures in the neighbourhood of 880°K. According to these authors, 
when correction to the incoherent approximation is applied using 
formula (42) of Kothari and Singwi (1955, hereafter referred to as 1), 
the discrepancy increases still further. In view of this we thought it 
necessary to re-examine this problem carefully. 

There can arise three possible sources of error in the calculation of the 
scattering cross section: (i) the use of the Debye frequency spectrum ; 
(ii) assumption of a constant Debye temperature © ; and (iii) an incorrect 
evaluation of the correction to the incoherent approximation. We shall 
now examine each of these separately. 


(i) Squires (1956) has recently calculated the frequency spectrum 
of aluminium from its elastic constants and has shown that so far as 
the total scattering cross section is concerned one gets almost the same 
results whether one uses the exact spectrum or the one given by Debye. 
This, therefore, rules out (i) as a significant source of error. 

(ii) Owen and Williams (1947), from their x-ray diffraction studies in 
aluminium, have concluded that the Debye temperature © decreases 
with increasing temperature. Variation of O with temperature calculated 
from their experimental results is shown in fig. 1. It is seen that © 
changes by about 12° in going from 300° to 800°K. Taking account 
of this variation and using .eqn. (40) of I, which has been derived on 
the incoherent approximation, we have recalculated the scattering cross 
section as a function of temperature for 8 A neutrons. In fig. 2, curve A 
represents the scattering cross section derived by assuming a constant 
value of @, whereas curve B represents the same quantity calculated by 
taking account of proper variation of @. At 800°K, a 12° change in © 
produces nearly a 12% increase in the cross section. 

(iii) It appears that formula (42) of I does not give the correction to 
the incoherent approximation correctly, and therefore, we have not used 
this in the present calculations. We evaluate the one-phonon cross 
section both by summing over the reciprocal lattice vectors (Weinstock 
1944) and also by replacing the sum by an integral (eqn. (30) of I with 
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Debye temperature as a function of temperature of aluminium. Circles 
denote the values computed from the experimental data of Owen and 
Williams (1947). 
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Scattering cross section of 8A neutrons in aluminium as a function of 
temperature. Curve A: multi-phonon cross section calculated on the 
incoherent approximation using a constant value of 9 (@=394°K). 
Curve B: same as curve A but using values of © as given by fig. I. 
Curve C: obtained by adding to curve B correction to the incoherent 
approximation. 
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J=1). The difference of these two values gives us the correction to the 
one-phonon incoherent cross section. At 800°K this correction is about 
+4, and will decrease rapidly with the order of phonon process. And 
since the contribution of the multi-phonon processes to the total cross 
section is nearly 30°, at this temperature, the correction to the incoherent 
multi-phonon cross section will not alter the final results by more than 
one per cent. We may, therefore, justifiably assume that the correction 
to the one-phonon incoherent cross section also represents the correction 
to the total incoherent cross section. Further, the correction has been 
evaluated for a fixed value of @. The effect of the variation of @ on 
the correction being negligible, has not been considered. 

It may be mentioned that in evaluating the sum over the reciprocal 
lattice vectors we have used formula (51) of Weinstock (1944) but unlike 
him we have retained the Debye-Waller factor within the sign of 
integration over his variable A, and have carried out the integration 
numerically. 

If we now add to the cross section represented by curve B of fig. 2, 
the correction discussed above, we get curve C which gives the final 
values of the calculated cross section. 

Apart from the statement of Zimmerman and Palevsky (1955) the 
only published value of the total cross section o,—o,-+o,, for aluminium 
in the temperature range considered is at 800°K (Hughes and Harvey 
1955) and is 1°83 barn. If we subtract from this the value of the 
absorption cross section, o,, for 8A neutrons, o,=1:02 barn, the 
scattering cross section is o,=0-81 barn. The calculated cross section 
(curve C of fig. 2) at this temperature is 0-79 barn. Thus the present 
theory and experiment agree to within 3% at this temperature. It may 
be mentioned that even if we take the correction to the incoherent 
approximation properly but neglect the variation of 6 with temperature, 
the calculated value of the cross section at 800°K is 15° lower than the 
observed value. In other words at high temperatures the effect of the 
variation of © with temperature on the cross section is more important 
than the correction to the incoherent approximation. 

Although the theory has been found to agree well with experiment 
at 800°K, it would be desirable to have more experimental data available 
on the scattering of cold neutrons in the high temperature range. 
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A Search for Natural Radioactivity in Vanadium 


By R. N. Grover} and D. E. Warr t 
Department of Natural Philosophy, The University, Glasgow, W.2 


[Received January 23, 1957] 


Mass spectrometer measurements by Johnson (1952) have indicated 
that the odd—odd nucleus °°V has an energy of 2:39-+-0-12 mev available 
for decay to *°Ti and 1-19-+0-12 mev for decay to *°Cr. Its spin has 
been measured to be 6 (Kikuchi et al. 1952, 1953, Baker and Bleaney 
1952). -decay should occur to the ground state of ®°Cr with spin 
change AJ=—6, since its first excited state is at 9-5 Mev (Endt and 
Kluyver 1954). Only one excited state of °°Ti exists below 2-5 mev, 
namely a 2+ level at 1-58 mev (Pieper 1952). Hence an electron capture 
transition with 4/—4 accompanied by a y-ray of energy 1-58 Mev might 
be expected. 

Unsuccessful attempts to detect activity have been made by Sheline, 
(1954), Strome (1954) and Selig (1954). The longest minimum half-lives 
have been quoted by Heintze (1955). However, the ferrovanadium 
source used was stated to be contaminated. Also in deducing the 
minimum half-life for K-capture no allowance appears to have been 
made for the facts that firstly, only a small fraction of the 4-5 kev x-rays 
could escape from a thick source (27°; from a 10 mg/cm? source) and 
secondly, in only 25% of the K-capture transitions are x-rays emitted, 
Auger electrons, which would be absorbed in the source, being emitted 
in the remaining 75°. Hence the minimum K-capture half-life which 
can be deduced from Heintze’s experiment is considerably shorter than 
the quoted value of 3x 101" yr. 

A search for B-rays was carried out using the proportional counter, 
the apparatus and its underground site having been previously described 
(Glover and Watt 1957). The source of vanadium pentoxide was of 
average thickness 10 mg/em?. For the background a similar source of 
chromium oxide was used. The difference between source-+ background 
and background counting rates over a period of 8 hours was +-113-+- 137, 
Assuming a detectable limit of thrice the statistics and an isotopic 
abundance of 0-25°, (White et al. 1956), then, taking into account f-ray 
absorption in the thick source, the half-life for emission of B-rays of 
energy >50 kev must be >2:-4x10' yr. In case the f-rays were 
emitted predominantly with very low energy as in the §’*Rb spectrum, 
the measurements were extended so that pulses of energy > 1-5 kev 
could be detected, but no activity was observed. 
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Using the 10 mg/cm? source the K x-ray region was closely examined, 
the difference between the source-+ background and background counting 
rates over Shr being +45--47. Assuming the same detectable limit, 
then, allowing for absorption in the source of the x-rays which would 
be produced in only 25° of the transitions, the K-capture half-life must 
be >4x 108 yr. 


‘| 


133 MEV 


COUNTS PER CHANNEL 


The y-ray spectrum of °°V. Each cross denotes the average excess counting 
rate over a group of ten channels centred about its position. 


A search for y-rays was undertaken with the scintillation spectrometer 
shielded by 4 in. lead lined with 4 in. steel. The 2in. Nal crystal was 
surrounded by 110g of vanadium pentoxide and by the same mass 
of chromium oxide for the background run. The source-+ background 
and background counting rates above an energy of 200 kev were con- 
sistent over several periods adding up to 10 hr the totals being 28503-+ 169 
and 27644-++166 respectively, an excess of 859+237. Since this excess 
counting rate was only 3% of the total background, the statistics were 
very poor as can be seen from the figure. There is no appreciable peak 
near 1-58 Mev but it seems significant that the spectrum drops sharply 
near 1-33 Mev, the energy at which the Compton edge would be expected. 
Also the result obtained from averaging the excess counting rate in each 
group of ten channels (see figure) suggests that the counts below 1-33 Mev 
are spread evenly over the spectrum, as would be expected for a Compton 
distribution. Failure to detect B-activity eliminates the possibility of 
B-contamination. The «-activity of both the vanadium pentoxide and 
chromium oxide was found to be only 2 c.p.m., the minimum level of 
impurity we have found to be attainable. Although tests have been 
carried out on the oxides of six different elements which gave this «-rate 
of 2 c.p.m., no y-activity has been found to be associated with the 
z-activity. Therefore, it is unlikely that the activity is due to any 
impurity. By surrounding the crystal with a similar mass of potassium 
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bicarbonate in the same geometry, the y-activity was estimated to 
correspond to an electron capture half-life for °°V of (4:0-L1-1)x 10" yr. 

Comparison of this half-life with the lower limit attainable in the 
proportional counter search for K x-rays shows the impossibility of 
detecting these x-rays using a thick source of unenriched material. It is 
estimated that if the search for x-rays was carried out using a minimum 
of 250mg of vanadium pentoxide, 50 times enriched in ®°V, then a 
significant counting rate should be obtained. 
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A Search for Electron Capture in '7°Lu 


By R. N. Guovert and D. E. Warrt 
Department of Natural Philosophy, The University, Glasgow, W.2 


[Received January 23, 1957] 


RECENT investigations by Arnold and Sugihara (1953), Arnold (1954) 
and Dixon et al. (1954) have conclusively confirmed the 6~-decay scheme 
of naturally radioactive 17Lu first suggested by Goldhaber and Hill 
(1952) on the basis of unpublished work by Scharff-Goldhaber (1952). 
Since 176Lu is an odd—odd nucleus and the middle member of the 
isobaric triplet (17®Yb,», 17%Lu,,, 17°Hf,,), decay by electron capture 
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might also be expected. However, it was not observed by Scharff- 

-Goldhaber (1952) nor by Arnold (1954), the latter concluding that its 
intensity was certainly <10% of that of the B-rays. From an examination 
of the L x-ray peak, Dixon et al. (1954) tentatively proposed that 32-1 hs 
of the disintegrations took place by electron capture. The ground state 
spin of 176Lu being 10, it is unlikely that electron capture occurs to the 
ground state of !7*Yb, since its spin is zero. Hence a search was under- 
taken for y-rays or conversion electrons which could be attributed to an 
electron capture branch. 
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_ The y-ray spectrum of '’$Lu. 


Both the scintillation and proportional counters were located 90 ft 
underground as described previously (Glover and Watt 1957). The 
y-ray spectrum from 500 mg of lutecium oxide placed on top of the 
2in. Nal crystal is shown in the figure. The source counting rate was 
700 c.p.m. compared with the background of 80 c.p.m. The energies of 
the y-rays were accurately determined to be 206-L5 and 311-7 kev in 
good agreement with the values of 203 and 306 kev found by Arnold 
(1954). The peak at 523+-10 kev was due to simultaneous detection of 
these two y-rays by the crystal. Similarly that at 366 kev corresponded 
to addition of the 311 kev y-ray and a 55 kev x-ray. The addition peak 
for the 206 kev y-ray and an x-ray would occur in the valley between 
the two y-ray peaks and was obscured in the present spectrum. The 
incompletely resolved peak close to the K x-ray peak was due to 89 kev 
y-rays which were not internally converted. There was, therefore, no 
indication of any peak not attributable to the known y-radiation, 
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120 mg of Lu,O; was spread over an area of 800 cm? of the proportional 
counter lining giving a uniform source of thickness 0-15 mg/em?. The 
anticoincidence Geiger ring was not required. The pulses were fed along 
a 400 ft cable to the Hutchinson—-Scarrott kicksorter. There was no 
indication of any electron peaks other than those known to be due to 
internal conversion of the 89 kev y-ray. 

The half-life was estimated from the observed y-ray spectrum. Allow- 
ance was made for the Compton distribution due to the 206 kev y-ray 
(see figure). Then the remainder of the counts in the spectrum represent 
the total number of 311 kev y-rays detected by the crystal either alone 
or in coincidence with a 206 kev y-ray or an x-ray. A source of potassium 
chloride was placed in the same geometry and its spectrum taken. 
Knowing the specific y-activity of potassium and calculating the 
average thickness of the crystal, the geometry was deduced using 
a value of y,,, the mass absorption coefficient in NaI for 1-46 Mev y-rays, 
interpolated from the tables given in the Appendix to Siegbahn (1955). 
Hence, using the value of y.,, for 311 kev y-rays obtained from the same 
tables, the half-life was calculated to be (2:1+0-2)10!% yr. This is in 
excellent agreement with the value of (2-15+0-10) x 10!° yr. found by 
Arnold (1954) also from examination of the y-ray spectrum, but not 
with that of (4:6-0-3) x 101° yr. obtained by Dixon et al. (1954) from 
examination of the f-spectrum. The half-life was also calculated from 
our observed f-spectrum, but an «-rate of 40 c.p.m. was obtained from 
the source. In the course of investigations on the natural radioelements 
with Z<80 it has been found necessary to check the purity of all sources, 
measurement of their «-activity, which is due to slight contamination 
by natural radioelements of the uranium and thorium series, being the 
most sensitive method. High «-rates have often been found to be 
associated with contamination, either 8 or y, and hence the half-life 
calculated from the 17*Lu B-spectrum (2:8 x 101° yr.) cannot be considered 
significant. 

It is concluded that there is no evidence of y-rays or conversion electrons 
which could be associated with an electron capture branch and that the 
branching ratio of (83-+-1)% proposed by Dixon et al. (1954) might rather 
be regarded as an upper limit. There appears to be a marked disagreement 
between the half-lives reported from y- and f-ray measurements. Un- 
fortunately for the reason given above, it was not possible in the present 
investigation to determine the f-ray half-life. 
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On the Surface Free Energy and Specific Heat of a Metal 


By R. STRATTON 
Metropolitan Vickers Electrical Co. Ltd., Manchester, 17 


[Received January 3, 1957] 


Tue surface contributions to the lattice and free electron specific heats of 
metals were derived by Montroll (1950) (to be referred to as M) in a paper 
which has only recently come to the author’s attention. The lattice 
contribution has however been calculated on a different basis, Stratton 
(1953 a) (to be referred to as S), and the electron contribution can be 
similarly derived by extending to non-zero temperatures the calculation, 
Stratton (1953 b) (to be referred to as I) on the surface free energy of a 
metal in the normal state at absolute zero of temperature. 

Since our method differs in some important details from that employed 
by Montroll the calculations are briefly presented although the results 
differ from those given by Montroll only by numerical factors which are 
not large enough to produce any changes of magnitude. The opportunity 
will also be taken to correct an arithmetical error in the derivation of the 
ratio (w) of the surface free energy difference (Af) between the super- 
conducting and normal states, to the corresponding bulk term, Stratton 
(1953 c) (to be referred to as II). 

(a) Comparing the expressions for the total number of elastic modes 
with angular frequency less than w given in M (eqn. (3)) and S (eqn. (68)) 
it will be seen that the term 10,~1, or #(C,-?+ 20-2), in the former is 
replaced by ($+ Y)?—J,;)C,-? in the latter (C, and CO, are the velocities of 
longitudinal and transverse elastic waves respectively, Y, and J, are 
functions of Poisson’s ratio o defined in §). 

The difference arises because periodic boundary conditions are assumed 
in M so that the elastic waves can be separated into longitudinal and 
transverse modes whereas it is shown in 8 that for a free surface the 
normal modes have both transverse and longitudinal components. 
(Similar conclusions would hold for a fixed surface.) The results in M 
can thus be used if O, is replaced by 


102 (5—60)0, 
HFY?—I, 8(1—o)G+¥e—h) 
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For the three values of o considered in 8, c=0, }, 4; L(c) is given by 
0-51, 0-69, 2-04, respectively. In particular the ratio between the bulk 
7* and surface 7? term in the specific heat will be multiplied by the 
factor L(c). 

(6) The density of electron states whose wave number & lies in the 
range k, k-+-dk is given by (I, eqn. (20a)) 


plk)dk=(V/n®)[1+-d/L,|[1—m/(4kb,)Wedk . 2. . (2) 


where V is the volume of a thin flat plate of thickness 2, and d is the 
distance separating the position of the lattice and potential barrier 
surfaces. (It has also been assumed that the potential barrier height is 
infinite since, as the calculations in I show, a finite barrier causes a 
correction to the surface free energy of only a few per cent.) The second 
term in square brackets represents a correction due to states with the 
z-component of the wave number vector equal to zero, which are included 
in integrals over wave number space with weight } but should be excluded. 
The corresponding expression for the density of states given, without 
derivation, in M differs from eqn. (2) in that d is zero and there is a 
positive sign inside the second square bracket. 
At a temperature 7’ the total number of electrons is given by 


vy % p(k)dk 

5 . / 0 1-+-exp [((H—C)/RT] 
where ¢ is the Fermi energy, # is the electron energy and k is Boltzmann's 
constant. Substituting for p(k) from eqn. (2) and expanding the integral 
in the usual way, neglecting terms of se ae gives 


rl SCEY hie eC) ]] 


where k,, is the maximum wave number hun eee to ¢ and the 
suffix 0 denotes values at 7=0. Contrary to the assumption that d=0 
in M it has been emphasized in I that d must be chosen so as to make £,, 
(or ¢) surface independent, i.e. the term involving V/L, in eqn. (4) must 
vanish. Hence 


(3) 


one L477 8) (RL Cg)2) a see 3 | oe AD) 
and 
eee 2 Fe) og) Wee ee! ve ns is o> 0 (9) 
Since £,,, in eqn. (5) still involves 7’ we have by using eqn. (6) 
==(37/Sk,.6)[1—(rt/12)(RT/E)7).  . 7) 


Jn a similar manner the internal energy 


el is Ep(k)dk 


=| ———__, (8) 
Jo 1texp (H—£)/kT 


turns out to be 


u= ent] {1+ peda =~). eta, = sae 9) HH 


where ¢ and d have been replaced from eqns. (6) and (7) respectively. 
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Thus the specific heat is 


1dU i Tw! ; 

Cy= Nav k— 9 ee ~|'+a 40 |: . . . . (10) 
The value of C', found in M, expressed in the same form as eqn. (10), has 
the factor 1/8=0-125 instead of 7/40=0-175. (N.B. The result in M 
can be found by putting 65=?k4,92/2m= (#?/2m)(372N/V)9 and L-1=S/V 
where S is equal to the surface area). 

(c) The last term in the equation immediately following eqn. (54) im IT 
should be omitted. This requires the following replacements of numerical 
factors : 47 by 4x in eqn. (55), ta by 47 in eqn. (56), 7 by $7 in eqn. (57) 
and 1-9 by 0-24 in eqns. (65) and (67). Thus the values of 4/f and » for 
mercury become 0-07 erg cm? and 1-2 x 10-7 cm respectively. 
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A New Method for Decorating Dislocations in Crystals of 
Alkali Halidest 


By D. J. Barper, K. B. Harvey and J. W. MircHeno 
H. H. Wills Physical Laboratory, University of Bristol 


[Received January 4, 1957] 


§ 1. InTrRopucTION 


LinEAR imperfections within transparent crystals which could, with 
little doubt, be recognized as dislocations were first observed in crystals: 
of silver chloride and silver bromide by Hedges and Mitchell (1953). 
Amelinckx et al. (1955) and Amelinckx (1956) then showed that disloca- 
tions could be made visible in crystals of sodium chloride by heating the 
crystals in contact with sodium—a method for rendering imperfections 
visible which had previously been employed by Rexer (1931, 1932 a, b). 
Boesman et al. (1956) next found that a trace of silver Pade had to be 
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present in the alkali halide before lislocations could be reliably decorated 
by the Rexer method. This method does, however, suffer from the 
disadvantage that the crystals have to be heated, in contact with sodium, 
to about 100°c below the melting point so that annealing structures 
alone can be observed. Decoration of dislocations at considerably lower 
temperatures has since been achieved by Van der Vorst and Dekeyser 
(1956) who have heated mixed crystals of sodium chloride with a few 
per cent of silver chloride in hydrogen at 500°c for periods of 24 hours or 
more. Silver then separates along the dislocation lines near the surfaces 
of the specimens. 

The purpose of this note is to describe another method for decorating 
dislocations in crystals of alkali halides at relatively low temperatures. 


§ 2. EXPERIMENTAL MeETHOD 


The method consists in heating crystals of alkali chlorides and bromides 
with the corresponding gold halide (HAuCl,, HAuBr,) in evacuated, 
sealed-off tubes. Specimens are cleaved from large single crystals 
and deformed by compression or bending. They are placed in a hard 
glass tube together with a small quantity of the gold halide and the tube 
is evacuated and sealed off. The gold halide is next melted at the lowest 
possible temperature and the tube is slowly rotated until the specimens 
are covered with a very thin layer of it. The tube is then heated to a 
constant temperature in the range between 500°c and 750°c for a period 
of up to four hours and allowed to cool to room temperature. The 
specimens are not coloured when the tube containing them is removed 
from the furnace at the high temperature. A deep blue-green or blue- 
violet coloration develops during cooling to room temperature. After 
being taken from the tube, they are cleaved into thin plates which are 
mounted in Canada Balsam on 3 in. x 1 in. slides for microscopic examina- 
tion. It has been found that the dislocations are decorated with gold by 
this method to a depth below the surface which depends upon the 
particular alkali halide and the time of heating. The definition of the 
dislocations depends upon the rate of cooling. Under favourable con- 
ditions, the decoration is quite continuous along the dislocation lines and 
there is very little separation of particles of colloidal gold elsewhere. The 
mechanism of decoration, which is of considerable interest, will be dis- 
cussed in a more detailed paper. It appears that the configuration of 
dislocations which is present in the crystal at the end of the period of 
heating in the furnace is made visible, with little modification, by the 
precipitation of gold along the dislocation lines during cooling. 


§ 3. EXPERIMENTAL OBSERVATIONS 


Our principal aim in this work was to develop a method for decorating 
dislocations on slip planes in slightly strained cleavage bars of potassium 
bromide. This aim has been accomplished. It is hoped that the new 
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observations will facilitate the interpretation of the far more complex 
behaviour of dislocations on slip planes in crystals of silver bromide (see 
Mitchell 1957). 

Figure 1 (Plate) is a photomicrograph (1250) showing individual 
dislocations on slip planes in a crystal of potassium bromide. The 
crystal was very slightly deformed by bending and then heated with gold 
bromide in a sealed tube at 550°c for four hours. It was then cleaved into 
plates, 0-4 mm in thickness, with surfaces normal to the axis of bending. 
The dislocation lines were found to be decorated throughout the crystal 
by the precipitation of gold uniformly along them. The points in the 
photomicrograph are defined by the passage of dislocations through the 
focal plane of the microscope. They represent approximately equally 
spaced edge dislocations which are parallel to each other on any 
particular glide plane. Phenomena due to the climb of dislocations are 
not observed at this temperature. Some relaxation must occur, however, 
because the strain birefringence associated with the deformation dis- 
appeared when a specimen was annealed at a much lower temperature. 
Almost identical observations have been made with similarly deformed 
crystals of sodium chloride heated with gold chloride for four hours at 
a temperature of 650°C. 

The climb of edge dislocations in crystals Br sodium chloride begins 
about 600°c and many interesting details of the behaviour of individual 
dislocations and of arrays of dislocations during climb have been observed. 
These properties cannot be adequately illustrated with photomicrographs 
and three dimensional wire models which can be photographed are being 
made. 

Figure 2 (Plate) (x 1250) shows two walls of dislocations in a crystal of 
sodium chloride which was heavily deformed by bending and then 
polygonized by heating at 700°c for four hours. The lines at the top 
and bottom of the photomicrograph are the traces of glide planes upon 
which individual dislocations cannot be resolved. Similar configurations 
have been observed by Amelinckx in crystals of sodium chloride, con- 
taining a trace of silver chloride, which had been deformed and heated in 
contact with sodium. 

Specimens of sodium chloride, lightly deformed by bending and then 
heated with gold chloride for four hours at 650°c have been found to 
contain many groups of closed dislocation loops with a common axis 
coinciding with one of the twelve 110) directions of the crystal. 
A typical system is illustrated in fig. 3 (Plate) (x 4650) in which (a), (6) 
and (c) were taken with the focal plane above, coinciding with, and below 
the axis of the system. The trace of a slip plane appears in the left hand 
part of the field. We assume, as a working hypothesis, that these are 
derived from segments of edge dislocations by climb as discussed by 
Bardeen and Herring (1952). We believe that the source is associated 
with a cavity in the crystal which, in fig. 3 (b) (Plate), appears as a dark 
spot near the centre of the system with a loop expanding from it on the 
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left-hand side. Many thousands of these systems have been examined, 
some with one or two closed loops, and many with twenty or more 
loops. 

A close relationship has been observed between these systems of closed 
loops and helical dislocations which appear, at least superficially, to be 
of the same type as those first observed in crystals of calcium fluoride 
heated in sodium and described by Dekeyser and his co-workers. A 
typical helical dislocation in a crystal of sodium chloride, which had beer 
heated to 650°c for four hours, is illustrated in figs. 4 (a), (b) and (c) 
(Plate) (3150). This helical dislocation appears to have been derived 
from a system of closed loops associated with a central source (shown as 
a dark spot in fig. 4 (6) (Plate) as a result of interaction with a dislocation 
which extended from the slip band, the trace of which is shown in the 
left-hand part of the field, to an adjacent slip band, which was located at 
a distance of approximately 8 times the width of the photomicrograph 
totheright. The resulting dislocation was followed continuously through- 
out its length. It was approximately straight from the slip band to the 
point where it entered the field on the right. It then followed a helix, 
the turns of which first contracted in radius to the position of the assumed 
assumed source, and then expanded again. The large individual loops 
of the helix, which are clearly shown in fig. 4 (6) (Plate), are consistent 
with inherent instability. 

The decoration of dislocations with gold is such that they may be 
studied with high power microscope objectives in the bright field illumina- 
tion provided by a large aperture condenser. This has not previously 
been possible with crystals of alkali halides. Figures 3 and 4 (Plate), 
which were taken with Kodak Panchromatic Microfile film using bright 
field red illumination, give some indication of the detail which can be 
observed. 


$4. CONCLUSION 
This method will supplement those already employed by Dekeyser and 
his co-workers. It allows observations of the dislocations associated with 
deformation structures to be made and permits details of the behaviour 
of dislocations during climb to be deduced. The fact that decoration is 
continuous along the dislocation lines is a.particularly valuable feature. 


Note added in proof—tIt has been found since the submission of this paper 
that dislocations in crystals of rock salt may be very conveniently decorated 
with gold simply by heating them to between 600 and 709°o for four to six 
hours in a sealed tube containing a spiral of thin gold wire and chlorine at a 
pressure of 300 mm Hg (measured at room temperature). The gold precipitates 
along the dislocation lines at about 400°c during the cooling of the specimen 
and makes the configuration of dislocations which was present at that 
temperature visible. 
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Diffusion of Zinc in Zinc Oxide 


By J. P. Roperts and C. WHEELER 
Royal Aircraft Establishment, Farnborough, Hants. 


[Received February 26, 1957] 


In the absence of large single crystals, tracer measurements of diffusion 
in zinc oxide have been confined to polycrystalline material (Miller1941, 
Lindner et al. 1952) and to vapour exchange experiments on small crystals 
(Secco and Moore 1955, Miinnich 1955, Spicar 1956). The results show 
considerable dispersion. 

We have concluded from work referred to below that diffusion of zine in 
dense polycrystalline zine oxide is predominantly along grain-boundaries 
up to at least 1300°c. This supports the previous inference of Lindner 
et al. (1952) and suggests a reason why, apart from any question of non- 
stoicheiometry, agreement between data from polycrystalline material 
evaluated only in terms of uniform diffusion and data from single crystals 
is not to be expected. 

We determined by a sectioning method the profile of diffusion of zine-65 
in sintered material (density «93° %, grain-size 10-10? 1, purity usually 
>99-9%) at 800-1300°c in oxygen and in argon. The tracer was initially 
contined to a surface layer (~1,). After heating, the diffusion zone was 
progressively ground away. The residual gamma activity (A,) of the 
specimen was measured at various depths (a) (Roberts et al. 1956). The 
figure shows typical results. 

Constant ¢ log 4,/dx, as observed here, is to be expected if grain- 
boundary diffusion is predominant because the penetration of tracer 
then approximately obeys a simple exponential law, as predicted (Fisher 
1951, essentially supported by Whipple 1954, Borisov and Lyubov 1955) 
and as observed in several investigations of metals. Confirmation of 
predominant grain-boundary diffusion in the present material was 
obtained by autoradiography. 
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_ From @ log a/dx (where a is the activity per unit thickness at x), which 
here equals 0 log A,/dz, Fisher’s theory allows evaluation of Dead’, 
where D, 6 and PD’ are respectively average values of the volume diffusion 
coefficient, grain-boundary width and grain-boundary diffusion coefficient. 
(The shape of our curves near x=0 was governed by surface roughness, 
which the friability of the material prevented us from reducing consistently 
below about +10 uw.) For diffusion in 1 atmos. oxygen, observations 
up to the present give 
D4/8D’=104*? exp e esa reeleisled 


We are investigating the reason for the appreciable scatter, also found 
with polycrystalline material by Lindner ef al. (1952), and the result will 
be published later. The range of D1?/5D’ values obtained is however no 
greater than one order of magnitude at all temperatures used (the 
reversed signs of the uncertainties in the above equation will be noted) 
and it is possible to make certain deductions now. 
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Diffusion of zine-65 in polyerystalline zine oxide. 


No measurements on single crystals in oxygen have been reported, but 
in | atmos. oxygen it is likely that 
—89 kcal mole -)em? baie 
RT 


‘approximately. (The most probable value of activation energy given 
‘by Spicar (1956) from crystals in zinc vapour is combined with 
D~10-13 em? sec-! at 1300°c, deduced from our autoradiography.) 
. With the usual assumption that 5=5 x 10 § cm at all temperatures, 
—69+12 keal me) Bnd coe 
rye 


D0 Vexp ( 


D/=108*3 exp { 
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giving D’ within one order of magnitude between 800 and 1300°c. The 
ratio of grain-boundary/volume activation energy is 0-77--0-13, which is 
reasonable. 

D12/8D’ in 1 atmos. argon (<10-? mm Hg oxygen) (see figure) was 
smaller than in | atmos. oxygen by up to an order of magnitude between 
800 and 1300°c. The fractional excess zine content of some of our 
sintered material at 1300°c, measured by a chemical method developed 
in this laboratory (Allsopp 1957), changed from 3 x 10~¢ in | atmos. oxygen 
to 210-5 in 1 atmos. argon, suggesting a connection between D1?/5D’ 
and composition. Since composition might be expected to affect D more 
than D’, the present results imply that wnder our conditions D may 
decrease with increasing excess zinc. This picture is less easily reconciled 
with the commonly assumed model of zinc oxide containing interstitial 
zine than with a model involving only anion and cation vacancies. 
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Circular Polarization of y-rays : 
Further Proof for Parity Failure in 8 Decay 


By H. ScHorrer 
Cavendish Laboratory, Cambridget 


[Received March 14, 1957] 


Leb and YANG (1956) suggested several experiments for testing the 
conservation of parity in weak interactions. Two of these have been 
performed (Wu et al. 1957, Garwin et al. 1957{) and have shown that 
parity is not conserved. Results of a third experiment (thought im- 
practicable by Lee and Yang) are reported here. They confirm the 
expectation that the y-rays emitted after B-decay at an angle 6 relative to 
the 6-particle should show circular polarization proportional to cos 0. 


SE EE 
+ Permanent address; Physikalisches Institut, Erlangen, Germany. This 
work was performed during a visit to the Cavendish Laboratory. 
} The author would like to thank Dr. D. H. Wilkinson for providing pre- 
prints of these papers and for other valuable information. 
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Circular polarization of y-rays can be detected through the Compton 
scattering by polarized electrons in magnetized iron. Of the three 
possible experimental arrangements | transmission (Gunst and Page 1953), 
back scattering (Clay and Hereford 1952) and forward scattering 
(Wheatley et a/. 1955)|, forward scattering seemed most favourable here. 

The scattering cylinder (see figure) was thick enough (4 in.) to make 
the scattering by the magnetizing coil negligible. The mean scattering 
angle was about 55° which is the optimum for a quantum energy of 
1-28 Mev. The effect is not smeared out appreciably by the geometry 
used. No influence of the stray field on the electrons was observed. 
B-y coincidences were counted on a standard fast-slow circuit (resolving 
time 0-lj sec). Pulse height selectors in the slow circuits allowed dis- 
crimination against back-scattered y-rays, annihilation radiation and 
slow £-particles. Random coincidences could be counted by inserting a 
delay of about 0-5 sec in the A-line. 


# and P expressed as parts per hundred 


E 2 . E 

— for e ee 1 for 
§—180° Exper. | theor. 9~90° 
© 4b 16. b964 30-41 .'| 40-67 
Co | 10.36 | + 4] 20-07 | ~2/2| 10.49 
ee Ee ae eee SeRO 
#Na | 9.52 | 16 | £008 | 71/2] .o-76 
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The magnetic field was reversed every 5 minutes in order to eliminate 
drifts in the coincidence resolution. Any variations in detection sensiti- 
vity—in particular those caused by the reversal of the magnetic stray 
field (change of single counting rates <1°,,)—were eliminated by forming 
R, the number of coincidences divided by the product of the single counts. 
Let us define the polarization effect H=2(R,— R,)/(R,+R,) where k, and 
R, refer to the magnetization pointing toward and away from the source, 
respectively. The observed effects at 180° are shown in the first column 
of the table. Co and ?2Na were chosen because they have allowed 
£-transitions with 4./=1, which should give maximum effect. Co emits 
two y-rays in a 4-2-0 cascade; in that special case the polarization of 
the second y-ray is the same as that of the first (Tolhoek and Cox 1953). 
Hence the results can be discussed as if (o®—like ?2Na—emitted only 
one y-ray. 

The polarization of the y-rays is shown in the second column, computed 
from H=2fP\c,,/o,|, where f is the fraction of oriented electrons in the 
iron (0-08 at the measured flux density of 20 000 gauss) and o, and oy, are 
the polarization-dependent and independent parts, respectively, of the 
differential Compton cross section (c/a) =9-52 for the geometry used). 
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As a check, some measurements were done at 90°. A small effect (of 
the order of the statistical error) was found (column 5); actually the 
average angle was somewhat larger than 90°. On counting random 
coincidences, no effect greater than the statistical error (~1%) was 
observed. 


Sel f-counter (plastic scintillator), Se2 y-counter (Nal(Tl)—crystal), L1 and 
L2 Perspex light guides, S source, M cylindrical magnet with magnetizing 
coil C, Pb lead absorber. 


Calculations by Skyrme (private communicationy), based on the two- 
component theory give P=«(v/c) cos 6, where « depends on the type of 
the 6- and y-transition (and is zero for pure Fermi interaction). Theo- 
retical values of « are given in column 4 of the table. 

The experimentally determined sign for ®°Co fits the experiment of 
Wu et al. The different signs for Co and Na prove that the anti- 
neutrino emitted together with a negatron has the opposite screw sense 
from the neutrino which accompanies the positron decay. The experi- 
mental values of « are somewhat higher than the theoretical ones but 
the deviations are within the experimental errors. 
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Note added in proof :—In the meantime, 24Na (allowed transition, 4./ == 0) 
has been investigated; only a small asymmetry H=--0- 56% -+0-38%, 
was found. The thoery (K. Alder, B. Stech and A. Winther, pirate 
communication) predicts a large asymmetry except for practically 
iis Fermi or Gamoy—Teller transitions; our result means that 


2 2 2 2 2 
C7?|M ex P(C_?|M ar P+C?|M >|?) 
is either >99%, or <2%. In cases where both interactions are present 


in comparable amounts it would be possible to find the relative sign of 
C,and Cy. I wish to thank D. J. Chandlin for helpful discussions. 
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The Weak Lines in the Radium D Gamma Spectrum 


By A. C. L. Barnarpt and B. C. Haywoopt 
University of Birmingham 


[Received January 15, 1957] 


Prior to 1939 the decay of radium D (?!°Pb) was regarded as a simple 
and well understood f-decay to the 46-5 kev level of *1°Bi (the By, 
transition). However Lee and Libby (1939) showed that only about 
three-quarters of the disintegrations follow this mode, and this prompted 
re-examination of the RaD y-spectrum. Results reported between 1939 
and 1954§ showed a number of weak lines, summarized by Wu ef al. 
(1953) as follows. 


Energy (kev) 46-5 42-6 37+0-5 | 3140-8 | 23-2+0-6 | 16-10-4 | 7-30-7 
Intensity % | 2-8-40-6 | 0-2+0-1 | 0-2+0-1 0-4-4 0-2 ~l a ~10 
Wu succeeded in showing that the lines at 7-3, 16-1 and 23-2 kev 


originated from secondary, atomic processes, but no conclusion was 
offered on the 31, 37 and 42-6 kev lines. Ignoring these, a decay scheme 


es at A.E.I. Ltd. Research Laboratory, Aldermaston Court. 
* Now at A.E.R.E., Harwell. 
S § The summarized results of 16 papers in this period are tabulated by 


Fink et al. (1956). 
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was postulated in which (92-45)% disintegrations go through the 
46-5 kev level, a ground-to-ground f-transition (8) transition, maximum 
8 energy ~65 kev) being suggested as accounting for the remaining 
(8--5)% disintegrations. Fink et al. (1956) are repeating some measure- 
ments with the object of showing that B,, accounts for all disintegration, 
but Stanners and Ross (1956) claim to have detected the 8) transition in 
(15-5+-3-5)%, of the disintegrations. 
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Dependence of gamma spectrum of radium D on chemical form of source. 
(Curves uncorrected for counter efficiency.) 


-The weak y lines have been the subject of conflicting reports in which 
both energies and intensities have varied considerably. For example, 
Ewan and Ross (1952) put an upper limit of 0-1°% on the intensity of 
the 31 kev line reported earlier, whereas Damon and Edwards (1954) 
report a definite line at 31 kev. Fink ef al. (1956) suggest that the radiation 
in the 31 kev region originates from the selective action of an aluminium 
absorber on the B-bremsstrahlung from Rak (if present) but the effect of - 
this is to introduce only an extremely broad, low hump starting at 20 kev 
and merging into the 46-5 kev peak. The reports suggest that a more 
definite peak has been seen by some workers. In. another paper 
Fink and Robinson (1956) suggest that the 42 kev line could be caused by 
Compton scattering of the 46-5 kev quanta. 

The comparatively intense line which had been reported at 7 kev 
was shown both by Wu ef al. and by West (quoted in Jaffe and Cohen 
1953) to arise from fluorescence induced in nearly copper, iron or nickel 
(e.g. brass counter walls) by radiation from the source. This effect is 
now well known, but the similar effect, fluorescence induced in atoms 
present in the source (either as impurity or as a chemical constituent) 
has not been mentioned in the literature. 
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Electrons from the conversion of 46:5 key quanta are emitted in 
(85--5)% of the disintegrations, about 75°, of these conversions being 
in the L-shell, and unconverted quanta of this energy in about 7°% of 
the disintegrations (Wu et al. 1953). This is much smaller than the 
conversion coefficient derived theoretically by Gell-Man et al. (1952) but 
corrections to the theory by Sliv (1951), Sliv and Listengarten (1952), 
Rose (1955), Ward and Mihelich (1956) and Church (1956) give better 
agreement. The f, transition will produce low energy electrons in some 
cases, but it seems that 8 or y radiation of energy greater than 30 kev 
is emitted in about 90°, disintegration. Such radiation could induce 
K-fluorescence in atoms up to Sm (atomic number, Z=—62) but the 
intensity would be greatest below Z=54, the highest element that can 
be excited by the L-conversion electrons. 

This effect in source atoms was noticed in this laboratory in 1954, 
and is illustrated in the figure. Bias curve | was obtained with an 
argon-filled, thick window (1:5 mm glass) proportional counter, and a 
solid external RaD source in the form of lead chloride, mounted on a 
250mg ecm copper tray. Bias curve 2 is for a lead iodide source, under 
identical conditions, and shows the iodine K—x-ray very prominently. 
The intensity of this line is about 20° of that of the 46-5 kev line, 
using theoretical values for the efficiency of the counter at the two 
energies. 

In the literature the chemical form of the source and the chemical 
separation procedures are frequently inadequately described when the 
detection of a very weak line is claimed. It is not therefore possibie 
on the basis of published information, to say whether or not some of the 
radiations reported could arise from fluorescence in a chemical component 
or impurity in the source, but the general inconsistency might well arise 
from a factor such as this, not usually specified and probably varying 
from experiment to experiment. In particular a line in the region of 
31 kev could be caused by the fluorescence of a Ba (Z=56) impurity, 
although this could not be excited by the L-conversion electron group. 
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The Axial Ratio of Zinc, and of the Eta and Epsilon Brasses 


By F. R. N. NaBarro 
University of the Witwatersrand, Johannesburg 


[Received February 2, 1957] 


Some years ago, Nabarro and Varley (1952) discussed the equilibrium of 
a hexagonal metallic lattice, taking into account only the energy of central 
forces between each ion and its twelve nearest neighbours, and the Fermi 
energy of a gas of two free electrons per atom. They calculated the 
energy as a fraction of a, the distance between nearest neighbours in 
the basal plane, and b, the distance between nearest neighbours in neigh- 
bouring hexagonal planes. It was shown that any reasonable central 
force always leads, on this model, to the existence of a mechanically 
stable close packed structure with a=b. Some central forces lead to a 
second mechanically stable structure, which is always prolate, with a<b. 
A detailed discussion was given of a particular case in which the prolate 
structure had b=1-2a, and its binding energy exceeded that of the 
close packed structure by 0-23 kcal/mol. 
It was pointed out that the axial ratio of zinc, and of the 7 brasses, 
differed from the value corresponding to close packing by an amount °* 
much greater than that usually observed to arise from Brillouin zone : 
effects, and it was suggested that these alloys, including pure zinc, were 
prolate structures which could not be considered as being formed by 
perturbation of the close packed structure. This view is supported by 
the fact that the increase in mean atomic volume in’ passing from the € to 
the 7 phase much exceeds the increases in passing from « to f, from 
B to y, or from y to « (Owen and Pickup 1933). On the other hand, 
the « brasses, which form when zine is alloyed with 2-3°%, of copper 
(2-7% at 380°C), may be regarded as perturbed close packed structures. 
Since the ionic radii of zinc and copper are not very different, it seems 
reasonable to represent the addition of copper to zine in our model by a 
reduction in the number of electrons per atom, with no change in the 
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central force. With two électrons per atom, the central force considered 
by Nabarro and Varley leads to a prolate structure. In the absence of 
an electronic contribution, the only mechanically stable structure is close 
packed. We therefore expect that a reduction in the number of electrons 
per atom will cause the energy of the prolate structure to increase relative 
to that of the close packed structure. Below a certain electron/atom 
ratio the prolate structure becomes thermodynamically unstable, while 
at a still lower electron/atom ratio the prolate structure is not even 
mechanically stable. 

Computations using the central force law assumed by Nabarro and 
Varley have been completed by Mr. G. Hawkes and Mr. and Mrs. R. 
Innes, and the resulting mechanically stable structures are shown in the 
following table, which uses the arbitrary units of the earlier paper. 


Close packed structure Prolate structure 
Electrons 
per atom eas Binding 5 b Binding 
energy energy 

2-00 1-0826 0-243178 1-0000 1-2000 0-245245 
F-99 1-058 0-248903 0-9988 1-1980 0-250791 
1-98 1-029 0-254941 0-9974 1-1958 0-256365 
1-97 1-0140 0-261300 0-9965 1-1932 0-261962 
1-96 1-0126 0-267718 0-9953 1-1908 0-267588 
1-95 1-0114 0-274165 0-9935 1-1896 0-273239 


The prolate structure is the more tightly bound when the electron/atom 
ratio exceeds about 1-962, and less tightly bound when this ratio is smaller. 
Thus the model predicts a sudden change of axial ratio as the concentra- 
tion of copper increases. The two phases have the same energy at a 
concentration of 3-8°%, of copper, and the transformation is characterized 
by a slight increases in a and a large decrease in b. In fact, the two- 
phase e-n region in brass extends from 2-7 to 13-2°% copper, and the 
transformation 7-<« is characterized by a slight increase in a and a large 
decrease in b. Numerical agreement is not good, largely because the 
central force assumed leads to an unreasonably large axial ratio for pure 
zine. Moreover, our model, which neglects Brillouin zone effects, can 
give no deviations from exact close packing in the « phase, gives the 
wrong sign for the small variation of a with composition in the 7 phase, 
and does not give the observed rapid variation of 6 with composition in 
the » phase. 

It seems clear that the changes of lattice parameters within either the 
< or the 7 phase can be accounted for only by taking into consideration 
the effects of the Brillouin zone, and probably other factors. On the 
other hand, the Brillouin zone effects will not account for the transition 
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from ¢ to a highly prolate structure, although they may modify the 
behaviour of this transition considerably. While the model discussed 
here is primitive, and arbitrary in its details,.it predicts a transformation 
of the type observed, and thus illustrates a phenomenon which is likely 
to be important in the theory of the systems Zn Cu and Cd Ag. 


REFERENCES 


NapaRro, F. R. N., and Varuey, J. H. O., 1952, Proc. Camb. Phil. Soc., 48, 
316. 
Owen, E. A., and Pickup, L., 1933, Proc. Roy. Soc. A, 140, 179. 


ERRATUM 


In the February issue of the Philosophical Magazine, Vol. 2, No. 14, 
p. 245, the name of the author was incorrectly printed. It should read 
O, KRISEMENT. 
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‘A 


Edge dislocations on slip planes in a crystal of potassium bromide made visible 
; by decoration with gold. 


Two dislocation walls in’a polygonized crystal of sodium chloride. 
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(b) 


(c) 


A source of closed dislocation loops in a crystal of sodium chloride. 


Fig. 4 


A helical dislocation in a crystal of sodium ‘chloride. 


Notr.—In figs. 3 and 4, (a), (b), and (c) were taken with the focal plane of the 
microscope, above, coinciding with, and below the axis, of the system or of the 
helix respectively. 


